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X XVII. BAND ZWEITES HEFT 1959 


Die Beschleunigungsanderung 
| II. Mitteilung 
Von W. Meyer zur Capellen 


1. Einleitung. In Fortsetzung der I. Mitteilung werden im Folgenden wichtige Sonderfalle und 
Anwendungen auf Rader- und Kurbeltriebe behandelt. 


2. Besondere Punkte der Ebene. Wir wenden uns zunachst einigen Sonderfallen und zwar besonderen 
Punkten sowie Sonderheiten der ganzen bewegten Ebene zu. Fiir Punkte auf der Kreisungspunkt- 
kurve k,, verschwindet das zweite Glied in i, aus! I (7c), wahrend fiir Punkte auf dem Wendekreis 
K,, das erste Glied in §{,, aus I (7c) und das zweite Glied in St, aus I (7b) verschwindet. Dort ist 
also R, = db,/dt. Im Schnittpunkt von k, und K,,, d.h. im Ballschen Punkt, verschwindet auch 
R,,, mit anderen Worten: Die drei geometrischen Orter, d. h. Kreisungspunktkurve k,, Wendekreis 
K,, und Tangentialkreis Kp schneiden sich auBSer im Momentanpol P noch im Ballschen Punkt U, 
wobei tg y,, = — m/l, ist. 


4 y 4 Polbahnnormale 


Wr. 


Normalkre/s 


Tangentialkreis 


‘Wendekreis 


ae 
Polbahntangente 


Abb. 1. Zusammenhinge zwischen Wendekreis, Tangentialkreis und Normalkreis. 


Ferner ist die Tangentialkomponente R, des Ruckes auf dem Wendekreis wie bei der gerad- 
linigen Bewegung gleich der Ableitung der (Tangential-) Beschleunigung. Auch gibt es auf dem 
Wendekreis einen Punkt, fiir welchen R, noch verschwindet, also 6, ein Extremum hat. Dies ist 
der Schnittpunkt S von Wendekreis und Normalkreis. Setzt man in I (27a) oder in I (18) die Glei- 
chung r = D sin y des Wendekreises ein, so folgt fiir den entsprechenden Strahl der Winkel 

3 D/m* 
tg?s = to DIR, , (1) 


wobei in den Umformungen der Faktor 1 + 3 D/l, =1+D(c’ —#) =1+ D/R,—1= D/R, 
nach I (13) geschrieben werden konnte. 


1 T bedeutet I. Mitteilung. 
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Bei der betrachteten Konfiguration gibt es also sieben wesentliche Punkte (um es nochmals zu 
sammenzufassen): den Momentanpol P, den Wendepol W, den Beschleunigungspol J, den Ruckpo. 
Q, den Tangentialpol Wy, den Ballschen Punkt U und den Punkt S mit extremaler (Tangential- 
Beschleunigung; vgl. Abb. 1, in welcher J nicht eingetragen ist. 


3. Sonderfalle im Gesamtzustand der hewegten Ebene. a) Kinematisch. (Tabelle1.) Bemerkens 
werterweise liegen dann, wenn die Winkelbeschleunigung der gesamten, allgemein bewegten Eben 
ein Extremum hat (Fall 2), keine wesentlichen Besonderheiten vor. 


Tabelle 1} 


Merkmal or) Bemerkungen 
a eee 
| | | tg Yo = Tr Tangentialkreis 
| | o | leich Gerade, d. h. 
1 ees erie ew n[2 W x a symmetrisch c Gn tjoaenek Strahl Yo 
| ee a zum Strahl | | (Senkrechte zu {tp) 
| | | PU | s 
| ieee 3 
ees oe _3p2 iL Ve | cal : w? __| keine wesentlichen 
vik 2° R he 3g Besonderheiten? 
| D 
| Normalkreis gleich 
; | | Gerade, d. h. iden- 
= 8 
3 [ae ee ak n/2 \Wr tisch Strahl go +3 
| | (Richtung von Sip) 
1 Striche heiBen: keine Besonderheiten. 
2 der geometrische Ort fiir b, = Extremum gemaB I (27) wird eine Kurve dritter Ordnung. 
Tabelle 2 
Nr.) Merkmal | m|m* i, k,, ka T N %o U, ey 
Polbahn- Polbahn- 
normale | normale | a 
fies __|___ und Kreis)und Kreis) 2/9, = = | | dure 
a cant Goh vom yon | ay Se arts v2) .— Des ae 
| Durch- | Durch- U=W stati 
_messer 1 messer Ig nat 
| | | 
| ae Polbahn- 
| | tangente D Ro = 
Bap ane al (yee _ und Kreis| “& oA | IR. 
Aes a el vom m* | 0 | (P) 0 er Kard 
Durch- | | lag: 
| messer m | 
Sn UEC a a a a a 
| | Polbahn- Ro: 
| tangente | Rul 
: pal Leen ool enc Kcreis | | Atkg = 
6 20+7'=0 co -D vonei = aay — —_— 3Deo?} — | Um 
| Durch- | | keh: 
messer m | Kard 
| | lag 


3 


gehért ebenso zur Tafel I. 
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b) Geometrisch. Beidenin Tafel 2 angegebenen Sonderfallen handelt es sich, abgesehen von Nr. 7, 
m die wichtigsten Entartungen der Kurven k, und k,,, wovon auch bereits an anderen Stellen)2 
rebrauch gemacht wurde. Dariiber hinaus sind auch Kombinationen der einzelnen Faille méglich, 
yie aus den folgenden Beispielen zu erkennen ist. 


4. Radertriebe. 1. Beispiel: Allgemeiner Umlaufradertrieb. Beim einfachen Umlaufradertrieb 
Abb. 2) rollt der Kreis k, vom Halbmesser Ryg = R, auf dem Kreis k, vom Halbmesser R, ab, und 
ler Steg werde mit konstanter Winkelgeschwindigkeit angetrieben. Dann ist die Winkelgeschwin- 
ligkeit @ =, des rollenden Kreises, d. h. der allgemein bewegten Ebene auch konstant, es ver- 


Polbahnnormale 


SEX 


y, Tangentialkreis 


NN WwW 4 ey Extrema von by 


VM iN OE, hh 
O 


“oe 
Q 


Abb. 2. Zum Ruck am Umlaufradertrieb mit n = R,/R, = 3. 


a ear: 
i N yy 
g Vy 
\ Se 


Polbahntangente 


rmalkre/s 


chwinden ¢ und é. Ferner bleibt der Wendekreisdurchmesser D konstant, d.h.? es gilt #’’ = 0. 
omit liegt eine Kombination der Falle 1), 2) und 4) bzw. 7) vor, welche eine Reihe von Verein- 
achungen bringt. So wird y, = 7/2 und R,, =0 sowie N = — 3 D?/I,. Der Tangentialkreis fallt mit 
ler Polbahnnormalen zusammen (Abb. 2). Setzt man R, =n R,, so wird 


n 
cee yeaa 
ind : 
Pi De 
=D )=Diz 1); 
eh. 
n—2 n(n — 2) 
eS et (n — 1? ie 
Iso 


R,. = No’ und R,, = Nw? cos p 


wnabhangig von r, speziell gleich Null fiir p = 7/2 (entarteter Tangentialkreis). 
Als Besonderheit muB noch die Entartung der Kurve fiir ein Extremum von b, gemaf I (27 b) 
1ervorgehoben werden. Da é und 1/m* verschwinden, gilt jetzt (mit R,, = R,) 


«w? D? sino (— see sing] =0 


1 yg]. FuBnote 2 von S. 58 der IJ. Mitteilung. 

2 vel. W. Meyer zur Capellen, Konstruktion 8 (1956) S. 510 sowie 9 (1957) S, 344. : 

8 W. Meyer zur Capellen, Arch. Getriebetechnik 4 (1936) S. 44; W. Meyer zur Capellen, Arch. Getriebe- 
echnik 4 (1936) S. 577; W. Meyer zur Capellen, Forschung Ing.-Wes. 24 (1958) H. 6. 


6* 
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oder 
r=R,sing und o=0; 


d. h. der geometrische Ort zerfallt in den ,,Halbrollkreis“‘, d. h. in einen Kreis vom Durchmesser R, 
welcher die Polbahntangente in P berihrt, und in die Polbahntangente (insgesamt Kurve dritte: 
Ordnung). Im ersten Fall liegen die zugehérigen Kriimmungsmittelpunkte auf dem Kreis vom 
Durchmesser R,, der ebenfalls in P die Polbahntangente berihrt, und im zweiten Fall ist immer F 
der Kriimmungsmittelpunkt (Abb. 2 fiir n = 3). 

Fiir n = 00 (Orthozykloide) wird D = R, und N sowie R,,, d.h. Jt, = 0. Der eben genannte 
Kreis wird zum Wendekreis. 

2. Beispiel: Die Kardanbewegung. Fiir n = 2 (Kardanbewegung) ist D = 2 R= Ky tera 
N =— R, (Normalkreis gleich ,,Halbriickkehrkreis‘‘ und Q der Mittelpunkt von K,). Auch 146i 
sich die Behauptung iiber die Lage der Punkte, fiir welche momentan 6, ein Extremum hat, unmit- 
telbar priifen. 


Polbahntangente Q 


Normalkreis 


_, geom. Ort fiir 
Extrema von b, 


Poloahnnormale | 
Abb. 3. Zum Ruck am Kandankreispaar (n = 2). 


Unter Benutzung der auf Seite 55 der I. Mitteilung angegebenen Gleichung fiir die Ellipse, wort 
jetzt ~« =wt den Winkel des Steges bedeutet (Abb. 3), sowie der Formel fiir v wird auch 


v= S22 /@ + B?) + (a? — B) cos2.a 


und daraus 


= —_ = oes Oj == Be 

‘dt y2 (a? + b?) + (a2 — b2) cos 2.x 
Differentiiert man nochmals, so fiihrt die Bedingung db,/dt = 0 auf die in cos 2 « quadratische Glei 
chung 


a? — b? sin 24 


2 (a? + 5?) cos 2 « + (a? — b?) (cos? 2a +1) =0 


mit der Lésung 


a—b 
cos 2q = — ey, > 
: : ie b : 
(Die zweite Lésung cos 2 « = — : Be ,— | scheidet zunachst aus). Liegt der bewegte Punkt A i 


nerhalb des kleinen Kardankreises, des Rollkreises, so ist A*A = a, B*A = b, R. = = Awe 4 
1 ert il : ¢ ‘ 
3 (a + b) und p = MA = ~ (a — 6); d.h. die Bedingung lautet cos (%7 — 2 x) = p/R,; also lee 


| 


| 
| 
| 
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wie aus Abb. 3 unmittelbar zu ent- %, 
nehmen, der Punkt A mit extrema- ri 
‘ enge BA 
ler Beschleunigung tatsachlich auf 
08 


dem oben genannten Extremalkreis!. 


Liegt der bewegte Punkt A, 
auBerhalb des Kardankreises, so 


mea A,— a, B*A, = b, aber 
ik ———— 
R, = = A*B* == (a—b) und 


p = MA, = — (a + 8); also gilt 
dann cos (7 — 2a) = R,/p, somit 
liegt A, auf der Polbahntangente, 
wie oben behauptet?. Wiirde man 
hier die Strecke b, da entgegenge- 
setzt wie vorher, negativ einsetzen, so 
hatte die zweite Lésung fiir cos 2 « 
auch ihren Sinn. Der Extremal- 
wert selbst wird b, = + (a — b) w®, 
also im ersten Fall gleich F 2 pw?, 
im zweiten Fall gleich +2 Rw? 6, 
=+R/o? = + Do? =+5,, wie 


auch zeichnerisch leicht einzusehen, 


da ja (vgl. oben), 4 M, = by, wenn 
«© = 1 gezeichnet wird. 


5. Kurbeltriebe. Die folgenden 
Beispiele betreffen vier allgemeine 
und vier besondere Falle. 


1.Beispiel: Kurbelschwinge. 
Es soll bei einer Kurbelschwinge 
(Abb. 4) eine der Stellungen,in denen 
die Winkelbeschleunigung des Glie- 
des 3 und damit die Tangentialbe- 
schleunigung des Punktes B ein Ex- 
tremum hat, ermittelt werden. Die 
Beschleunigungen des Punktes B 
wurden nach den oben angegebenen 
Konstruktionen ermittelt, und dar- 
aus wurde geschatzt, wo ein Ex- 
tremum auftritt. In der Nahe dieser 
Stellung wird nun der Ruck unter 
Benutzung der in Abb. 4a und b er- 


mittelten Werte fiir die Geschwin- Pee 
digkeiten und Beschleunigungen 7 
bestimmt. Dann ist zunachst, da Abb. 4. Zum Ruck an einer Kurbelschwinge in allgemeiner Lage. 


: Besti der Geschwindigkeit b) Besti der Beschleunigu 5 
D4t-acm, gezeichnet rde und a) Bestimmung der Geschwindigkeiten, b) Bestimmung der Beschleunigungen 


~) = konstant ist, Ji, = —v w,also Ry + acm, und entgegengesetzt gerichtet zub 4. Ferner gilt fiir 
den Ruck 3iz des Punktes B zunachst Gleichung I (2b), und die Normalkomponente betragt 


Ripa = 3 Upa& mit & als Winkelbeschleunigung der Koppel, wobei ¢, = bp 4/AB = tg v, 


1 Dadurch sind die an anderer Stelle gefundenen Ergebnisse in einem breiteren Zusammenhang geriickt; 
vgl. W. Meyer zur Capellen, Werkstatt und Betrieb 91 (1958) S. 723; ferner Maschinenbautechnik 7 (1958) S. 404, 
und Vortrag auf der VDI/AWF Fachtagung Getriebetechnik in Konstanz: Uber raumliche Kurbelschleifen 
und ihre gleichwertigen ebenen Getriebe, VDI Getriebeheft Bd. 29 (1958) S.91. In der hier betrachteten Stellung 
wird bei der elliptischen Schleife mit dem Steg Null (gleichwertig dem Kreuzgelenk) das Ubersetzungsverhalt- 
nisi = 1. 

2 Wendet man auf die Bahn des Punktes A den Satz von der doppelten Erzeugung der zyklischen Kurven 
an, so liegt der Punkt A im zweiten Getriebe wieder auf dem ,,Halbrollkreis“; vgl. FuBnote 3 von Seite 75; 
ferner W. Meyer zur Capellen, Konstruktion 8 (1956) S. 268. 
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Abb. 4,c) Bestimmung der Ruckkomponente R,pA=3 pa ° & 4) Bestimmung der 
Ruckkomponente Rnp = 3p. 6s, e) Bestimmung des Ruckes nach dem Satz von 
Euler, £) die Koppelebene in ihrer Gesamtkonfiguration. 


ist (Abb. 4c). Somit kann, wie dort gezeigt, sofort Ripa abgegriffen werden. In gleicher Weisel 
kann R,z = 3 vg €, ermittelt werden (Abb. 4d). Man tragt nun an die Spitze des in B ange-| 
tragenen Vektors i (Abb. 4e) den Vektor Rt, p 4 parallel zu BA an, ebenso von B aus den Vektor 


i el i ot 
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tp auf. Die Senkrechte zu B,B im Endpunkt der Spitze von Si,,p und die Senkrechte zu AB im 
Endpunkt der Spitze von i, p 4 treffen sich in X und liefern damit die Tangentialkomponenten Sit, z 4 


Seer ° . ° 
und i, 7 von Itz = BX. Hierinist nun i,g = rpg é; —rg ow’ =rgé, — b, pW, wobei der letztere Betrag 
Yo = bp tg B, mit tg 0, = vg/rg leicht zeichnerisch gewonnen werden kénnte. 
Abb. 4f zeigt die Kurbelschwinge innerhalb der Gesamtkonfiguration. 
Wenn nun ¢, ein Extremum, also &, gleich Null sein soll, so ist R,, gleich dem Wert y,. Tragt 


_man dann den Quotienten 


ao et 
TBS Ya 
tiber « auf (Abb. 5a), so ergeben die Schnittpunkte der Kurve y mit der Parallelen Y = — 1 zur 


o-Achse diejenigen Kurbelwinkel a, fiir welche ¢, ein Extremum hat; vgl. a. Abb. 5b hinsichtlich des 
Verlaufes von ¢, in Funktion des Kurbelwinkels «. Wohlgemerkt benétigt man die Konstruktion 
des Ruckes nur in wenigen Stellungen in der Nahe des Extremums. Das Ergebnis der ahnlichen 
Untersuchung fiir die Koppel zeigen Abb. 5c und d. 

2. Beispiel: Schubkurbel. Bei einem geschrankten Schubkurbelbetrieb ist in Abb. 6a 
nach dem gleichen Verfahren der fiir den Ruck erweiterte Satz von Euler benutzt worden, und zwar 


ee ee a 8) 
~ 


| 
| 
| 
| 
| 
| 
— 
a ee 
| 
| 
| 
: 
| 
| 


Abb. 5. Zur Ermittlung des Kurbelwinkels o des Antriebes der Kurbelschwinge nach Bild 4 fiir die Stellung maximaler Beschleunigung. — 


Ri B 
a) der Quotient y = Pn oe (x) fir die Abtriebsbewegung; b) die Bewegungsgesetze des Abtriebes; c) der Quotient « = ae = h(a) 
3 Ta 
fiir die Koppelbewegung; d) die Bewegungsgesetze der Koppel. 


in der Stellung, in welcher die Beschleunigung b des geradgefiihrten Punktes B ein Extremum hat. 
Die Stellung konnte auf Grund der Mafbeziehungen A = ale = Kurbellange/Koppellange und ¢ = 
e/e = Schrankung/Koppellange einer anderen Arbeit entnommen werden,1 und die zeichnerische 


1 W. Meyer zur Capellen u. Mitarbeiter, Bewegungsverhiltnisse an der geschrankten Schubkurbel, For- 
schungsbericht Nr. 449 Opladen 1957. 


80 W. Meyer zur Capellen: Die Beschleunigungsanderung. II Ingenieur-Archiy | 


Lésung bestatigt die dort gefundenen Ergebnisse. Gegentiber Beispiel 1 vereinfacht sich das Ver- | 
fahren ein wenig, da die Normalkomponente R,, des Ruckes fiir B immer verschwindet. 

Fiir den Bewegungszustand der Koppelebene ergeben sich in dieser Stellung folgende Besonder- | 
heiten: B ist der Ballsche Punkt U, d. h. B liegt auf dem Schnittpunkt von Wendekreis und Tangen- | 
tialkreis. Da R, = db,/dt — v3/o? = R verschwindet, weil @ = ©0 ist und db,/dt = db/dt gleich Null | 
ist, liegt B auch auf dem Normalkreis, d. h. aber Ruckpol und Ballscher Punkt fallen zusammen, wie | 
in Abb. 6b angegeben ist. 


| 
| 
| 
| 


Abb. 6. Der Ruck einer Schubkurbel in der Stellung fiir die maximale Beschleuni d i 
a) die Bestimmung der Geschwindigkeit, der Beschleunigung und des Rucker det aE a ner 
Koppelebene in ihrer Gesamtkonfiguration. : 
3. Beispiel: Totlage. Wenn eine Kurbelschwinge sich in der Totlage befindet, und die Koppel 
senkrecht zur Schwinge steht (Abb. 7), so hat, wie sich aus dem Satz von Euler leicht ergibt}, die Win- 


1 Vgl. FuBnote 3 von Seite 75, 
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kelgeschwindigkeit ~@, der Koppelebene ein Extremum, d. h. es ist €é, = € = 0. Dann zeigt die An- 
wendung des Satzes von Euler fiir den Ruck, daB bp = b,, und damit ¢, ein Extremum hat. Dies folgt 
auch aus allgemeinen Erérterungen: ¢ =0 und m= g (= 0) bedeutet die Kombination der Falle 
1 und 7 aus den Tafeln 1 und 2. Der Tangentialkreis geht in die Polbahnnormale iiber, Riy= Rip 


Abb. 7. Die Kurbelschwinge in der Totlage (x ABB, = m/2). 


_verschwindet, und da R,, = db,/dt — v?/o? = db,/dt ist, weil ja v im Augenblick verschwindet, so 
mu db,/dt auch verschwinden, also b, , = b, ein Extremum haben, und zwar, wie nach dem Satz von 
Euler leicht einzusehen ist, von der Gréfe by (1 + a/c). Wird B,B unendlich groB, so haben wir den 
zentrischen Schubkurbeltrieb in der auBeren Totlage. 


Steht bei einer Kurbelschwinge in der inneren Totlage die Schwinge senkrecht zur Koppel, so 
_ergeben sich die gleichen Dinge, nur hat das Extremum den Wert — b, (1 — a/c). 


| 4. Beispiel: Doppelkurbel. Dazwischen den Winkelgeschwindigkeiten der Glieder 1, 3 und 
64 = Odie Beziehunga.s) = @1) + 3, besteht und w,) = konstant ist, so muB G4) = 5, d.h. Eg) = &5), 
also in der betrachteten Stellung auch é,, gleich Null sein. Halt man also jetzt das kleine Glied a 
_ fest, so mu bei der entstehenden Doppelkurbel (Abb. 8), bei welcher die Abtriebskurbel sich auf dem 
| Steg befindet und die Koppel senkrecht zu diesem steht, die Winkelbeschleunigung der Koppelebene 
und die Winkelgeschwindigkeit des Abtriebgliedes! ein Extremum haben. Dies folgt auch unmittel- 
bar: Da A, mit dem Momentanpol zusammenfillt, ist vy = €@1)9 = 49, d. h. © = Wy = Wop. Fer- 
ner liegt A als Punkt des Gleichenkreises K, auf der Polbahntangente, d. h. es ist A = G, und der 
Kreis tiber PA als Durchmesser ist der Gleichenkreis, d. h. PA = (—) g. Dieser Kreis geht auch 
durch B, also hat B keine Tangentialbeschleunigung, es ist ¢, = 0. Da aber A und B auch auf der 
| Kreisungspunktkurve k, liegen miissen, entartet diese in den Gleichenkreis und die Polbahntangente, 
mes ist g — m = (—) PA, 1 = oo, d.h. 1, = — D: Wir haben die Kombination der Faille 1, 6 und 7 
aus Tafel 1 und 2. Hiernach wird T = ne =—m= PA und gy =2/2, d. h. der Gleichenkreis ist 
0 

auch der Tangentialkreis. Nun gilt fiir den Ruck des Punktes A als Endpunkt der Kurbel doch tan- 
gential R, 4=— PA w’, da €, = 0, andrerseits aber als Punkt der Koppelebene R,4 = PA (&,— 3), 
) 


und weil beide Werte einander gleich sein miissen, folgt mit @ =w, =, auch ¢, = 0. Nach 
dem Satz von Euler fiir den Ruck ist aber R, 4 = R,, + 39 & PA, und muB dieser Wert, da A auf 
dem Tangentialkreis liegt, gleich Null sein. Nun ist gema® Tafel 2, Nr. 6, doch R,, = + 3 D (w,)°, 
d.h.es muB ¢ = ¢, = —w* D/m werden, wie ja auch aus der Definition des Gleichenkreises (m = g) 
oder aus dem Satz von Euler fiir die Beschleunigungen, angewendet auf Punkt B, folgt. Der Durch- 
-messer des Normalkreises ergibt sich beilaufig aus Nr. 6 der Tabelle 2 zu N = 3 D (Abb. 8), wobei 
der Wendekreis mit der Hartmannschen Konstruktion zeichnerisch ermittelt wurde, allerdings 
aus der Euler-Savaryschen Formel fiir das Punktepaar B, B, auch rechnerisch bestimmt werden 


_konnte. 


1 Vgl. FuBnote 3 von Seite 75, zweite Arbeit. 


| 
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5. Beispiel: Parallellage. Die Kurbelschwinge befinde sich in der Parallellage, d. h. P sei 
unendlich fern (Abb. 9). Dann ist w, = 0 und, wie an anderer Stelle gezeigt wurde}, aber auch 


Ke=kersk 
Abb. 8. Die Doppelkurbel in der Steglage Abb, 9. Die Kurbelschwinge in der Parallellage 
(x ABB, = n/2, Rg 2 a/3). (x ABB, = 1/2,0 = P). 


leicht mit dem Satz von Euler fiir die Geschwindigkeiten nachzuweisen ist, ¢, = 0. Damit gilt | 
einerseits 
Ie = Vat Ries t Riza 
wobei 
R, = (—) ao, Ripa 3 O,€,¢— 0 
und 


Riza = ¢(&,— 02) e#), da id WE 
Andrerseits gilt ,=R,p+ Vip, wobei R, zg =3,€,B,B=0, da e, =0, und beilaufig R,p = | 
B,B (é;—@?) ist. Es muB also R,p4=R,,g, d.h. é, = 0(und R,z = Ry) sein. Ebenso wird, da hier 


R,= R,p und 
aa elas an 
1 B,B b 
gilt, 
: [CaP 
a ea 3 b2 ) 1 


gefunden. Im iibrigen weist die in Funktion des Kurbel- 
winkels « aufgetragene Kurve des Koppelwinkels y in der 
betrachteten Stellung ein Extremum auf (da = @, = 0), 
ebenso die Kurve der Winkelbeschleunigung y, = €,, bzw. hat 
die Kurve fiir w, hier einen Wendepunkt (da @, = &, = 0). 

6. Beispiel: Doppelkurbel. Betrachtet man die Re- 
lativbewegung von Glied 3 gegen Glied J, d.h. halt man | 
das kleine Glied fest, so erhalt man eine Doppelkurbel in | 
der Stellung, in welcher Koppel und Steg einander parallel 
sind und die Abtriebskurbel senkrecht zum Steg steht (Abb. | 
10). Aus ahnlichen Uberlegungen wie beim Ubergang von | 
Beispiel 3 zu Beispiel 4 folgt dann, daB «, = 0 und e, = 0 
ist.2 Mit Hilfe der Satze vom Ruck 1aBt sich die letztere 
Besonderheit auch direkt nachweisen. So muB® fiir den | 


Abb. 10. Die Doppelkurbel in der Parallellage 


(x ABB, = 7/2,0 = P). Punkt B die auf S. 54 der I. Mitteilung angegebene Winkel- 
beziehung tg (x + 7/2) = 3 tg y erfiillt sein. 
Im iibrigen steht nach der Bobillierschen Konstruktion der Polbahntangente diese senkrecht 


zu PA, d.h. wir haben die Kombination der Falle 1, 6 und 7 aus Tafel 1 und 2, wobei noch] = — PA 
und 1, = — PA, wird. Da Koppel und Steg parallel sind, also der Pol P,, unendlich fern liegt, ist 


1 Vg). FuBnote 3 von Seite 75, zweite Arbeit, und FuBnote 2 von Seite 75. 
2 Vgl. FuBnote 2 von Seite 77, zweite Arbeit. 
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nach dem Satz iiber die Relativpole @,, = 0, d. h. w, = wg, woraus sofort mit &=€, = 0(R,p = Rina 
= 0) und é, = é, = 0 der Ruck des Punktes B zu Ry = —w’ ByB = — ow? b und damit 


somali f)- fa 


mit ¢c = ABundd = A,B, gefunden werden. 
7. Beispiel: Kurbelschleife. Bei der zentrischen schwingenden Kurbelschleife nach Abb. 11 
ist der Ruck fiir die Schleife nach den vorhergehenden Konstruktionen (vgl. Abb. 3 der I. Mit- 


v, 


Abb. 11. Der Ruck einer zentrischen schwingenden Kurbelschleife in ihrer Stellung fiir maximale Abtriebsbeschleu- 
nigung. — a) Bestimmung der Geschwindigkeiten, b) Bestimmung der Beschleunigungen. 


teilung) ermittelt worden und zwar in derjenigen Stellung, in welcher die Winkelbeschleunigung 
& = & = €, = &; ein Extremum hat, also é = 0 wird, wobei der Antriebswinkel « sich aus einer qua- 
dratischen Gleichung! ergibt. Nach der iiblichen Geschwindigkeits- und Beschleunigungskon- 
struktion (Abb. 11a und b), bei welcher 6, — 6, — b, =b,¢ + b,¢, und b,¢ = v3/ByA — BA te? 9, 
6, =2,,,tg 0 geschrieben wurde, bestimmt man den Ruck Ji, = Ni,,, = — V4 w® + acm gezeichnet, 
ebenso die Komponente R,; des Fiihrungsruckes ip, d.h. R,- = (—) 3 vf 7 = 3 vf tg v, wobei 
& = b,,/B,A = tg v gesetzt ist (Abb. llc) und dann die Komponenten des Zusatzruckes gemaB I 
(33c), d. h. in Abb. 11d die Komponente R,,,, =3 wsb, = 3b, tg & senkrecht zur Fiihrung gemaB 
Reusi = 30% X b,, wobei tg # = v,/B,A ist, in Abb. 11e schlieBlich R,,,. = 3 Ef v, = 3 v, tg v, eben- 
falls senkrecht b, gem4B ft,,,,9 = 3 €, X V, und in Abb. 11f ferner 9,,,,; = (—) 3]w7 v, =(—) 3b, tg? 
entgegengesetzt ),, daraus 3It,,, als geometrische Summe (Abb. 11g). Manschreibe I (33 b) in der Form 
Ty =— Np — Mavs = Mies = ote 
— > > 
und trage an A (Abb. 11h) nacheinander an: 41 = ty, = fi, 12 = — hi, ,, 23 = — R,,,. Dann 
—> —> 
trifft die Senkrechte durch 3 zur Fihrung diese in 4, und es ist 4 4 = ff, sowie 43 = fi, ;. Da nach 


1 W. Meyer zur Capellen, Werkstatt und Betrieb 89 (1956) S. 581. 
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— 

I (26b) Jiyp = e x BA — w# vs ist und é; = 0 sein soll, mub also 43 = — OF De, d. h. der Betrag 
gleich v; tg? # sein, wie in Abb. 11i konstruiert und wie sich durch Vergleich beider Strecken zeigt. 
Fs sei ausdriicklich betont, daB es sich um ein Erlauterungsbeispiel handelt, da ja die Lésung der er- 
wahnten quadratischen Gleichung schneller zum Ziel fiihrt. 


Abb. 11 c) Bestimmung der Ruckkomponente Rif =30 ef» d) Bestimmung der Ruckkomponente R 
stimmung der Ruckkomponente R,,,.. = 3 v, eps f) f 
vektor {it 


3 zus1 = 3b, we, e) Be- 

est!mmung der Ruckkomponente R,,,,, = 3v,@?, g) der Ruck- 

Op h) die Bestimmung des Ruckes nach dem erweiterten Satz von Euler fiir die Relativbewegung, i) Be- 
stimmung der Ruckkomponente Rif = 15 oe. 


8. Beispiel: Elliptische Schleife. Ebenso wie im vorstehenden 7. Beispiel bei der zentri- 
schen schwingenden Kurbelschleife (Abb. 11) kann auch bei einer zentrischen, schwingenden ellip- 
tischen Schleife+ die Stellung des Antriebes gefunden werden, bei der die Winkelbeschleunigung 
des Abtriebes ¢ = €, = é, ein Maximum aufweist, ¢ = é, = és also demnach Null sein muB, Bild 12. 
Die Durchfiihrung der Konstruktion, also die Ermittlung der Geschwindigkeiten, der Beschleuni- 
gungen und der Komponenten des Rucks erfolgt in gleicher Weise wie fiir die Kurbelschleife nach 
dem 7. Beispiel, wahrend diesmal Geschwindigkeit, Beschleunigung und Ruck des absolut bewegten 
Punktes 4 nach Abb.4 der I. Mitteilung (vgl. dort Ziff. 3b) in gewohnter Weise gefunden wurden. 


1 Vgl. FuBnote 1 von Seite 77, erste und dritte Arbeit, 
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Abb, 12. Der Ruck einer zentrischen schwingenden elliptischen Schleife in ihrer Stellung fiir maximale Ab- 
triebsbeschleunigung. a) Geschwindigkeit, Beschleunigung und Ruck des Antriebes (durch ein Kardan- 
kreispaar verwirklicht), b) Bestimmung der Geschwindigkeiten, c) Bestimmung der Beschleunigungen. 
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Als eine der Méglichkeiten wird die Bahn des Punktes A, des Antriebes, durch das Kardankreispaar 
verwirklicht (Abb. 12a). 

In Abb. 12 b wurden die Geschwindigkeiten, in Abb. 12 die Beschleunigungen mit b, = 20,.. X wy 
ermittelt. Abb. 12d bis f zeigt die Komponenten des Ruckes t,,,, wahrend in Abb. 12g und h die 


Abb.12.d) Bestimmung der Ruckkomponente Rus 1= 36, © fs e) Bestimmung der Ruckkomponente R, 2 
a f 4 2Us 
= 3 uv, e- f) Bees der Ruckkomponente Riss = 30, w?, g) Bestimmung der Ruckkomponente 

nf 3 vp eps h) Bestimmung der Ruckkomponente arf = 3 oF oF » i) die Bestimmuug des Ruckes 


nach dem erweiterten Satz von Euler 


ur die Relativbewegung. 


Komponenten des Ruckes {ty dargestellt sind. Im Besonderen wird hier St, ; = (— : 

da das Getriebe in der Stellung éf = Null gezeichnet ist. uf \ Me eae 
Um diese Stellung in Anlehnung an das 2. Beispiel in Ziff. 4, also den zugehérigen Winkel « fiir 

das Kardankreispaar, zu finden, wird ahnlich wie dort in der Nahe des erwarteten Maximums fir 
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é der Wert y = y,/y,mity, = Ry =r; (€f —@#) und y, = rpw? (rp = Strecke AB,) iiber dem An- 
triebswinkel « aufgetragen. Der Schnittpunkt der Kurve y = ¥1/Yo, mit einer Parallelen Y = — 1 
zur Abszisse liefert, wie leicht einzusehen ist, denjenigen Winkel «, bei dem é; = Null, also die Win- 
kelbeschleunigung des Abtriebes ¢, ein Extremum aufweist (Abb. 13). Die angegebene Konstruktion 
des Ruckes zur Auffindung der Stellung fiir maximale Abtriebswinkelbeschleunigung ist nur in der 
Nahe des erwarteten Maximums fiir ¢ notwendig, zumal y = R, five wo} sehr rasch sein Vorzeichen 
wechselt, der gesuchte Winkel « also sehr gut einzugrenzen ist. In Abb. 12i wurde fiir ein Minimum 
von é die Ruckkonstruktion an der elliptischen Schleife durchgefiihrt, wie auch der zugehérige 
Winkel a in Abb. 13 angegeben wurde. Diese geometrische Konstruktion fihrt immer und schneller 
ans Ziel als analytische Verfahren, welche hier in diesem Sonderfall der elliptischen zentrischen 
Schleife fiir die Bedingung é = 0 auf eine Gleichung 10-ten Grades in cos « fiir den entsprechenden 
Antriebswinkel « fihrt. 

6. Zusammenfassung. Es wurden fiir die Beschleunigungsinderung oder den Ruck St = db/dt unter 
Beschrankung auf die ebene Bewegung bei der Zweikurvenfiihrung (Satz von Euler) und bei der 
Relativbewegung Formeln und Satze aufgestellt, ferner wurde der Zustand der bewegten Ebene 
hinsichtlich des Ruckes in Erweiterung des Beschleunigungszustandes untersucht und auf weitere 


Abb. 13. Die Bewegungsgesetze des Abtriebes der elliptischen Schleife 
nach Abb, 12. 


geometrische Orter hingewiesen. Die Beispiele dienten zur Erlauterung, vor allem aber zur Ermitt- 
lung der Extrema der Beschleunigungen in Getrieben. Da die Zeichenarbeit zur Ermittlung der Be- 
schleunigungen selbst schon gréfer ist als die zur Geschwindigkeitsbestimmung, nimmt es nicht 
wunder, daB diese Arbeit bei dem Ruck als der zweiten Ableitung der Geschwindigkeit noch wachst 
— ,,Differentiieren rauht auf‘; aber die erforderlichen Konstruktionen benétigt man nur in weni- 
gen Stellungen, und sie fiihren schnell zum Ziel.? 


(Eingegangen am 21. April 1958.) 
Anschrift des Verfasssers: Professor Dr.-Ing. Walther Meyer zur Capellen, Aachen, Technische Hochschule. 


Say ee : 4 : : ie 
1 Bei der Abfassung der Arbeit, insbesondere bei Herstellung der Abbildungen unterstiitzte mich se 
Herr oat ing. K. A. richer und es sei ihm an dieser Stelle fiir seine Miihen herzlichst gedankt. 
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Die Lenkstabilitat des luftbereiften Kraftwagens 
gegentiber kleinen Stérungen 


Von G. Mitterlehner 


1. Einleitung. Die Stabilitat des fahrenden Kraftwagens gegen die von der Fahrbahn her und 
durch Luftkrafte bewirkten Stérungen ist ein selten behandeltes Problem. Wie in der gesamten 


Dynamik des Kraftwagens ist es auch in diesem Fall schwierig, zu allgemein giiltigen Lésungen zu — 


gelangen. Der Grund dafiir ist, daB das Kraftfahrzeug mit seinen zahlreichen Massen, Federn und 
Dampfungselementen ein System von vielen Freiheitsgraden darstellt, das durch vereinfachende 
Voraussetzungen auf einige wenige Freiheitsgrade reduziert werden muf, bevor eine Berechnung 
gelingen kann. Durch verschiedene Voraussetzungen entstehen verschiedene Modelle des Kraft- 
wagens, von denen jedes aus der Vielfalt der am Kraftwagen beobachteten Erscheinungen nur 
wenige wiedergeben kann. Trotzdem wurden damit brauchbare Ergebnisse erzielt, die die Ent- 
wicklung beeinfluBt haben. 

So hat S. de Lavaud! in mehreren gedankenreichen Arbeiten die Wirkung der gyroskopischen 


Kopplung zwischen dem ,,Flattern‘‘ der Vorderrader um ihre Lenkzapfen (Abb. 1) und dem ,,Tram- _ 


peln“ der gemeinsamen Vorderachse um eine Achse in der Fahrtrichtung aufgezeigt und der Einzel- 
aufhangung der Vorderrader den Weg gebahnt. 

rabrmciaiang H. Fromm?  beriicksichtigte auferdem das 

,ochieben“ der Vorderachse normal zur Fahrt- 
richtung, das infolge des Flatterns der Vorderrader 


Wegfall der starren Vorderachse wesentlich ver- 
bessert worden, doch tritt bei sehr hohen Geschwin- 
digkeiten weiterhin Flattern der Vorderrdder auf. 

Die Eigenschaften des Luftreifens wurden 


T. E. Schunck* die Stabilitat eines Systems, das 
neben dem ,,Schieben“ des gesamten Fahrzeugs 
das ,,Gieren“’ um die Hochachse des Fahrzeuges 
beriicksichtigt. Dabei wurden auch die Luftkraf- 
te in die Rechnung einbezogen. Die Bewegung 
Abb. 1. Spreizung und Nachlauf der Lenkzapfen der Vorderrdader. des Kraftwagens auf einer Kreisbahn wurde von 

denselben Verfassern und danach ausfihrlicher 
von Paslay und Slibar® behandelt, die auch den zusatzlichen Leistungsbedarf beim Durchfahren 
einer Kurve berechneten. Mit der Stabilitat eines Fahrzeugs mit in der Geradeaus-Stellung blockier- 
ter starrer Lenkung gegeniiber seitlichen Windsté8en hat sich vor allem W. Kamm‘ befaBt. 


7% 


2. Die Wirkung einer Unebenheit der Fahrbahn. Bei Uberfahren eines kleinen Einzelhindernisses 
der Fahrbahn durch ein Vorderrad des Kraftwagens entsteht im allgemeinen ein Moment der StoB- 
kraft, das die durch die Spurstange verbundenen Vorderrader um ihre Lenkzapfen zu verdrehen 
sucht. Dieselbe StoBkraft auf den Lenkzapfen iibertragen hat Momente um alle drei Tragheitshaupt- 
achsen des Fahrzeuges, von denen je eine mit der Langs-, Quer- und Hochachse zusammenfallen 
soll. Der Einfluf des Moments um die Langsachse wird durch die Federung des Vorderrades ver- 
kleinert, und der Einflu8 des Moments um die Querachse ist von vornherein unbeachtlich. Das 
Moment der StoBkraft um die Hochachse des Fahrzeugs verhalt sich zu dem um den Lenkzapfen 


* Sensaud de Lavaud, Les vitesses critiques d’une voiture automobile und ,,Causerie scientifique“, beide 
im Selbstverlag des Verfassers Paris 1927. 

* Becker, Fromm, Maruhn, Schwingungen im Automobillenkungen, Berlin 1931. 

° L. Huber, Die Fahrtrichtungstabilitat des schnellfahrenden Kraftwagens. Deutsche Kraftfahrtforschung, 
Heft 44, Berlin 1940. 

* P. Riekert und T. E. Schunck, Ing.-Arch. 11 (1940), S. 210. 

° P. R.Paslay und A. Slibar, Ing.-Arch. 24, (1956) S. 412. 

6 W. Kamm, Automobiltechn. Z. 56 (1954) S. 117. 


Lenkzapten entsteht. Die Fahreigenschaften sind durch den — 


von L. Huber*® untersucht. Unter Verwendung — 
seiner Ergebnisse berechneten P. Riekert und | 
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les Vorderrades wie die jeweiligen Kraftarme und ist infolge der Spreizung (Abb. 1) des Lenkzapfens 
mindestens 20 mal so grofs wie das letztere. Die Drehmasse des Fahrzeugs um seine Hochachse ist 
dagegen etwa 1000 mal so groB wie die des Vorderrades um den Lenkzapfen, so da® nur das Ein- 
schlagen der durch die Spurstange verbundenen Vorderrader mit seinen Folgen und nicht der redu- 
fierte exzentrische Sto um die Hochachse gefahrlich werden kann. 

Das Fahrzeug wird als lineares System mit drei Freiheitsgraden-Flatterschwingung der 
Vorderrader, seitliche Schiebung und Gierschwingung um die Hochachse des gesamten Fahrzeugs — 
angenahert dargestellt. Das Ziel ist das Aufsuchen der Stabilitatsgrenzen und eine Berechnung 
des sog. Nachlaufes der Vorderrader (n in Abb. 1), welcher ausreicht, um ein gegebenes Fahrzeug 
yegen Stérungen zu stabilisieren, die von kleinen Unebenheiten der Fahrbahn ausgehen. 

Das Wanken des Aufbaus um seine Langsachse soll als vierter Freiheitsgrad des Systems in 
einer spateren Untersuchung mitbehandelt werden. 


3. Die Stellung der Lenkzapfen der Vorderrader und der Schriglauf der Luftreifen. Der Lenkzapfen 
eines Vorderrades steht weder parallel zur Radebene, noch liegt er in der lotrechten Meridianebene 
des Rades (Abb. 1). Die erste Eigenschaft, die sk 
Spreizung, halt das Moment der oben erwahnten [x 
StoBkraft um den Lenkzapfen klein, kann es aber 
wegen der Reifenbreite nicht ganz zum Verschwin- \ 
den bringen. Die zweite Eigenschaft, der Nach- 
lauf — der Schnittpunkt des lotrechten Raddurch- 
messers mit der Fahrbahn folgt in der Fahrtrich- 
tung dem Schnittpunkt des Lenkzapfens mit der 
Fahrbahn nach — bewirkt ein Moment der beim 
Kinschlagen des Rades auftretenden Seitenkraft S 
(Abb. 2) um den Lenkzapfen, das das Rad in 
seine Geradeaus-Stellung zuriickzufiihren sucht. y=158¥2 
Nachlauf und Spreizung werden gewohnlich als 
Winkel angegeben, im folgenden wird jedoch als 
Nachlauf n der Kraftarm der Seitenkraft S um / 
den Lenkzapfen bezeichnet. il 

Nach ausgedehnten Versuchen kam L. Huber 
zu dem Ergebnis, daB ein luftbereiftes Rad eine xg / eee ee 
Seitenkraft nur dann von der Fahrbahn auf das / V 
Fahrzeug iibertragen kann, wenn die Radebene y, 
mit der fortschreitenden Geschwindigkeit des Rades 
den sogenannten Schraglaufwinkel « einschlieft. / g 
Abb. 2 zeigt den Zusammenhang zwischen Sei- 
tenkraft und Schraglaufwinkel bei verschiedenen 
Radlasten. Die Fahrgeschwindigkeit hat auf diese 
Ergebnisse nur wenig EinfluB, doch scheint bei 
kleinen Schraglaufwinkeln die Seitensteifigkeit des HEE RA RE 
luftbereiften Rades mit wachsender Geschwin- ita neticdetiack<r. Hien, > 
digkeit abzunehmen. 


750 


100 


a 
0 on 70° 157 


4. Voraussetzungen. Im folgenden wird die Seitenkraft der Radlast N und dem Schraglaufwinkel 

proportional angesetzt: 

S =k, IN j= hoe 
Die entsprechende Gerade fiir N = 200 kp in Abb. 2 zeigt, dafs der daraus entstehende Fehler bis 
mu « = 5° tragbar erscheint. Damit sind die Grenzen gegeben, innerhalb der das Problem linear 
behandelt werden kann, wenn auch bereits in diesem Bereich eine bessere Anpassung an die Ver- 
suchsergebnisse wiinschenswert ware. 

Wegen der Beschrankung der Schraglaufwinkel kénnen die Spuren der Rader auf der Fahrbahn 
nur schwach gekriimmt sein, so da die Spuren des rechten und des linken Vorderrades einander 
gleich angenommen werden kénnen. Dasselbe gilt fiir die Spuren des rechten und des linken Hin- 
terrades. 

Wegen der Kleinheit der Winkel kénnen die Komponenten normal zur Fahrtrichtung der in 
den Radebenen wirkenden Krafte, d.i. der Schubkrafte der Treibrider und der Rollwiderstande 
aller Rader vernachlassigt werden. Durch die Drehung der Rader um ihre Achsen entstehen beim 

7 
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Gieren des Fahrzeugs um seine Hochachse und beim Einschlagen der Vorderrader Kreiselmomente | 
um die Langsachse des Fahrzeugs. Sie tragen zum Trampeln der Achsen, dem sogenannten »Shim- | 
my“ bei, das beim Ubergang von Hochdruck- auf Niederdruckreifen zusammen mit dem Flattern| 
der Vorderrider erstmals starker hervortrat und von S. deLavaud und H. Fromm ausfihrlich behan- ) 
delt wurde. Ihr Einflu8 auf die Radlastverteilung ist wegen der Kleinheit der Drehmassen der: 
Rader und ihrer Drehgeschwindigkeiten um Hochachse baw. Lenkzapfen bei nicht zu hohen Fahr- 
geschwindigkeiten gering und kann in dem gewahlten System nicht aufscheinen. Dies ist einer der} 


Griinde, die fiir eine Erginzung des Systems durch das Wanken als viertem Freiheitsgrad sprechen. | 


Abb. 3. GrundriB des Fahrzeugs. 


Die Radlast wurde zunichst fiir alle Rader einem Viertel des Fahrzeugsgewichts gleichgesetzi | 
und danach der EinfluB des Luftwiderstandes auf die Verteilung der Radlasten auf Vorder- und 
Hinterrader beriicksichtigt. Eine andere Aufteilung der statischen Radlasten bringt keine grund- 
satzlichen Schwierigkeiten, macht aber die Rechnung durch das Auftreten zahlreicher weiterer 
Konstanten unibersichtlich. 


5. Verwendete Bezeichnungen. Abb. 3 zeigt einen GrundriB des Fahrzeugs und der Spuren der | 
Rader auf der Fahrbahn, wobei die Kriimmung der Spuren und die GréBe der Winkel iibertrieben 


eingezeichnet sind. 

Es bedeuten: 

C(m1s?] =m/4k = I,/l?k den ,,Tragheitsbeiwert“; 

g[ms~?] die Erdbeschleunigung; 

G [kp] das Fahrzeuggewicht; | 
i [1] — /— 1 die imaginare Einheit; 
I, [kgm*] die Drehmasse eines Vorderrades um seinen Lenkzapfen; 

I, [kgm?] die Drehmasse des Fahrzeugs um seine Hochachse; 

ky [1] die Seitensteifigkeit eines Rades bei 1 kp Radlast; 

k [kp] = k,G/A4 die Seitensteifigkeit eines Rades bei der statischen Radlast, 
k,[kpm] den Riickstellbeiwert der Lenkung; 

1 [m] den Radstand; 

m [kg] die Fahrzeugmasse ; 

n [m] den Nachlauf der Vorderrader; 

N [kp] die Radlast; 

s [m] die Spurweite; 

S,, S. [kp] die Seitenkraft am Vorderrad, bzw. am Hinterrad; 

t [s] die Zeit ; 

v[ms—'] die Fahrgeschwindigkeit ; 

x [m] die Schwerpunktskoordinate in der Fahrtrichtung vom ortsfestem System; 
y [m] die Schwerpunktskoordinate der Seitenschiebung von ortsfesten System; 


¥1>¥_[m] die Koordinaten der Bodenberithrungspunkte von Vorder- bzw. Hinterrad. 
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¥y% [1] die Schraglaufwinkel des Vorderrades, bzw. des Hinterrades; 


y [1] die Spreizung; 

5 [1] den Nachlaufwinkel; 

p [1] den Lenkwinkel ; 

p (1) der Gierwinkel; 

y [s-1] die Kreisfrequenz ; 

F(s-*] das Argument der charakteristischen Gleichung ; 

-(v) [1] = AN/N die Anderung der Radlast durch den Luftwiderstand. 


Die Differentiation nach der x Koordinate wird durch Striche, die Differentiation nach der 
Zeit durch Punkte bezeichnet. Der Schnittpunkt des lotrechten Raddurchmessers mit der Fahr- 
bahn wird Bodenberiihrungspunkt genannt. 


6. Die Bewegungsgleichungen fiir das Fahrzeug. In Abb. 3 sind die Radspuren als Funktionen 
y(x) des Fahrweges eingezeichnet. Da man wegen der vorausgesetzten Kleinheit der Winkel 
sina + a, cosa % 1 und wegen der schwachen Kriimmung der Radspuren 1/9 ~ y’’ setzen kann, 
ergeben sich der Schwerpunktsatz in der Richtung der Kurvennormalen und der Momentensatz um 
die Hochachse zu 


—my" vw = 2(S, + S.) =2k (a, + a), — I,y =2(S,—S,) 1/2 = kl (o,—a,). (1) 
Die Neigung der Radspuren in den Bodenberiihrungspunkten der Rader ist y,, bzw. y;. Damit 
wird 


, ¥, see , — 
0 a ee pe Be Oy = ¥g— > 2. (2) 


Wird im Argument der Weg durch die Zeit ersetzt und y,, bzw. y, durch die Seitenschiebung des 
Schwerpunkts y und die Gierwinkel y ausgedriickt, so folgt annahernd 


ney rs typi, Joey +s/2—ylf2, 
Neyt+yle, Yet y— pl 
sowie RO Net we Vg VG OY SY Ye, 


Damit wird aus (1) 
e PANS = aja ; ; 
sm j= 2k LF —2y +), ae Babs (ID +9). (1a) 


Unter der Annahme, daB die Fahrzeugmasse in zwei Punkten zwischen den Radpaaren konzen- 
triert werden kann!, wird I, = m/I?/4. 
Mit 
Fhe UO 1 m Bat! 
4k PR 4,64 hey 


vereinfachen sich die Bewegungsgleichungen ae 


5S ae eh ead ee 
Gy ep y ytt=0, Clo+ [pt+e=0, (Lb) 


worin als Konstante nur die Fahrzeuggeschwindigkeit, der Radstand und der ,,Tragheitsbeiwert“ 
C vorkommen, der durch die Reifeneigenschaften bestimmt ist. 

Da die Gleichungen (1b) den Lenkwinkel y enthalten, tritt zu ihnen noch die Bewegungsglei- 
chung der Vorderrader um ihre Lenkzapfen, die zunachst méglichst einfach an gesetzt wird. 

7. Elastische Bindung der Vorderrader an ihre Geradeaus-Stellung. Unter der Annahme statio- 
narer Vorwartsbewegung des Lenkzapfens und bei Vernachlassigung der Normalbeschleunigung 
kommt mit der Drehmasse I, eines Vorderrades um seinen Lenkzapfen zu (Ib) noch 


21,6 +k py =0. (Le) 


Da darin die anderen Koordinaten y und y nicht vorkommen, beeinfluBt zwar die freie unge- 
dampfte Schwingung der Vorderrader die Bewegung des Fahrzeugs, letztere hat aber in erster 
Naherung keinen Einflu8 auf das Lenksystem. Erst die spater vorzunehmende Kinfiihrung des 
Nachlaufes wird diese einseitige Kopplung in eine beiderseitige verwandeln. Tatsachlich tritt 
allerdings schon durch die Reibung im Lenkzapfen eine weitere Kopplung hinzu. 

— 2 
1 Wenn man diese Annahme weg!a8t und allgemein I, = wu we » 4 + 1 setzt, andert sich die Form der 


Endergebnisse nur wenig (s. dort). a 
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Wenn die drei Gleichungen (1b) und (Ic) miteinander vertraglich sein sollen, mu8 ihre Koeffi- 
zientendeterminante verschwinden. Mit den Ansatzen 
aoe Vy Oe ee 
kommt 
| sh i ae ly 
OP = eo 9 | | 
1 = 
0 [ctor +o 1 0 (3) 
0 0 [2 I, 0% + k,] 
Daraus folgt die charakteristische Gleichung sechsten Grades in o 
a, o® + a, o® + a, o* + a, 0% + a,0? + a,0 + 4, = 9 (4) 
mit den Koeffizienten 
2Cl l 
a = C1.21,, “1=— a1 a,= Clk, +52 1,, | 
2Cl ly. mad (5)) 
a;/= a a= 7k, d= 4.0% 
Dieses System erfiillt nicht die Hurwitzsche Stabilitatsbedingung durchweg positiver a;.. 


Streng genommen ist die gesamte Bewegung des F'ahrzeugs weder stabil, noch stabilisierbar, da in| 
(Ib) bzw. (3) bei den Koordinaten y und yp die Riickstellglieder fehlen. Das bedeutet, da weder der! 
Schwerpunkt des Fahrzeugs an die x-Achse in Abb. 3, noch die Richtung seiner Laingsachse an die! 
x-Richtung als Mittellagen gebunden sind. Ware das anders, so ware das Fahrzeug nicht lenkbar, 
da ja dann auch ein Einschlagen der Vorderrader durch den Lenker diese Bindungen nicht auf- 
heben kénnte. 

Die Riickstellzahl der Lenkung k, ist sehr klein bei groBer Elastizitat des Lenkgestanges oder 
direkter Lenkung, und kann bei starrem Lenkgestange und hoher Lenkiibersetzung ziemlich grof 
werden!. Man soll daher das Riickstellen der Lenkung méglichst nicht der veranderlichen Elastizi-| 
tat des Lenkgestanges tiberlassen. 


8. Zur Definition der ,,Lenkstabilitat. Das Auffinden simtlicher Lésungen der charakteri- 
stischen Gleichung (4) ist hier auf ganz elementarem Wege méglich. Aus (1c) folgt, daB die Eigen- 
frequenz der Lenkung vy = yk,/2 I, auch eine Eigenfrequenz des Systems (3) sein muB, so daB 


oq IF 
Oia = ae il fe (6) 
das erste Lésungspaar ist, da in algebraischen Gleichungen komplexe Lésungen stets konjugiert 
auftreten. 
Da man auBerdem Gleichung (4) durch o? kiirzen kann, kommt dazu die Doppelwurzel o;,, = 0. 


Wenn man schlieBlich aus den bereits gefundenen Lésungen die Linearfaktoren bildet und durch 
deren Produkt (o — o,) (o —o,)o* die Gleichung (4) dividiert, so kommt 


1\2 
(Cot, )2t,1=0 
und damit die Doppelwurzel 


il 
O34 = CH . 


Die allgemeine Lésung des Systems (3), z. B. fiir die Koordinate y angeschrieben, lautet dann 
i wes = ae 1 = 
7 Ty 2 Ty cae Cate Sr Gok ; 
y = be + bye + be + b,te +6,+ bet. (7) 
Die Glieder mit dem reellen negativen Exponenten klingen rasch ab und brauchen ebenso wie 


das absolute Glied nicht beriicksichtigt zu werden. Die ersten beiden Glieder ergeben eine unge- 


dampfte Schwingung mit der Kigenfrequenz der Lenkung und das letzte Glied eine Seitenschiebung 
mit konstanter Geschwindigkeit. 


Fir die Koordinate p gilt das entsprechende, hier stellt das letzte Glied eine Drehung des 
Fahrzeugs um seine Hochachse mit konstanter Drehgeschwindigkeit dar. Man erkennt in den 


1 Giche Bubuote 1 von Seite 88. 
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sechsten Gliedern den EinfluB des zu Beginn der Arbeit erwahnten exzentrischen StoBes auf den in 
der Fahrtrichtung starr mit dem Fahrzeug verbundenen Lenkzapfen. So klein dieser direkte EinfluB 
auch ist, so kann er nur durch das Eingreifen des Lenkers ausgeglichen werden, dem dazu geniigend 
Zeit zur Verfiigung steht. 

Ein zweiter Drehsto8 um die Hochachse erfolgt, wenn das Hinterrad dieselbe Unebenheit iiber- 
fahrt. Auch auf guter Fahrbahn iiberfahren alle vier Rader stindig sehr kleine Unebenheiten, so 
daB sich die erhaltenen Drehimpulse nahezu aufheben, bzw. nur langsam summieren. Auf jeden 
Fall verhindert die groBe Drehmasse des Fahrzeugs um seine Hochachse eine schnelle Anderung 
der Fahrtrichtung. 

Viel gefahrlicher ist der indirekte Einflu8 des DrehstofSes um den Lenkzapfen eines Vorderrades, 
da bei der geringen Drehmasse des Lenksystems sein Einschlag verhaltnismafBig groB wird und 
dadurch die groBen Tragheitskrafte der gesamten Fahrzeugmasse ins Spiel bringt. 

Ein Fahrzeug soll demnach schon dann als stabil angesehen werden, wenn der durch eine kleine 
Stérung verursachte Einschlag des Lenksystems aperiodisch oder unter gedampfter Schwingungen 
guriickgeht. AuSferdem diirfen die durch den Einschlag der Lenkung ausgelésten Bewegungen 
des ganzen Fahrzeugs, die gleichzeitig und in ahnlicher Weise ablaufen, keine zu groBen Amplituden 
erreichen. Es ist leicht genug, die erste Bedingung allein zu erfiillen, doch ist die Erfiillung der 
zweiten Bedingung praktisch viel wichtiger, wie der folgende Abschnitt zeigt. 


9. Berechnung der stationiren Amplituden. Nachdem in Gleichung (7) fiir die Stabilitatsbe- 
trachtung nur die rein imaginaren Liésungen iibrigbleiben, fiihrt das Fahrzeug nach einem sehr 
kurzen Einschwingsvorgang harmonische Schiebe- und Gierschwingungen mit der Kreisfrequenz 


des Lenksystems » = |/ k,|2 I, aus. Ihre Amplituden yo, y) und Phasenlagen ®, y gegeniiber dem 
Lenksystem kénnen mit den Ansatzen 

P=Ysin(yt), y=yosn(vt—D), yw=wposin(vt— x) 
durch Koeffizientenvergleich errechnet werden. Nach EHinsetzen in die Bewegungsgleichungen 
(1b) kommt: 
| — Cr® yo sin (vt — D) += 9yq cos (vt — D) — y sin (vt—y) + sin vt =O, 
1 (8) 
— C1 yp sin (vt — x) +—- 7 Yo cos (vt — x) + Yosinvt =O. 
Dain beiden Gleichungen die Summe des Glieder mit sin vt, bzw. der Glieder mit cos vt fiir sich 
verschwinden mu, folgt fiir die Koeffizienten: 


— Cr yy cos D + —v yp sin B — yy cos y + =0, (1) 
) Cory, sin D + —vyy cos D + yysing =0, (IT) 
) — C1 yp cos 1 + 9 yo sin x +H = 0, (IIT) 
Cl v% yy sin x + ~ vy cos x = 0. (IV) 

Aus (IV) ergibt sich ; 
Sy Pe arena: (9) 


{ 


Wenn man (III) und (IV) quadriert und addiert erhalt man die Amplitude des Gierens 
——— 10) 
Yo ~ Ty W(Crve +1 ( 
Sie ist fiir groBe v von der Geschwindigkeit praktisch unabhangig (Abb. 4); groBer Radstand und 
hohes » sind giinstig, da der Radstand hier die Frequenz nicht beeinfluBt. 
_ Durch Quadrieren und Addieren von (I) und (II) ergibt sich mit cos y aus (9) die Amplitude der 
Seitenschiebung des Schwerpunkts zu 


ae 4 Pw + Boru, (11) 


-worin 
u=([(Cyv)?+ 1), w=1—IPC 
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ist. Thr in Abb. 5 eingezeichneter Verlauf kann, wie folgt, gedeutet werden: Ohne Schraglauf der 
Rader mii®ten bei einer bestimmten Frequenz der Lenkung die Amplituden der Vorderrader und 
damit auch die des Schwerpunkts yp etwa proportional mit der Fahrgeschwindigkeit wachsen. Kin 
ahnlicher Verlauf zeigt sich aber nur bei langsamer Schwingung des Lenksystems und kleiner Fahr- 
geschwindigkeit, wahrend bei gréBerem v und »v infolge der steigenden Massenkrafte groBe Schrag- | 
laufwinkel auftreten, die die Wellen der Radspuren und der Schwerpunktsbahn stark abflachen. 

Es zeigt sich hier deutlich, da eine geniigend 
hohe Frequenz der Lenkschwingung viel wichtiger 4 Yo/% [m/rad] 
ist, als ihre starke Dampfung. Trotz des unge- @ 
dampften Lenksystems ist das Verhalten des F'ahr- 
zeugs als durchaus stabil zu bezeichnen,'solange 
y > 6,28 s ist. Der Grund,{warum die: ela- 


Yo|%o 


v=4715" 


V=62057 


u[m/sek] 


0 10 20 50 0 10 20 80 
Abb, 4, Gieramplituden bei ungediimpfter Lenkschwingung. Abb. 5. Amplituden der Seitenschiebung bei ungedampfter 
Lenkschwingung. 


stische Bindung der Vorderrader an ihre Geradeaus-Stellung nicht zur Stabilisierung verwendet 
wird liegt wohl darin, daB eine zur Erzielung eines hohen y geniigend harte Feder das starke Kin- 
schlagen der Rader bei langsamer Fahrt zu sehr erschwert. 

Fiir Sonderzwecke, wie z. B. die Sicherung eines extrem schnellen Fahrzeugs gegen seitliche 
WindstéBe, diirfte die elastische Bindung in Betracht kommen, da in diesem Fall der Nachlauf 
zumindest zu Beginn der Stérung, eine recht bedenkliche Wirkung hat, wie noch spater naher | 
ausgefiihrt werden soll. 

Der Vollstandigkeit halber sei hier noch der Phasenwinkel zwischen Lenk- und Schiebeschwin- 
gung angefiihrt: 


__ Gy £/x(4 + C2 y2 aL k) 
vD Vv 
@y—K ~~ 5 


te Oe (12) 


mit 


4 
~ 92 (402 w + 2 »? u) 


K 


Das doppelte Vorzeichen in (12) fiihrt nicht auf zwei verschiedene Amplituden. Die Phasenver- 
schiebung zwischen Schieben und Gieren ist in beiden Fallen die gleiche, es andert sich nur die 
Reihenfolge. 


10. Der Nachlauf der Vorderrader. Werden die Lenkzapfen so geneigt, daB8 ihre Mittellinien 
die Fahrbahn vor den Beriihrungspunkten der Vorderrader treffen, so erhalten die Seitenkrafte S 
auf die eingeschlagenen Rader Momente um die Lenkzapfen, die dem Einschlag entgegenwirken. 
Abb. 6 zeigt eine vereinfachte Darstellung dieses Sachverhalts. 
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Fiir ein Vorderrad wird das Riickstellmoment um den Lenkzapfen wegen (2) 


nS) = nko =n bly, AT + o\=nk 


hve [Bae 
; (V+ 59)-v +9]. (13) 
Mit I, p + nka, = 0 und (13) kommt das System 


x Ubi ye 
CH +—5 eee | 


y (14) 
Ot Lit yb eee 
dessen Koeffizientendeterminante 
Co? + 2 0 —1 = 
0 JCro* +o} sel 00, 9a (15) 
a ree — ot 


sein mu. Daraus folgt die charakteristische Gleichung sechsten 
Grades ino 


dy o® — a, o® ae a, o4 ae as o3 ae a, o2 a azo a a, = 0 ; (16) Abb. 6. Schema der ect des Nach- 


laufes 
mit 
IE XG Ak ee 
~ 2 P ead Se 2 pa ah 
% Bias 4 v nk’ EP eI Bear git (17) 
I 
ele Gia G6 a. 5 a, = 0". 


Da man (16) durch o? kiirzen kann, kommt zunichst wieder die Doppelwurzel o,,, = 0, die in 
der allgemeinen Lésung des Systems (14) auf ein absolutes und ein mit der Zeit linear anwachsendes 
_Glied fiihrt. Das erstere ist uninteressant und das letztere wurde weiter oben als die Folge des 

exzentrischen StoBes auf das ganze Fahrzeug gedeutet, die trotz ihrer Kleinheit nur durch den 
Lenker behoben werden kann. Bei der Beurteilung der Lenkstabilitat miissen daher diese beiden 


Lésungen ausgeschlossen werden. 
) Auf die durch o? dividierte Gleichung. (16) 


a,o* + a,0° + a,0* + a,0 + a,=—0 (18) 


| kann man das Stabilitatskalkiil nach Hurwitz anwenden, das auBer durchweg positivem a; in diesem 
Fall verlangt, daB 


Ar= ae 0, A, = 050 ApS) a, a, “a, |\0 
| OQ a,. 

und A, =a,A, > 0 sein miissen, wenn jede Wurzel von (18) einen negativen Realteil haben soll. 
Dann setzt sich die allgemeine Lésung des Systems (14), die ahnlich wie (7) angesetzt werden kann, 
aus vier mit der Zeit abklingenden Gliedern zusammen. 


Unter der Annahme, daB sich das System an der Stabilitatsgrenze befindet, d.h., daB rein 
imagindre Wurzeln auftreten, muB die Summe der Glieder mit den ungeraden Potenzen von o in 
(18) fiir sich verschwinden. Es wird a, 0% + a,0 = 0 und daraus mit (17) 


x ink 
O12= 41 4 i/et. (19) 


Mit einem dieser Werte fiir o muB jetzt auch die Summe der Glieder mit den geraden Potenzen von 
-o und des absoluten Gliedes verschwinden. Es muf 


Ig [nk\2 Palo nk e. 
ork (he (ors tk) Gt +e~0 
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n= PRL TB te ) ) > ( ) 


eine Bedingung, die mit A, = 0 identisch ist. Wenn auBerdem A, > 0 und A, > 0 erfiillt sind, | 
kann man erwarten, daB die beiden restlichen Wurzeln negative Realteile haben, d. h. daB durch 


(20) die Stabilitatsgrenze in der n —v Ebene festgelegt ist. Die Bedingung A, = a,> 0 fiihrt mit 
(20) in (17) auf C—1/2v? > 0, der Nachlauf muB also stets positiv sein. Die Bedingung Age a, dy — 
ay 43 > 0 fiihrt auf C> 0. Da in (17) a, = C ist, sind beide Bedingungen schon in der Bedingung 
durchweg positiver a; enthalten. Wegen o,,, = iv ergibt sich mit (20) in (19) die ungedampfte 
Figenfrequenz v an der Stabilitatsgrenze als Funktion der Fahrgeschwindigkeit v zu 


a 

| (oe eee 21 

Var Ce en 

Die Stabilitatsgrenze (20) in der (n, v)-Ebene erreicht bei v = 1/2 C die Abszissenachse, dort wird 
auch vy = 0. Fir kleinere v ist die Abszissenachse selbst Stabilitatsgrenze. 


Die beiden restlichen Wurzeln von (18) ergeben sich, wenn man (20) in (18) einsetzt und 
durch das Produkt der zu o,, gehdrigen Linearfaktoren 


nk 
(¢—=0;)(¢—¢,) =o Wy 
dividiert, zu 
1 ‘ 2 1 
| eee 
wie erwartet, mit negativem Realteil. 
Aus (20) wird fiir v = 0 


at 
Tr Cine 


(20a) 


Mit diesem oder einem gréBeren Nachlauf ist die Lenkung bei allen Geschwindigkeiten stabil. Fir | 


I, = 0,06 kg m?, C= I/kopg = 0,02 ms?, J= 2m und k = 1000 kp wird h = 0,006 m = 6 mm. 
Schon der recht kleine Nachlauf von 6 mm (das sind fiir ein Rad von 600@ etwa 6 = 1°10’ nach der 
iblichen Bezeichnung) verhindert angefachte Schwingungen. Da einerseits der ausgefiihrte Nach- 
lauf nach de Lavaud 2° bis 4° betragt und andererseits nach Becker, Fromm und Maruhn die 
Wirkungslinie der Seitenkraft nicht durch den geometrischen Bodenberithrungspunkt des Rades, 
sondern etwas dahinter vorbeigeht?, ist der mechanisch wirksame Nachlauf der Kraftwagen jeden- 
falls gréBer als nach (20a) errechnet wurde. 

Die in Abb. 7 als Funktion von n und v eingezeichnete Stabilitatsgrenze trennt das Gebiet der 
gedampften von dem der angefachten Schwingungen, ohne etwas iiber die zugehérigen Amplituden 
auszusagen. Im vorigen Abschnitt wurde gezeigt, daB die GréBe der ungedampften Schiebeampli- 
tuden mit sinkender Frequenz stark zunimmt, was auch fiir die gedampften gelten muB. 

Ks ist daher nachzupriifen, ob im stabilen Bereich in Abb. 7 schwach gedampfte Schwingungen 
mit niederer F'requenz auftreten, die unmittelbar nach der Stérung auf groBe Amplituden der Sei- 
tenschiebung fiihren. Starke Dampfung bei der gleichen Frequenz ist nicht so bedenklich, da 
durch sie der Kinschlagwinkel y der Vorderrader beim Uberfahren der Stérung verringert wird, 
dem die Schiebeamplituden proportional sind. Deshalb ist auch ein aperiodisches Zuriickgehen des 
Kinschlagwinkels nicht gefahrlich, wenn die sich daraus ergebende Kursanderung zusammen mit 
der durch den exzentrischen Sto um die Hochachse bewirkten vom Lenker behoben wird. Nieder- 
frequenten Schwingungen des Lenksystems und damit des ganzen Fahrzeugs kann der Lenker je- 
doch nicht leicht entgegenwirken. Das neben dem Verlangsamen der Fahrt beste Mittel, eine flache 
Kurve zu fahren, ist meist nicht anwendbar. 


11. Die konjugiert komplexen Lésungen der charakteristischen Gleichung. Da nur die periodi- 
schen Lésungen von (18) interessieren und diese konjugiert auftreten, entstehen nur zwei Schwin- 
gungen mit den Frequenzen y,,, und », ,, denen die paarweise konjugiert komplexen Lésungen 0, » 
und o;,, der Gleichung vierten Grades entsprechen. Der Weg, diese Gleichung fiir beliebig vorge- 
gebene Werte von n und v zu lésen und so ein Bild von der Verteilung der Frequenzen , , und ¥5 , 


1 Bir w+ 1 ist die rechte Seite von (20) mit dem Faktor 1/u zu multiplizieren; (19) bleibt unverandert. 
Die Bedeutung von y ist aus FuBnote 1 auf Seite 91 zu ersehen. 
2 Der dadurch bewirkte sog. ,,Reifennachlauf wachst mit zunehmender Radlast, bzw. mit abnehmendem 
Innendruck des Luftreifens an und kann bei Niederdruckreifen die GréGe von einigen Zentimetern erreichen. 
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und ihrer Dampfung im stabilen Bereich zu erhalten, ist jedenfalls sehr langwierig. Es soll daher 


versucht werden, der Gleichung (18) ahnlich wie beim Aufsuchen der Stabilitatsgrenze eine bestimmte 
Art von Lésungen vorzuschreiben 


und als Nebenbedingung dafiir 2, : 3 
einen Zusammenhang zwischen /\ { 
nund v zu erhalten. Dadurch kén- 74-90 
nen unter Umstinden innerhalb iA 
‘des stabilen Bereichs in Abb. 7 
eine oder zwei Linien festgelegt Eon 
werden. Auf diesen Linien kénnen — oe V49=155" 
wieder Punkte zu vorgegebenen ”% 
Frequenzen aufgesucht und die 
zur gleichen Frequenz gehérigen 
-Punkte miteinander verbunden 
werden. Allerdings ist ein solches 
Verfahren nur bei Gleichungen 
vierten Grades mit einer gewissen 
Symmetrie der Koeffizienten an- = 
wendbar, wie sie (18) tatsachlich V42=10 ae 
aufweist. 

Zunachst sollen zweizweifache, 
konjugiert komplexe Wurzeln als 
Lésung vorgeschrieben werden. 3 oe 

Dann isto,,,=03,4=%, t1v.Das_ 5 = SS SO Habililttsgrenze 

Produkt der ersten zwei Linear- ‘ 

faktoren wird [o — (x, +1 7)] y 

[o— (x, —iv)] = o?—20x%,+ 

x? +? und wenn man diesen Aus- S aay 

-druck quadriert und nach Poten- 
zen von o ordnet, kommt 


v[m/sek] 
ee 
of —4x, 0° + (627 +27) 0? 0 70 20 30 40 
3 2 Abb. 7. Stabilititsgrenze und Linien konstanter Eigenfrequenz abhangig 
=k (x7 35 a a ) o on Fahrgeschwindigkeit und Nachlauf. 
4 (2 +r =0. 
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Man kann nun durch Koeffizientenvergleich mit der durch a, dividierten Gleichung (18) x, und 
berechnen. Da aber mehr Gleichungen zu erfiillen als Unbekannte vorhanden sind, erhalt man als 


Nebenbedingung eine Beziehung zwischen den Koeffizienten. d.h. hier zwischen den darin auftre- 
-tenden Parametern n und v. 


_ Mit (17) wird 
ee ee ee 
a a Ay 1? 1 Ga 
a, nk 1 . ea sae oO nk I 
Ge ls Ce Ese aE asec ae Tp 2C?e?}? 
SSN ae 
OW 2 
a3 1 wh 2 2 
eR Ce genie) 8 
Gg i mk _ (424 422 
ine CLL, eg 
Da 
1 nk 1 nk 1 — th aalhrGa 
2 oe age oe eeeeaile = 
cs Nae eA erh Te Cont, (3 Cv Ip 2 
sein muB, so folgt = 
1 nk n? k? ? : 
ee n =——- = const. (20a) 
ah ee ee Clk 


Das ist in Abb. 7 eine waagerechte Gerade und Asymptote an die Stabilitatsgrenze (20). Auf 
ihr hat Gleichung (18) die Wurzeln 


1 .j/nk ae 1 Va 91 
O,2=o6= — gen ti | Peas Korie CE AG (243) 
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Die Kreisfrequenz v wird Null fiir v = /1/8 C, fiir kleinere v existieren solche speziellen Lésungen 
nicht. 

Danach soll untersucht werden, ob und fiir welche Werte von n und v Liésungen von (18) existie- 
ren, die aus zwei konjugiert komplexen und einer reellen Doppelwurzel bestehen. Das Auftreten | 
der letzteren zeigt den Ubergang von zwei konjugiert komplexen zu zwei verschiedenen reellen 
Wurzeln an. Die zugehdrige Kurve in Bild 7 bedeutet mechanisch die ,,Schwingungsgrenze* 
zwischen gedampfter Schwingung und aperiodischem Ubergang. 

Es soll 6,5 =%, Liv, 634 =x, sein. Das Produkt der Linearfaktoren ist [o — (x, +177)] 
[o — (x, —i»)] (o — x, )®, und nach Entwickeln und Ordnen nach Potenzen von o wird 


ot — 2 (x, + x) 08 + (x2 + v2 + 08 + 4 x, x9) 0? — [2 x, 228 +2 x, (xt + *)] o + (xf + 7”) 2g = 0 
Der Koeffizientenvergleich mit (18) ergibt 


=e t+r4 28+ 4%, %,, 


az 1 nk 2 
a Cee ag heel ae 
a 2 2) 42 
Oni ee 
Mit 
1 
eee ae 


wird aus der zweiten, bzw. dritten Gleichung 


{I 2 Xo 
pasate ee pee ee Ny 2 — pee 
City 1a ey Pe eee 


2 2 x9 1 nk 
Sie es 


Wenn man die letzte Gleichung mit Cv erweitert und von der vorletzten abzieht, folgt 


w+t2xn,vCv=0, x,=—1/2Cv und nach Einsetzen dieses Wertes fiir x, in die vorletzte 
Gleichung 
Tike eto 
; fe / = 2 C22 
sowie endlich 
ee eae 
LT OG a 


Die Werte fiir x,, 7 und x, in die vierte Gleichung eingesetzt ergeben als Nebenbedingung eine Be- 
ziehung zwischen den Parametern n und v 
if Cl 


ne 4 C80? (C1— 4 C20’) ° 


(20b) 


Das ist die Gleichung der ,,Schwingungsgrenze“ in Abb. 7. Auf dieser Kurve hat (18) die Lésungen 


: eee 8 Civi— Cl 
oO — = — —— 
tome toi | eS 
l 


O34 =X = 


(21) 
Cr 


Das Minimum dieser Kurve liegt bei v = y1/8 C, wo sie von der Geraden nach (20a) berihrt wird. 
In diesem Punkt (D in Abb. 7) hat die Gleichung (18) eine vierfache negativ reelle Wurzel. Die Real- 
teile der Lésungen sind auf der Kurve und Gerade von vornherein gleich und in Punkt D verschwin- 
den wegen vy = 0 die Imaginarteile, die in Anzahl und GréBe auf beiden Linien verschieden sind. 
Auf dem rechten Ast der Kurve ist y reell, dort gibt es auSer der negativ reellen Doppelwurzel zwei 
konjugiert komplexe Wurzeln, d. h. dieser Ast ist die Schwingungsgrenze fiir eine der beiden Kigen- 
frequenzen und zwar fiir y; ,, wahrend Y,,. dort noch auftritt. Auf dem linken Ast wird y imaginar, 


und die zwei konjugiert komplexen Wurzeln verwandeln sich in zwei negativ reelle. Dieser Ast ist 
die Schwingungsgrenze fiir 7, .. 
+ 
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Zu einer bestimmten Higenfrequenz y findet man den auf der Stabilitatsgrenze (20) liegenden 
Punkt mit 
l 


eee ee nt 
res 4h Cie 2) 
den auf der Geraden (20a) liegenden Punkt mit 
l 
ee aie 
‘Sa CQ—sCh Pee 
und den auf der Schwingungsgrenze liegenden Punkt mit 
SS GO SE /4+ C2 yA 
ae <I? iin Se. Ge) 


Auf den beiden ersten Linien sind die beiden Kigenfrequenzen gleich y = y, 5 = v3, 4, auf dem rech- 
ten Ast der Schwingungsgrenze ist y = 7, ,, v3 4 = 0. Zwischen Punkt D und ‘dem Schnittpunkt von 
Stabilitatsgrenze und Abszissenachse konnte die Schwingungsgrenze mit dem oben geschilderten 
Verfahren nicht als Funktion von n und v dargestellt werden. Sie fallt dort fiir beide Frequenzen 
71,, und yz , zusammen, so da durch den K oeffizientenvergleich nur zwei Unbekannte zu bestimmen, 
dafiir aber zwei Nebenbedingungen zu erfiillen sind. 

In Abb. 7 wurden die Linien gleicher Eigenkreisfrequenz fiir y,,. = 0, 5, 6, 7, 10 und 15 s— einge- 
zeichnet, nachdem zur Kontrolle fiir einige beliebig gewahlte Punkte in der n-v Ebene die Wurzeln 
der Gleichung (18) mit einem von A. Leonhard! angegebenen Naherungsverfahren berechnet 
wurden. Die Linien », , = const verlaufen bei héherer Geschwindigkeit praktisch parallel zur Ab- 
szissenachse. Die Linien v, , = const breiten sich strahlenférmig im stabilen Bereich zwischen Sta- 
bilitatsgrenze und Schwingungsgrenze aus. Fiir ein gegebenes v nimmt im stabilen Bereich mit 
wachsenden n die eine Kigenfrequenz y,, und ihre Dampfung zu, die andere Eigenfrequenz 1,,, 
nimmt dagegen ab. Fiir ein gegebenes n nimmt mit wachsenden v die Dampfung allgemein ab und 
beide Eigenfrequenzen nehmen geringfiigig zu. 

Daraus folgt, daB der Nachlauf der Vorderrader etwas grifer als der Wert nach (20a) n= 4 J,/ 

Clk ausgefiihrt werden soll. Bei kleinerem n geht die Lenkstabilitat bei héheren Geschwindig- 
keiten verloren, wahrend sie durch erhebliche VergréBerung von n zumindest nicht wesentlich ver- 


bessert werden kann. 
=) Z 1 
Sy Cle Cie? 


Nach (21) ist 
d. h. ein groBer Radstand ist ungiinstig, da er die Eigenfrequenzen des Fahrzeugs herabgesetzt und 
die Amplituden des Schiebens und Gierens erhéht, wie man erkennt, wenn man mit y aus (21) in (11) 

bzw. (10) eingeht. 
Mit (20a) und C = 1/k, g, k = ky G/4 kann als Richtwert fiir den Nachlauf zunachst 
cui 
aaa (20c) 

angegeben werden. In diesem Ergebnis ist jedoch der Luftwiderstand noch nicht beriicksichtigt, 
was bei hohen Geschwindigkeiten unerlaBlich ist. 


12. Der Einflu8 der Luftkrifte. Es wurde vorausgesetzt, da die resultierende Luftkraft bei 
Windstille der Schwerpunktsgeschwindigkeit des Fahrzeugs entgegengerichtet und deren Quadrat 
proportional sei von der Form W;,, = kz, v®. Um zu einem mittleren Wert fiir ky; zu kommen, wurde 
angenommen, daB bei v % 50 km/h der Luftwiderstand die GréBe des Rollwiderstands erreiche. 
Damit wird W,, =k, v? = G fp und mit G = 800 kp, fp = 0,02 wird k, = 0,08, was eher etwas 
zu hoch gegriffen ist. ; 

Mit dem Schwimmwinkel f zwischen Schwerpunktgeschwindigkeit und Fahrzeuglangsachse 
baw. genauer mit der Summe des Schwimmwinkels f und des Gierwinkels y multipliziert, hat der 
Luftwiderstand eine Komponente quer zur Fahrtrichtung. Da in der Einleitung die Schraglauf- 
winkel und damit auch f und yp klein vorausgesetzt wurden, bleibt die Komponente von Wy, sehr 
klein gegeniiber den Seitenkraften der Rader. Dasselbe gilt fiir den Kraftarm eines etwa auftreten- 
den Moments des Luftwiderstandes um die Hochachse des Fahrzeugs. Bei Windstille kénnen diese 
Einfliisse kaum zur Anfachung kleiner Schwingungen beitragen, was jedoch nicht gegen die grofe 
Gefahrlichkeit von seitlichen WindstiSen oder des Luftwiderstandes auf ein bereits aus der Bahn 
geratenes Fahrzeug mit groBem Schwimmwinkel spricht. 


1 A. Leonhard, Die selbsttatige Regelung, Berlin 1957. 


100 G. Mitterlehner: Die Lenkstabilitat des luftbereiften Kraftwagens Ingenieur-Archiy 


Der Luftwiderstand hat aber noch eine ganz andere Wirkung, die iiberhaupt nicht von den 
Amplituden bereits vorhandener Schwingungen abhingt. Da seine Wirkungslinie jedenfalls in 
einem erheblichen Abstand iiber die Fahrbahn verlauft, entlastet er die Vorderrader und belastet 
die Hinterrider in ahnlicher Weise wie die Tragheit beim Beschleunigen des Fahrzeugs. Dadurch 
setzt er die Lenkstabilitat wesentlich herab, wie der Grenzfall der vollkommen entlasteten Vorder- 
rider zeigt, die keine Seitenkrafte mehr aufnehmen kénnen. 

Die bezogene Anderung des Raddrucks eines Vorderrades durch den Luftwiderstand ist 

N—AN 0,25 G — 0,5 k, v? h/t 2h 

eo = =e ee i = 1— er Bal ubssiSe @) 
worin mit dem Abstand des Luftwiderstandes von der Fahrbahn hz = 0,75 m und den schon frither 
verwendeten Werten ¢ = 0,75- 10-4 v2 wird. Damit lauten die Bewegungsgleichungen des Fahr- 
zeugs unter Beriicksichtigung von (2) und (2a) und wegen k = ky N 


—my=2k[(1— Qa + (+ ag) =24( 5-29 +9 coo), 


(24) 


~ 17 = 25 [0 —Da— 1+ Dal =H Hp —2E SF + ly — E09). 


und mit C = m/4k = I,/l2k kommt nach Erganzung durch die Bewegungsgleichung der Lenkung 
das System 


3 a eas . 
CP yy (ee ed 


l 
v 


Clp+ip+26y—20-F +1 a (25) 
th l 


3 Lee ° 


dessen K oeffizientendeterminante 


Cotto een =e) 

a 8G Clot tlo+2¢] i= =) (26) 
1 l I 
red DE nE(Ga (wee 


sein soll. Die daraus folgende charakteristische Gleichung (16) vom sechsten Grade in o hat die_ 
Koeffizienten 


C4 aati ey OE il 
BS get ghey mga rayaet 4 awn ey 
I 
en Saale +26C— 2), Gs ie “ (l-p<¢)., ee 


a=C(1+@), a=—a,=0. 

Man kann die Gleichung daher wieder durch o? kiirzen, wobei die Bedeutung der Doppelwurzel 
05,4 = 9 schon bei der Besprechung der Gleichung (5) bis (7) erwahnt wurde. 

Ubrig bleibt die Gleichung vierten Grades (18) mit den Koeffizienten ay bis a, in (27), die wieder 
die Stabilitatsbedingungen [alle a; positiv, sowie A, > 0, A, > 0,4, >0 und A, > 0, siehe Glei- 
chungen (18) ff] erfiillen miissen, wenn (18) nur Wurzeln mit negativen Realteil haben soll. An der 
Stabilitatsgrenze sind auch rein imaginare Wurzeln zugelassen und dort mu wieder gelten: 


a, 0? + a,06 =0, ee eee) Ca (28) 
QP 


Mit der einen wie mit der anderen dieser Wurzeln mu dann auch der Rest der Gleichung (18) 
verschwinden, z. B. Mg oj + a, of + a, = 0, worin ay,-a, und a, aus (27) zu entnehmen sind. Daraus 
ergibt sich nach einigem Umformen die Stabilitatsgrenze als Beziehung zwischen den Parametern 
n, v und C(v): 

Pee th I 
= Coe (en) Ic a (4 0. 2% 


* Auch fiir (29) gilt sinngemaB FuBnote 1 auf Seite 96; (28) bleibt unverandert. 
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Durch Einfiihren des Nachlaufs n nach (29) in die Ausdriicke fiir die Koeffizienten (27) priift 
man, ob die anderen Stabilitatsbedingungen erfillt sind. 4, = a, > 0 fiihrt auf 


I 


d.h., daB nur der positive Ast der Kurve (29) Stabilitatsgrenze ist. Die Bedingung A, = a, a, — 
a, a; > 0 fiihrt auf C > 0. Die Gleichung der Stabilitatsgrenze (29) ist identisch mit J, = 0, weshalb 
auch A, = a, A, — 0 wird. 

Mit (29) in (28) erhalt. man den Zusammenhang zwischen der Frequenz der ungedampften 
Schwingung und der Fahrgeschwindigkeit an der Stabilitatsgrenze: 

2 l 
ene a eo). (30) 

Durch Umformen dieses Ausdrucks erhalt man fiir ein vorgegebenes y eine biquadratische Glei- 

chung fiir das zugehérige v an der Stabilitatsgrenze, da man ¢ = ¢ (v) beachten muB. 


[mm] 
Fits mm 


15 


ae 
Zale 0 

Ne “ade ent? 
—F granites 


Schwingungsgrenze Y, 4 =0 


unstabil 


v [m/sek] 


0 70 20 30 40 50 60 70 80 0 100 
Abb. 8. Stabilitatsgrenze und Linien konstanter Eigenfrequenz bei Beriicksichtigung des Luftwiderstandes. 


In Abb. 8 sind die Stabilitatsgrenze bei Beriicksichtigung des Luftwiderstandes nach Gleichung 
(29) und einige Linien konstanter Frequenz y,, iiber der Fahrgeschwindigkeit aufgetragen. Aller- 
dings wurden von den Frequenzlinien nur die Punkte auf der Stabilitatsgrenze (29) und auf dem 
rechten Ast der Schwingungsgrenze (20b) berechnet. Letzteres ist zulassig, da an der Schwingungs- 
grenze die Geschwindigkeiten so niedrig sind, daf dort der Luftwiderstand keine Rolle spielt. Des- 
wegen konnten auch die y, , = konst. Linien ohne Berechnung von weiteren Zwischenpunkten, die 
jeweils die Lésung einer Gleichung vierten Grades verlangt hatte, mit guter Naherung eingezeichnet 
werden. 

Bei hohen Geschwindigkeiten wird der stabile Bereich infolge der Entlastung der Vorderrader 
durch den Luftwiderstand stark eingeschrankt, wobei die Windschliipfigkeit des Fahrzeugs und das 
Verhaltnis der Héhe des Angriffspunktes der Luftkrafte zum Radstand das Ergebnis wesentlich 
beeinflussen. Auer fiir die bezogene Radlastanderung ¢ = 0,75 - 10~4 v? ist die Stabilitatsgrenze (29) 
auch noch fiir € = 0,5- 10-4 v? eingezeichnet, wobei die letztere Kurve einem geringeren Luftwider- 
standsbeiwert bzw. einer tieferliegenden Wirkungslinie des Luftwiderstandes entspricht. Fir 
v = 50 m/s, d. h. fiir eine Fahrgeschwindigkeit von 180 km/h, ergibt sich der Nachlauf an der Sta- 
bilitatsgrenze aus (29) mit € = 0,75 - 104 v? und den vorher benutzten Werten zu n = 7,2 mm; das 
entspricht fiir ein Rad von 600 @ nicht ganz 1,5 Grad Riickwartsneigung des Lenkzapfens. 


13. Die obere Grenze des Nachlaufs. Bei der technischen Ausfiihrung des Nachlaufs durch Nei- 
gung der Lenkzapfen der Vorderrader (Abb. 1) kann schon aus kinematischen Griinden n nicht 
sehr groB werden. Andernfalls wiirde beim Einschlagen der Vorderrader der Schwerpunkt des 
Fahrzeugs gesenkt und bei Geradeaus-Stellung der Lenkung ware zumindest das stehende Fahr- 
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zeug in einer labilen Gleichgewichtslage. Das kénnte allerdings bei einer Ausfiihrung nach Abb 6 
mit in der Fahrtrichtung nicht geneigten Lenkzapfen vermiedern werden. Es gibt aber noch einen 
zweiten wichtigen Grund, den Nachlauf nicht iiber das notwendige Maf hinaus zu vergréBern. 


Bisher wurde stets von einer Stérung durch eine Unebenheit der Fahrbahn ausgegangen, die 
zunichst die Vorderrader auslenkt und damit indirekt die Bewegung des Fahrzeugs beeinfluBt. 
Ein seitlicher WindstoB wirkt dagegen zunachst auf den Aufbau des Fahrzeugs und wird zum Teil 
iiber die Lenkzapfen auf die Vorderrader iibertragen. Daher entsteht infolge des Nachlaufs ein Mo- 
ment, das die Vorderrader nach der falschen, dem Wind abgewandten Seite einzuschlagen sucht, 
wihrend ein Einschlagen gegen die Windrichtung zur Erzielung der ausgleichenden Seitenkraft er- 
forderlich ware. Je gréBer der Nachlauf, desto empfindlicher reagiert das Fahrzeug auf seitliche 
WindstéBe. Nach Meinung des Verfassers laBt sich dadurch die bekannte Empfindlichkeit gegen 
Seitenwind gerade vieler schneller, windschnittiger Fahrzeuge zumindest teilweise erklaren. 


14. Der Auftrieb. Im allgemeinen wirkt die resultierende Luftkraft auf einen Kraftwagen nicht 
genau parallel zur Fahrbahn, so daB zum Luftwiderstand noch ein kleinerer Auftrieb hinzukommt. 
Dieser sei als Kraft definiert, die die Radlasten aller Rader des Fahrzeugs gegeniiber der statischen 
Radlastverteilung um den gleichen Betrag verandert, wahrend der Luftwiderstand eine Belastungs- 
anderung bewirkt, deren Betrag fiir ein Radpaar entgegengesetzt gleich dem fiir das andere ist. 
Damit kann jede im Windkanal gemessene Radlastverteilung in die Rechnung eingefiihrt werden. 
Vorausgesetzt ist allerdings, daB sich der Charakter der Umstrémung des Fahrzeugs nicht mit der 
Fahrgeschwindigkeit andert. Es mu8 nun ahnlich wie bei der Beriicksichtigung des Luftwider- 
standes die durch den Auftrieb bewirkte bezogene Radlastanderung £*(v) in die Gleichungen (24) 
und die zugehérige Bewegungsgleichung der Lenkung eingefiihrt werden. 

Um den Einflu8 des Auftriebs auf die Stabilitat ohne weitere Rechnung zu zeigen, soll hier der 
Luftwiderstand vernachlassigt werden. Dann geniigt es, in den Gleichungen (la) und in der zu- 
gehérigen Bewegungsgleichung der Lenkung (13) die ganzen rechten Seiten, die von den Seiten- 
kraften der Rader stammen, mit dem Faktor 1 — &*(v) < 1 zu multiplizieren. Statt dessen kann 
man sich auch die Masse, bzw. die Drehmassen auf den linken Seiten entsprechend vergréBert denken. 
In (20) wachsen damit der ,,Tragheitswert“ C und die Drehmasse des Vorderrades um den Lenk- 
zapfen I, im selben Verhaltnis, doch wird dadurch der Verlauf der Stabilitatsgrenze auch bei hohen 
Geschwindigkeiten kaum verandert. Da der Auftrieb die Bodenhaftung und damit die Stabilitat 
des Kraftwagens offenbar verringert, scheint das zunachst ein zweifelhaftes Ergebnis. Es muB hier 
jedoch daran erinnert werden, daB die Lenkstabilitat des Kraftwagens nicht nur von der Lage des 
betreffenden Zustands im stabililen Bereich der n — v Ebene sondern auch von der Kleinheit der 
durch die Stérung ausgelésten Schiebe- und Gierschwingungen abhiangt. In (21) bewirkt die schein- 
bare VergréBerung der ,,horizontalen“ Massen bzw. Drehmassen durch den Auftrieb eine Ver- 
minderung der Eigenfrequenz y an der Stabilitatsgrenze mit zunehmender Fahrgeschwindigkeit v; 
wenn der Auftrieb den Wert des Gewichts erreicht, wird v = 0; nach (10) und (11) wachsen dann 
Gier- und Schiebeamplituden ins unendliche. 

Bei einem Kraftwagen kann dieser Zustand nicht erreicht werden, da beim Abheben aller Rader 
keine Schubkraft zur Uberwindung des Luftwiderstandes verfiigbar ware. AuBerdem geht die 
Stabilitat infolge Entlastung der Vorderrader durch den Luftwiderstand schon bei viel niedrigerer 
Geschwindigkeit verloren. Aus dem obigen geht neuerlich hervor, daB nur der Teil der Stabilitats- 
grenze praktische Bedeutung hat, auf dem die Eigenfrequenz nicht unter einen bestimmten Wert 
sinkt. Der andere Teil stellt zwar ebenfalls eine Folge von Zustanden ungedampfter Schwingungen 
dar, doch nehmen dabei Schiebe- und Gieramplituden erhebliche Werte an. Da gewéhnlich eine 
ziemlich schmale F'ahrbahn eingehalten werden muf, kann dieser Teil nicht die Grenze eines prak- 
tisch stabilen Bereiches bilden. Dazu kommt, da8 fir groBe Gier-, Schwimm- und Schraglauf- 
winkel die getroffenen Voraussetzungen — Seitenkraft proportional Schraglaufwinkel, seitliche 
Komponente des Luftwiderstandes vernachlassigbar klein — nicht mehr gelten. In Abb. 7 und 8 
kommt ein solcher ,,ungiltiger“ Teil der Stabilitatsgrenze nur bei sehr niedriger Fahrgeschwindig- 
keit vor und ist dort praktisch uninteressant. Nur bei einem Flugzeug, dessen Auftrieb gro8 gegen- 
tuber dem Luftwiderstand ist, wird er bei hohen Geschwindigkeiten auftreten. Bei diesem iiber- 
nehmen schon lange vor dem Abheben vom Boden Heckflosse und Seitenruder die Stabilisierung. 

SchlieBlich soll nochmals auf die Reihenfolge der Wirkungen der Luftkrafte bei Windstille hin- 
gewiesen werden. Bereits vor Eintritt einer Stérung vermindern sie die Radlasten zumindest eines 
Radpaares und machen damit das Fahrzeug fiir die Stérung anfallig. Schon eine kleine Stérung 
kann dann infolge der verminderten Seitenkrafte der Rader zu niederfrequenten Schiebe- und Gier- 


i 
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schwingungen mit gréBeren Amplituden fiihren. Nur in diesem — recht gefahrlichen — Fall wird 
die Anfachung durch die seitliche Komponente des Luftwiderstandes, bzw. durch sein Moment 
um die Hochachse des Fahrzeugs, wirksam. 


15. Zusammenfassung. Bei der Definition der Lenkstabilitat und der Berechnung der statio- 
nairen Amplituden wurde gezeigt, daB nicht so sehr die Dampfung der Lenk-, Schiebe- und Gier- 
schwingungen eines Fahrzeugs wie geniigend hohe Eigenfrequenz des zugehirigen Systems fiir ein 
stabiles Verhalten absolut notwendig sind. Es wurde daher nicht nur die Stabilitatsgrenze errech- 
net, sondern auch die Verteilung der Eigenfrequenzen im stabilen Bereich untersucht und zwar zu- 
nichst ohne und dann mit Beriicksichtigung des Luftwiderstandes. Abb. 7 und 8 zeigen die Ergeb- 
nisse in der (n, v)-Ebene. Dort hat die Stabilitatsgrenze die Gleichung (20) bzw. (29). In (29) kann 
fiir héhere Geschwindigkeiten das zweite Glied in der eckigen Klammer gegeniiber dem ersten ver- 
nachlassigt werden, es kommt dann mit C = 1/k)g, k = k, G/4 


16 Ip 


n= eset (29a) 
worin 
2h 
C= Gp he 


mit v als gréBter erreichbarer Geschwindigkeit ist. Obwohl n hier der Nachlauf an der Stabilitats- 
grenze ist, diirfte es geniigen, diesen Wert durch entsprechende Riickwartsneigung der Lenkzapfen 
auszufiihren. Der wirksame Nachlauf ist gréBer, da nach H. Fromm die Seitenkraft in der Drauf- 
sicht eteas hinter der Radachse angreift. Die oben angefiihrten Griinde sprechen gegen eine erheb- 
liche VergréBerung des Nachlaufs tiber den Wert nach (29a) hinaus. 

Der Beiwert ky, d.i. die Seitensteifigkeit der Reifen kommt zwar in (29a) nicht vor, wohl aber 
in der Gleichung (30) der ungedampften Eigenfrequenz an der Stabilitatsgrenze. Wegen C = 1/k,g 
sieht man, daB groBe Seitensteifigkeit der Reifen diese Eigenfrequenz erhéht und damit fiir die 
Lenkstabilitat giinstig ist. Schwieriger ist es, den Einflu8 des Radstandes I zu beurteilen. Ein 
groBer Radstand vergréBert zwar den stabilen Bereich in der (n, v)-Ebene und verringert die schad- 
liche Anderung der Radlastverteilung infolge des Luftwiderstandes. Andererseits vermindert er die 
Eigenfrequenz an der Stabilitatsgrenze und begiinstigt so das Entstehen von gefabrlichen nieder- 
frequenten Schiebe- und Gierschwingungen mit groBen Amplituden. 

Bei Steigerung der Fahrgeschwindigkeit mit einem gegebenen Fahrzeug (n = konst.) kommt man 
zunichst noch im stabilen Bereich in Zustande immer geringerer Dampfung, wahrend die Eigen- 
frequenzen des Fahrzeugs sich kaum andern. Beim Uberschreiten der Stabilitatsgrenze ist, wie vor- 
her die Dampfung, danach die Anfachung nur sehr gering, und es wachsen deswegen die Amplituden 
einer durch eine kleine Stérung verursachten Schiebe- und Gierschwingung nur langsam an. Die 
Stabilitatsgrenze in Abb. 8 wird von der waagrechten Geraden n = konst. in einem spitzen Winkel 
geschnitten, so daB mit wachsendem v ein allmahlicher Ubergang aus dem stabilen in den unstabilen 
Bereich stattfindet. Ein aufmerksamer Lenker sollte die Gefahr rechtzeitig erkennen, bevor er 
einer gréBeren Stérung begegnet, wenn z. B. das Fahrzeug von einem seitlichen WindstoB getroffen 
wird. 

Das gewahlte Ersatzsystem fiir das Fahrzeug beriicksichtigt zwar den KinfluB des Schraglaufes 
der Luftreifen bei seitlicher Belastung und die Lenkbarkeit des Fahrzeugs, aber nicht den Einflu8 
der Federung. Eine folgende Arbeit soll diese Liicke ausfiillen. 


(Eingegangen am 16. Mai 1958.) 
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Der durchlaufende Balken auf iquidistanten elastischen Stiitzen* 


Von V. Bogunovie 


1. Einleitung. In der vorliegenden Arbeit wird das Biegungsproblem eines durchlaufenden Bal- 
kens auf elastisch senkbaren Stiitzen untersucht. Dabei wird angenommen, da die beiden End- 
stiitzen starr und simtliche Zwischenstiitzen elastisch, gleich steif und in gleichen Abstanden ange- | 
ordnet sind. Der Durchlaufbalken wird entweder durch Einzellasten, die iiber elastische Stiitzen | 
wirken, oder durch verteilte Last beansprucht. Die Lésung dieses Problems, d. h. die Ermittlung 
der Stiitzkrafte und Stiitzmomente, wird in Form eines trigonometrischen Polynoms angegeben, 
dessen Gliederzahl mit der Anzahl der elastischen Zwischenstiitzen tibereinstimmt. 


2. Elastisch gestiitzter Durchlauftriger durch Einzelkrafte belastet (Abb. 2). 

a) Durchbiegungen und Biegemomente an den Aquidistanten Stellen eines frei | 
aufliegenden Balkens. Wir betrachten einen frei gelagerten Balken von der Stiitzweite I, der 
durch m — 1 aquidistante irgendwelche Einzellasten beansprucht wird (Abb. 1). Bezeichnen P;, u; 
und M, die Einzellast, die Durchbiegung und das Biegemoment an einer Tragerstelle i im Abstand 
i I/m von der linken Stiitze, so lassen sich! mit der Abkiirzung 


[8 
fe = EI 
die Differenzengleichungen fiir die Durchbiegungen u, 
und Biegemomente M, 


8 
u,=— +B (Mia, +4M;+M;_;), 
(2k ott) ; (1) 
l 
ali] by A Pe (2) 


E id on . 


Abb, 1. Balken auf zwei Stiitzen mit Einzellasten. 
mit den entsprechenden Randbedingungen 
uy — U0; By Phy Bad Ue 
anschreiben, wobei J das unverdnderlich vorausgesetzte Tragheitsmoment des Balkenquerschnitts 


und & den Elastizitatsmodul des Materials bedeuten. 
Dabei bezeichnet der Differenzenoperator 6 iiblicherweise die zentrale Differenz 


Of: = heey Lee? 
t+ 
und entsprechende héhere Differenzen 
Of; = fii — 2h + fiers 
Oa) Oe 
wo f, den Funktionswert f; = f(«;) einer Funktion f(x) an aquidistanten Stellen x; = i 1/m bedeutet. 
Aus (1) und (2) 14Bt sich leicht die Differenzengleichung vierter Ordnung 


8 9 
Ont, a Jp (Pig ew Py yey) (§@ =1,2,...,m—1) (3) 
ableiten. 
Die Werte von aquidistanten Einzellasten P; kénnen, wie bekannt?, in Form eines trigonometri- 


schen Polynoms 
m—1 


P= Yon hay hei) (4) 
pai 
dargestellt werden, wobei die Koeffizienten durch die Summe 
m—l1 
Gg == >) Prinika (5) 
t=1 


gegeben sind. 


* Herrn Prof. J. M. Klitchieff in Beograd zum 70. Geburtstag gewidmet. 
1 F, Stiissi, Vorlesungen uber Baustatik, Bd. I, S. 55 u. 238. Basel 1946. 
® R. Courant, Vorlesungen iiber Differential- und Integralrechnung, Bd. I, S. 360. Berlin 1930. 
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Auf ahnliche Weise kénnen auch die Werte M, der Biegemomente an den aquidistanten Trager- 
stellen durch das trigonometrische Polynom 
m—l1 


M,= ¥ bsinika, (6) 
k=1 


largestellt werden, wo die Koeffizienten 6, vorlaufig noch unbekannt sind. 
Wenn wir die Ausdriicke (6) und (4) in die Differenzengleichung (2) einsetzen und beachten, daB 


Of, = — 2 (1—coskA) sinika (7) 


Z 


st, wobei wir die Bezeichnung 
f, =sinika 
infiihren, so erhalten wir ein System von m — 1 linearen Gleichungen fiir die m — 1 Unbekannten 


a...» b 


2 On_1* 


m—l m—1 
ns l ro ‘ 
a og eri ¢, sintk A (t= 1,2,...,m—1). (8) 


Dieses Gleichungssystem kann, wie bekannt, sehr leicht explizit aufgelést werden. Wir wahlen 
rgend eine ganze Zahl p unter den Zahlen 1,2,..., m—1, multiplizieren die erste Gleichung mit 
sin 1 p A, die zweite mit sin 2 p/,..., die letzte mit sin (m—1) p A, addieren sodann saimtliche Glei- 
chungen und ordnen die so entstandene Gleichung nach den Koeffizienten b, bzw. c,. Unter Be- 
niitzung der bekannten Formel 


m—l1 Bi 
fi = 
> sintkAsini pA = me sails a 


i=l 


(9) 


ergibt sich, daf alle Beiwerte von b, bzw. c, fiir k += p verschwinden bis auf k = p. Auf diese Weise 
erhalten wir 


b, = 5 (1—cos kA) gy. (10) 


Setzt man in (6) den Ausdruck (10) fiir b, ein, so erhalt man die Werte M; der Biegemomente an 
len aquidistanten Stellen eines durch Einzellasten P; beanspruchten Balkens 


m—l1 
Ay bse 
M; imal (l1—coskA)“* c, sinikd. (11) 


Auf ahnliche Weise lassen sich die Verschiebungen u; an den aquidistanten Tragerstellen durch 
las trigonometrische Polynom 


m—l 
Uj= >) 6, snk 1, (12) 
k=1 


larstellen, wo die Koeffizienten a, vorlaufig unbestimmt sind. Wenn wir die Ausdriicke (12) und (4) 
n die Differenzengleichung (3) einsetzen und beachten, daB 

Of, =4(1—coskdsinika, fir t4fitfii1 =2 (2 + coskd)sinikd (13) 
st, so erhalten wir ein Gleichungssystem von m — 1] linearen Gleichungen fiir die m — 1 Unbekann- 
MD O,, Bo, 0 «+5 On 4 


m—l1 m—1 


4 a,(1—cos kA) sinikA=—%fg >) o (2 +coshkA)sinikA ((=1,2,...,m—1). 
=1 


k=1 


Dieses Gleichungssystem kann auf dieselbe Weise wie das System (8) explizit gelést werden; so 
rgibt sich 


4 2-+ cosk A 14 

Opes JE oak sb ( ) 

Setzt man in (12) den Ausdruck (14) fiir a, ein, so erhalt man die Verschiebungen u; an den aqui- 
listanten Stellen eines durch Einzellasten P; beanspruchten Tragers 


m—1 
4 2+ coska hag 15 
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b) Balken auf aquidistanten elastischen Stiitzen. Esseiein an beiden Enden senkrechi | 
unverschieblich gelagerter Balken von der Lange! gegeben, der in den Aquidistanten Punkten 1,2,.. 1 
m — 1 elastisch gelagert und durch die Einzellasten P; beansprucht ist (Abb. 2), wobei lineare> 
Zusammenhang zwischen Verschiebung u; und Stiitzent 
kraft R; vorausgesetzt sei, derart, daB | 


u; =f, R; (LOE 


ist und mit f, die Verschiebung der elastischen Stiitzy) 
infolge R; = 1 bezeichnet wird. | 

Zur Bildung eines statisch bestimmten Systems be: 
seitigen wir die m—1 mittleren elastischen Stiitzenr 
In diesen Punkten 1 bis m — 1 bringen wir die statisck) 
unbestimmten Stiitzenwiderstande R; als Belastung any 
die mit den vorgegebenen Einzellasten P; die resull] 
Abb, 2. Elastisch gestiitzter Durchlauftrager tierenden Ejinzelkrafte | 


mit Einzellasten. | 


ergeben. Die vorgegebenen GréfSen der Einzellasten P; kénnen wir in Form (4) eines trigonometrifi 
schen Polynoms darstellen, dessen Koeffizienten durch (5) bestimmt sind. Auf ahnliche Weise lassen} 
sich auch die Stiitzenwiderstande R; in Form | 


m—1 
R,= DG, sint ks (183 


= | 


ausdricken, wo die Koeffizienten ¢,, vorlaufig noch unbekannt sind. Setzt man in (16) den Ausdrugl 


(18) fiir R; ein, so erhalt man | 


m—1 
Uf, ac. si Ui An (19) 
k=1 | 
Nach Einsetzen der Ausdriicke (4) und (18) in (17) erhalt man 
m—1 | 
0; = Dd (¢,— cy) sini kd. (20F 
k=1 


Setzt man nun die Ausdriicke (19) und (20) in die Differenzengleichung (3) ein, wobei an Stell 
von P; der Ausdruck (20) fiir Q; einzusetzen ist, so erhalt man ein System von m — 1 lineare 
Gleichungen mit m — | Unbekannten ¢,, 


m—1 m—l1 | 


A>: Cet (l—cosk A sinika=—° fy >) (¢, — ¢ 4) (2 + coskA) snikdA 
bel b= 


G=1,2, ... 7 
Lést man dieses Gleichungssystem auf dieselbe Weise wie das System (8), so erhalt man 
m® f, (1 —cosk A)?|-1 
as 4 fp 2+ coska | “ 
Setzt man in (18) den Ausdruck (21) fiir ¢,, ein, so ergibt sich 


m—l1 
3 fs (l—cosk A)? 
U pe; 4 fg 2+ coska 


k=1 


Nun sind die resultierenden Einzelkrafte Q; durch die Gleichungen (20), (5) und (21) bestimmt} 
und es ist leicht, die Stiitzmomente M; mit Hilfe (11) zu ermitteln, wobei an Stelle von c, die Dif+ 
ferenz c, — c,, eingesetzt werden soll. Auf diese Weise ergibt sich 


m—1 
Me a 


ai = 


(21) 


1+ 


I. csinikA. (22), 


4 fe2 “1 } 
1—coskd 4 mi pie c, sintk). (23)) 


c) Beispiele. 


1. Eine Einzellast iiber eine elastische Stiitze. Mit Riicksicht auf (9) ist leicht zu erkennen, daB deri 
Ausdruck 


m—1 
9 i) 
Diesen > sink pdAsinkid (24)} 
ta 
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fiir alle Tragerstellen i  p verschwindet und nur an der Stelle i = p den Wert P, annimmt und 
somit die vorgegebene Belastung (Abb. 3) in den diskreten Punkten 1, 2,...,m— 1 darstellt. 
Dementsprechend kénnen wir schreiben 


m—1 
P=) o emik/, 
k=l 
wo 
2 : 
tes P, sink pa (25) 


ist. Setzt man (25) in die Gleichungen (22) und (23) ein, so folgt 


m—l1 
: cae, m3 fs (l—cosk A)? 
entre, L | 4 fg 2+ coska 


if 


sinpkdAsinikd, 


—1l1 


Mi, = te 


es 


4 fe 2+ coska 
m3 f, 1—coskA 


B 
eee - is ~Sn-1 
ae 


[1 — eos ka + | sink pAsinika, 


7 L m-1 
Le ae 
cae So areal 
Abb. 3. Elastisch gestiitzter Durchlauftrager mit Einzellast. Abb. 4. Sinusférmige Verteilung der Einzellasten, 


wobei die Indexbezeichnung so durchgefiihrt ist, daB der erste Index den Ort, der zweite die Ursache 
angibt. 


2. Gleiche Einzellasten an sémilichen Auflagerpunkten der elastischen Zwischenstiitzen. Unter Be- 
riicksichtigung, daB 


Pi = P2=::-=P,,_,=P 
ist, folgt mit Hilfe der Formel 
m—1 
a sini k A = sin? we ctg ee 
aus (5) 
c, = ae sin® a ctg iM (26) 


Setzt man (26) in die Gleichungen (22) und (23) ein, so erhalt man 


io me: fe (la cosk Alay kA 0 Baia 
R= 7, PDI Dy SER Eg) ctg = sin” sintkd, 


m—l1 
PH 4 =a [BGG a 6 
M;=—5 > (1 — cos k A)? we 72 + 08k] sin k Asin? sini kA. 
k=1L 


3. Sinusformige Verteilung der Einzellasten. Die GréBe P; der Kinzellast an der Tragerstelle 1 


(Abb. 4) sei durch den Ausdruck 
P= Pisin tn) (27) 


gegeben, wo n irgendeine nicht durch m teilbare ganze Zahl bedeutet. Die vorgegebene Zahl n 
bestimmt eindeutig eine ganze Zahl j (> 0) durch die Bedingung 


2jm<an<2(7+1)m. 
Mit einer so bestimmten Zahl j dirfen wir 
n=2jm+y fiir n<(2j+1)m, (a) 


n=2(j+1)m—y fiir n>(2j+1)m, (b) 
g* 
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in der Form | 
P,;=+ Psiniva (283, 
schreiben; hier gilt positives Vorzeichen fiir den Fall (a) und negatives fiir (b). Aus den Gleichungex | 
(22) und (23) und mit Riicksicht auf (28) und (4) erhalt man ] 


3 athe a a | 
nak Piya 


-1 
M, = Ly, [!— 00894 + us aes siniva. 


m> f, 1—cosvA 


3. Balken auf elastischen Stiitzen mit kontinuierlicher Belastung. a) Durchbiegungen unc( 
Biegemomente an den Aquidistanten Stellen eines frei gestiitzten Balkens. — Did 
Durchbiegung w eines freiaufliegenden Balkens von der Lange Ji 


2 
96 ee Aaa 8 (Bee? | 
u=——Ofp > ja sin —— sin (29), 


fs gS 
NEES ceca the ees 
l dar 


hie gestellt werden. Aus dieser Grundlésung ergibt sich sofort dis 

Abb. 5. Balken auf zwei Stiitzen mit 2 ee. fs 5 ae 4 
kontinuierlicher Belastung. Lésung fiir jede beliebig verteilte Belastung, derer Intensitat q eine 
Funktion von é ist (Abb. 5). Dafiir mu8 man in den Ausdruck! 
(29) Q = qdé einsetzen und iiber die belastete Balkenstrecke integrieren. Auf diese Weise erhal1 
man 


6, < i | 
rere a pm he _ (Sy 
mit Abkiirzung 
1 
I = 1B 
An = | ©) sin 7 ae (31) 
0 
Wenn wir 
fa eee a 
m 


% , © 

u=—; Se >) Ansinika. (32) 

k=l 

Die Durchbiegungen u; an den Aquidistanten Tragerstellen lassen sich in Form eines trigonometri+ 

schen Polynoms | 
m—1 

96 as : 

us =—7 fe >} a, sini kA (33) 

k=1 | 


darstellen, wo die Koeffizienten a, vorlaufig noch unbestimmt sind. Wenn wir die Ausdriicke (32), 


und (33) gleichsetzen 
} ApsnikA= J’ asinikA, (i =1,2,...,m—1), (34) 
k=1 k=l 


so erhalten wir ein System von m— 1 linearen Gleichungen mit m — 1 Unbekannten Gy, Ag, +0050, 4 
Dieses Gleichungssystem kann auf dieselbe Weise wie das System (8) explizit aufgelést werden. Wi 

wahlen irgend eine ganze Zahl p unter den Zahlen 1, 2,..., m—1, multiplizieren die erste Gleichung; 
mit sin 1 pA, die zweite mit sin2 p2,..., die letzte mit sin (m — 1) pA, addieren sodann samt+ 
liche Gleichungen und ordnen die so entstandene Gleichung nach den Koeffizienten 4, baw. Oy, 
Unter Benutzung der bekannten Formel 


<> m/2 fiir k=2jm+p, | = 02 
>) sinikdAsinips= — m/2 fiir k=2jm—p, Shae ee tin 


| 
| 
} 
| 


i=l 0 fiir alle iibrigen Werte von k, 
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und mit Hilfe der Beziehung 
sin(2jm+i)kA=4+simikaA (j =0,1,2,...) 


ergibt sich sodann, daB samtliche Beiwerte von A, bz fii ] i 
A r¢ , baw. a, fir k ~ 27m + pbaw.k verschwin- 
den bis auf k = 2j m+ p baw. k = p. Auf diese Weise erhalten wir ‘i aes ng 


ON, = 2 Azim ~ O Atimaw (35) 
Unter der Voraussetzung, daB die Koeffizienten A, die Gleichung 
Ee ahs ry eae a Se 
2jmitk 1yr+1 (2j m -4 k)r A, (36) 
erfiillen, woj = 1,2,...;k =1,2,...,m—1;r=4,5,... sind, kénnen wir dann den Ausdruck 
(35) in der Form 
o = A 
schreiben, wobei die Abkiirzung nares ve 
Per =k SS) Qim+hy- (38) 
j=—o 


eingefiihrt ist. Die Reihe (38) 1a8t sich geschlossen auswerten, wie in Ziff. 4 gezeigt wird. Indem 
man (37) in (33) einsetzt, erhalt man schlieBlich den Ausdruck fiir die Durchbiegung 


m—1 
96 ee 
u, = — fa 2, Q,, 4, sniki. (39) 


Differentiiert man die Reihe (29) zweimal nach x, so erhalt man den Ausdruck fiir das Biege- 
moment 


9 foo} 
M=—EIu! =+ SB, sin ** , 
a 


oder an den aquidistanten Punkten 


21e eee 
M, ==5 >) Bysini kA, (40) 
k=1 

wo 
B, =K? A, (41) 
ist. Mit Riicksicht auf (41) und (36) erhalt man 
- bet 
Bajmsk = (_qy+1 Q@jm Eh? ME a 


Auf ahnliche Weise laBt sich auch das Biegemoment M, an den dquidistanten Tragerstellen in 
Form eines trigonometrischen Polynoms 


m—l1 
2 
M, Bet > Prr—2 B, sini kd 
k—1 


7 


ausdriicken. Mit Riicksicht auf (41) ergibt sich daraus 


= 
a's nike 
M; = ake > Pr,r—2 k2 A, sin U ka . (43) iy —| 
b=1 1 | 
Aas i Abb, 6. Elastisch gestiitzter Durchlauftrager mit 
b) Balken auf Aquidistanten elastischen kontinuierlicher Belastung. 


Stiitzen. Es sei ein an beiden Enden senkrecht un- 

verschieblich gelagerter Balken von der Lange | gegeben, der in den aquidistanten Punkten 1, 
2,...,i,...,m— 1 elastisch gelagert und durch eine kontinuierliche Belastung von der Intensitat 
4(€) beansprucht wird (Abb. 6). Wenn wir mit u; die Verschiebung der elastischen Stiitze i, mit R; 
den Stiitzdruck und mit f, die Verschiebung dieser Stiitze infolge R; = 1 bezeichnen, dann gilt die 
Beziehung (16). Die unbekannten Stiitzedriicke R;, lassen sich in Form eines trigonometrischen Poly- 
inoms (18) darstellen, dessen Koeffizienten vorlaufig unbestimmt sind. Aus (16) und (18) erhalt 


‘man die Verschiebung u,; (19). 
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Die Durchbiegung u;,) eines Balkens auf zwei Stiitzen an der Tragerstelle i, d. h. im Abstand 

i 1/m vom linken Auflager, ist infolge einer kontinuierlichen Belastung von der Intensitat q(&) dure) 

das Polynom (39) dargestellt. Die Durchbiegung u;, eines Balkens auf zwei Stiitzen an derselbes) 

Tragerstelle i ist infolge von vorlaufig noch unbekannten aquidistanten Stiitzkraften R; durch da 

Polynom (15) ausgedriickt, wobei der Koeffizient ¢,, an Stelle von c, zu setzen ist. 
Diese Verschiebungen miissen offensichtlich das Gleichungssystem 

Ug — Uz, = Uy; (i Leet (444) 

erfiillen. I 

Setzt man nun die Ausdriicke (19), (15) und (39) fiir die Verschiebungen in (44) ein, sq) 

ergibt sich ein lineares Gleichungssystem von m— 1 Gleichungen fiir die m — 1 Unbekannten) 


Cyy> Cora + + +o Om—a1y1 | 
m—l1 m—1 } 
4 2) ki Ais 96 at : 

> \f em te Ga eas |e siniki=— fs >) Py, A, Sint ki, (i = 1,2,...,m—l)) 
k=1 k=1 | 
(455 


Lést man dieses Gleichungssystem auf dieselbe Weise wie das System (8), so erhalt man 
_ 24 m3 [ m?> fs 2+ coskd |-1 uM 
mw | 4 fe ' (l—cosk mA Vhrecke 


Setzt man in (18) den Ausdruck (46) fiir c,, ein, so erhalt man 


Chi 


m—l1 


24m’ ee kaj at Y 
R=—~ >} E es Tea Pur A, Sint ka. (477 


at 
i 


Ist M;, das Moment an der Tragerstelle i infolge der m — 1 Stiitzkrafte R;, und setzen wir fii 
c, den Ausdruck (46) in (11) ein, so ergibt sich 


ae deity = 
Mya ee? Peleus “Ope Ap sini kA. (48) 
ee at | 


Ist M;, das Moment an der Stelle i infolge der gegebenen Belastung q(é) (43), so erhalt man dan}, 
das endgiiltige Stiitzmoment durch Superposition 


M; = Mio— Mix 
ay DSS wy m> fs 1 ki 2+ coskA|“ , 
ogee —, Sma Phr-2— [Ta gpa at Terearka Pky 


wobei zur Abkiirzung 


kt A, sinikd, (499) 


| Phe = Pn kt (50) 
gesetzt ist. | 

c) Der durchlaufende Balken auf aquidistanten starren Stiitzen. Setzt man f, = {| 
in (47) und (49) ein, so ergeben sich die Ausdriicke fiir die inneren Auflagerkrafte R, und die Stiitzy 
momente M; eines durchlaufenden Balkens auf aquidistanten starren Stiitzen 


m—1 

__ 24m8 (1 — cos k A)? Ra 

R, =| a ek Yer Apsinik A, (51) 
= } 


m—l1 
1192 2 2 oe 
hea Sy io We eee IG 7p (52) 
kw | 


GE ny Deora 


d) Die Belastungsfunktionen q(§), welche die Bedingung (36) erfillen. Fiir einiget 
in der technischen Anwendung sehr oft vorkommende Belastungsfialle, ist die Bedingung (36) erfiillt} 
Nehmen wir, z. B. an, daB die Belastungsfunktion in der Form eines Polynoms 


q(S) = qo P.(€) = % (ap + 4, + a, &% +++. + a, én) (53} 
gegeben sei. Nach Einsetzen des Ausdruckes (53) in (31), erhalt man durch partielle Integratio 1) 


A, = © {[P(0) — P(t) cos kx] — (Ez) [Pn(0) — Pa(l) cos kat] + +++, (54) 


Mit Riicksicht auf (54) ergibt sich, daB fiir n < 2, d.h. fiir die gleichmafige (n = 0) und fi " 
die hydrostatische Belastung (n = 1) die Bedingung (36) offensichtlich erfiillt ist. Auch wenn 
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n = 2 ist — solche Belastungsfalle kommen in der technischen Praxis sehr selten vor — kann die 
Bedingung (36) erfiillt werden, z. B., wenn die Koeffizienten des Polynoms die Gleichungen 


PPO) =0, PRI) = 0 


fiir 
p=0,2,4,:..,4—2, p=n—s[l—(—1)" 
erfiillen. 


Belastung g(E) Ax 
2ql 2k 
pm [oye 


% = as sinkpA sinksA 


17 
= ps 00s (k+1) 0 


& 
~|Ne 


Fiir eine Streckenlast gq), gemaB Tabelle (Abb. b), wo p und s ganze Zablen sind, ergibt sich aus 
(G1) 


A; = “aol “sink pAsinks A. 


Dieser Koeffizient erfiillt offensichtlich die Bedingung (36). In Tabelle sind einige in der technischen 
Praxis oft vorkommende Belastungsfalle angegeben, deren Koeffizienten A, die Bedingung (36) 
erfiillen. 


e) Beispiele. 

1. Der Durchlauftrager auf elastischen Stiitzen mit den in Tabelle dargestellten Belastungen. 
Setzt man in (47) und (49) den Ausdruck fiir yf, (r = 5) (s. Ziff. 4) ein, so lassen sich mit der Ab- 
kurzung 
ait ip 
OE fp 
die Ausdriicke fiir die Stiitzdriicke R; und die Biegemomente M, in der Form 


m—l1 


Py ; 5+ coska 
Rs =m 2 : An yee icak A A oak 


ctg ys sinik), 
m—1 

= 2 5 26 (1—coskA)—1 
LE Sina 2 ep 2+ coskA+ B(1—cosk A)? 
schreiben, wobei die Koeffizienten A, in Tabelle angegeben sind. Setzt man 6 = 0 in obige Aus~ 
driicke ein, so erhalt man die Stiitzkrafte und Stiitzmomente eines durchlaufenden Tragers auf 


starren Stiitzen 


ctg = sntikA 


m—l1 , 


_ 2 54 .tecoska kA. 
a EEN ac gp My cm 
ki 
m—1 Se 
= ae 5 eee pee ; 
Mam A, > cos ka 80 t kA 


| 2. Sinusférmig verteilte Belastung (Abb. 7). Fiir den betrachteten Fall ist die Belastungsfunktion 
q(é) durch den Ausdruck 


5, Wan 
Tae 6 OU = 


‘i 
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gegeben, wobei n irgendeine nicht durch m teilbare ganze Zahl bedeutet. Der erase fiir dii 
Verschiebung u;, eines Balkens auf zwei Stiitzen lautet infolge dieser Belastung bekannterweise | 


oes oe Yin = fod lsinind. 


wai m4 
Wenn wir das gleiche Verfahren wie im Beispiel 3 vor) | 
Ziff. 2 anwenden, ergibt sich i 


U;9 = + feqol sini vA. 


7 2 i m-1 Setzt man diesen Ausdruck in (44) ein, so erhalt mal 


N 
te tae + nach den im Ziff. 3b) beschriebenen Verfahren die Auss) 
Ll -| driicke fiir die Stiitzkrafte R; und Biegemomente M/ 
Abb. 7. Sinusférmig verteilte Belastung. 12 m? 2+ cosva 


R,=+ ae 16 ! meron sintvd, 


6 m? n? 7? 2+ cosvd cae 
Mi, = gato eas F |B (Leos) +5 ert sini, 


wobei positives Vorzeichen fir den Fall (a) und negatives fiir (b) gilt, [s. Beispiel 3 von Ziff. 2] . 
4, Geschlossene Auswertung der Reihensumme!. Durch gliederweise (r — 1)-malige Differenj 
tiation der bekannten Reihe | 
foe) 
etgx = J) (jm + x)" 
j=—o 


ergibt sich 


(org a9 = (— 47) Y (ja tay. (55 


j=—o | 


Setzt man x = k7/2m in (55) ein, so erhalt man nach einigen Umformungen, folgende, hier bis 
r = 5 angefiihrte, geschlossene Ausdriicke 


Re es ff 4 

Git Sane: oa ae J@ + c08 & 4) (1—coskA)°, 
eee rs 2 kh 
GV. = (=) (1—coskA)*, 3 =- ae ) (5 + cos kA) (l—cosk A)” ctg > | 

ha =2(> JQ cos k ay et os 


fiir die Reihe 
= 5 Qim+h. 


j=—o 


(Eingegangen am 16. Mai 1958.) 
Anschrift des Verfassers: Dr.-Ing. V. Bogunovic, Beograd (Jugoslawien), Palmotiéeva 31. 


1 J. M. Klitchieff, The Aeronautical Quarterly, 2, 1950, S. 157. 
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Torsionseigenfrequenzen von Kreisscheiben veranderlicher Dicke 


Von H. Wittmeyer 


1. Einleitung. Die Ermittlung der Torsionseigenfrequenzen von Kreisscheiben hat Bedeutung 
bei dem Entwurf von Dampfturbinenscheiben, wie Grammel! ausfiihrlich dargestellt hat. Dort 
werden auch Liésungsmiglichkeiten fiir die zugehérige Differentialgleichung angegeben. 

Nun kann die Differentialgleichung fiir die Torsionsschwingungen einer Kreisscheibe ohne 
Knotendurchmesser auch als Differentialgleichung fiir die Torsionsschwingungen eines Stabes 
aufgefaBt werden, wie weiter unten gezeigt wird. Daher liegt es nahe, die einfachen Formeln, 
die der Verfasser® in einer friiheren Arbeit zur Berechnung der Torsionseigenfrequenzen eines 
Stabes hergeleitet hat, auf die Berechnung der Torsionseigenfrequenzen einer Kreisscheibe zu 
tbertragen. Die auf diese Weise erhaltene spater angefiihrte Grundformel erméglicht in einfacher 
Weise die Berechnung der ersten drei Torsionseigenfrequenzen der Scheibe, falls sie an ihrem 
Innenrande starr eingespannt ist und keine trage Masse am AuBenrande tragt. Mit einer weiteren 
Formel kann berechnet werden, wie sich die Grundeigenfrequenz andert, falls man den Innenrand 
elastisch einspannt und den AuSenrand mit einer tragen Masse (Beschaufelung) belest. 


2. Herleitung der Niherungsformel. Im folgenden seien dieselben Bezeichnungen wie von 
Grammel benutzt, wobei der Index o und der Index a besagen mégen, daf die betreffenden GréBen 
am Innen- bzw. am Auf®enrand der Scheibe zu nehmen sind, namlich r der radielle Abstand eines 
Scheibenpunktes von der Scheibendrehachse; y(r) die Amplitude der sinusférmig von der Zeit 
abhangigen Torsionsschwingung eines Scheibenpunktes im Abstand r von der Scheibendrehachse 
um diese Achse, gemessen als Winkel im Bogenma; y(r) die Scheibendicke; « = a; die k-te Eigen- 
frequenz der Scheibentorsionsschwingung; y das spezifische Gewicht des Scheibenmaterials; g die 


_Erdbeschleunigung; G der Schubmodul des Scheibenmaterials; ¢ die Bettungsziffer, die sich wie 


folgt berechnet: Ist M, ein am Innenrand angreifendes Torsionsmoment, bezogen auf die Dreh- 
achse und #, die Drehung des Innenrandes um die Drehachse infolge des Momentes M,, so ist 


-e€=M,/(221r2 ¥) 0), wie man leicht aus der ,,Technischen Dynamik“, Seite 36, Gleichung (5) 


B. 
2 
achse) der tragen Randmasse (Beschaufelung) je Flacheneinheit des AuBenrandes. Die meist 
kleine dimensionslose Zahl § wird ,,Gewichtsfaktor der Beschaufelung“ genannt (siehe ebendort, 


und (8) folgert. Ferner ist 9 = >y a das Massentragheitsmoment (bezogen auf die Scheibendreh- 


§. 38). 


Vernachlassigt man den EinfluB der Zentrifugalkraft, was nach Grammel (ebendort S. 35) bei 
der Bestimmung der Torsionseigenfrequenzen zulassig ist, so lautet die Differentialgleichung fiir 


die Torsionsschwingungen ohne Knotendurchmesser [ebendort S. 33, Formel (9)] 


: 


/ 3 , 

Wi 4(Z42)¥ +Py=0, 

_wobei Striche Ableitungen nach r bedeuten, mit den Eigenwerten 
Aah, =22m\/% - (2) 

Fir den Innenrand (r = r,) gilt [ebendort S. 36, Formel (5)] 
Gry’ —eyp =0 (3) 
und fiir den AuBenrand (r = r,) [ebendort Formedl (6)] 

Gry’ —407eOy=0. (4) 


Durch Multiplikation mit geeigneten Faktoren lassen sich die Gleichungen (1), (3) und (4) nach 
_Einfithrung von = 470? auf die folgende selbstadjungierte Form bringen: 


(CyrFy) +uZyry=0, (5) 


1 Siche CG: B. Biezeno u. R.Grammel, Technische Dynamik, Bd. 2, S. 29ff., Berlin/Gottingen/Heidelberg 


1953. 


2 H, Wittmeyer, Forsch. Ing.-Wes. 24 (1958) S. 37. 
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—Gyry +teyry=0 (r=), (6). 
Cyr BY TOU) ST aia (7) 
Diese Gleichungen kénnen gema®B meinen friiheren Gleichungen? (1), (3), (8) und (9) als die’ 
Differentialgleichung nebst zugehérigen Randbedingungen fiir die Torsionsschwingungen eines | 
Stabes von der Lange r, — ry) mit der Langskoordinate r — ry aufgefaBt werden, wenn man diesem | 
Stabe die folgenden BaugréBen gibt: G y r° die Torsionssteifigkeit des Stabes, (y/g) y r? das Massen- | 
trigheitsmoment des Stabes je Lingeneinheit, ey) rj die Torsionsfederzahl der elastischen Ein- | 
spannung des Stabes bei r = rp (Moment je Winkel im Bogenmaf), y,r, 0 das Massentrigheits- | 
moment, bezogen auf die Torsionsachse des Stabes, einer Masse am Stabende (r = r,), 4 das Quadrat_ 
einer Torsionseigenkreisfrequenz des Stabes. 


Bei der Anwendung der Formeln aus meiner Arbeit (Abschnitt 3) auf diesen ,,Ersatzstab“, der) 
die gleichen Eigenfrequenzen wie die Scheibe hat, ergeben sich insofern einfache Verhiltnisse, als 
die Verdrehsteifigkeit dem Tragheitsmoment je Langeneinheit proportional ist. Kine Scheibe mit | 
y =konst.- r~® entspricht einem Stabe mit konstantem Tragheitsmoment je Langeneinheit und | 
konstanter Steifigkeit, eine sich weniger nach aufen verjiingende Scheibe entspricht einem Stabe} 
mit nach auBen zunehmendem Tragheitsmoment je Langeneinheit und zunehmender Steifigkeit. | 
Um insbesondere die letztgenannten, praktisch wichtigen Scheiben durch die Naherungsformeln | 
erfassen zu kénnen, wurde der Giiltigkeitsbereich der Formel (68) meiner friiheren Arbeit fiir die} 
Kennwerte a7, =7,(L7) und t, =1,(Lr) angegebenen Formeln entsprechend erweitert. Die. 
einfache Anwendung meiner friiheren Formeln auf den hier definierten Ersatztorsionsstab soll hier; 
nicht im einzelnen ausgefiihrt werden. Das Resultat ist in den beiden nachsten Abschnitten dar- | 


| 


gestellt. 


3. Torsionseigenfrequenzen einer Kreisscheibe mit fester Einspannung am Innenrand und mit 
freiem Aufenrand. Die Dickenverteilung der Scheibe sei durch y = y(r) mit viermal stetig dif- 
ferentiierbarem y(r) gegeben. Ist die Scheibe bei r = ry fest eingespannt und bei r =r, frei, so) 
berechnet man Niaherungswerte fiir ihre drei ersten Torsionseigenfrequenzen a,, 4, und a3 auf/ 

| 


folgende Weise: Mit 
ry =r) + 0,12 (r, — 19), To = 1) + 0,88 (r, — ro) (8)) 
bilde man 
pe es r? y¥(r1) : 
pep = IS V(r) ©) 
und weiter 
Op, = 0,25 + 0,076743 Lp + 0,0057825 L7 — 0,00007882 Lz, (10) 
Op, = 9,75 + 0,025581 Lp + 0,0091235 L% — 0,00010018 L,, (11) 
Ops = 1,25 + 0,015349 L, + 0,0056615 L7 — 0,00002878 L*. (12) 


Ly muB in dem Intervall 
$<IpS17 (13) 
liegen. Dann wird 


1 eG 
oy, FY Ore Ge am (k = 1, 2,3). (14) ) 


Ta) MY 

E alls 
fe re y(ry) rn y (73) 
Gr = lg = [r3, (r, Ay ae 


ist, erhalt man einen im allgemeinen besseren Naherungswert fiir o,, wenn man in (14) fir k = 1| 
die Ersetzung 


ae 1 
SLM fory , tm => (to +7) (14a) 


r 


Ory —> Ory —— Ory == 0,119 Gr (14b) ) 


vornimmt. Andert sich dabei die Frequenz a, wesentlich, so wird von einer Anwendung von (14) 
fir k > 1 und von (18) abgeraten. Ist (14a) nicht erfiillt, so kénnen die Torsionseigenfrequenzen | 
der betreffenden Scheibe mit den Formeln dieser Arbeit im allgemeinen nicht mit brauchbarer > 
Genauigkeit berechnet werden. Dies ist das Ergebnis einer Reihe von Beispielrechnungen. | 


4. Anderung der Grundtorsionseigenfrequenz der Kreisscheibe infolge elastischer Einspannung } 
des Innenrandes und Massenbelegung des Aufenrandes. Der Innenrand sei mit der Bettungs- | 
| 


4 Siehe FuBnote 2 Seite 113. 
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ziffer ¢ elastisch eingespannt und der AuBenrand mit einer Beschaufelung mit dem Gewichtsfaktor 
versehen. Die Tragheitswirkung der Einspannung sei vernachlassigt und die Beschaufelung als 
starr angesehen. Zur Frequenzberechnung ermittle man zunachst unter Benutzung von Ly, nach 


(9) den Wert 


t, = 1 — 0,330938 Lp + 0,0648817 L + 0,00006119 L*,. (15) 
Hiermit und mit 07, nach (10) bestimme man 
eS Ta p Grr, 
Aap, = — hay (ra — 1) 5 “ts 9 | . (16) 
Mit . 
Ory = Oy + Aaqy (17) 


ergibt sich schlieflich fiir die Grundtorsionseigenfrequenz a} der am Innenrand elastisch einge- 
spannten und am AuBenrand mit trager Masse besetzten Scheibe 


* 1 gG 
at Yar, er) ae (18) 
Damit der Fehler von aj héchstens etwa 1% (2%) groBer ist als der von a,, sollte die aus (16) 
zu berechnende GréBe | Aany,/a7, | den Wert 0,085 + 0,0056 Lp (baw. 0,127 + 0,0064 Lp) nicht 
uiberschreiten. Dies ist das Ergebnis einer gré8eren Anzahl von Beispielrechnungen. Hinsichtlich ¢ 
ist diese Bedingung erfiillt, wenn der Wellenhalbmesser nicht wesentlich kleiner als der Innen- 
halbmesser der Scheibe sowie die gréfte Scheibendicke ist, und wenn der kleinste Abstand eines 
-Schwingungsknotens der Welle von der Scheibe klein gegen den Scheibenhalbmesser ist. Sollte die 
_Bedingung wegen eines zu groBen f nicht erfillt sein, so empfiehlt es sich (insbesondere, falls Lp< 0 
he den Einflu8 der Zusatzmasse am AuBenrande mit Hilfe des Dunkerleyschen Verfahrens zu er- 
_fassen. 
) 5. Beispiele. Als erstes Beispiel fiir eine Frequenzberechnung wahlen wir wie bei Grammel 
_(ebendort Seite 38 unten) eine Stahlscheibe mit der Dickenverteilung y = c r~*, wo c eine positive 
_Konstante ist, mit starrer Kinspannung am Innenrand mit dem Halbmesser ry = 10 cm und einer 
Beschaufelung mit dem Gewichtsfaktor 6 = 1/20 am AuBenrand mit dem Halbmesser r, = 40cm 
Wie bei Grammel (ebendort Seite 34) sei /g G/y = 3,26- 105 cm/s. Nach (9), (10), (15) und (17) 
wird Lp = 1, op, = 0,25, t. = 1 und a}, = 0,2417, und damit wird die Grundtorsionseigenfrequenz 
) der Scheibe ohne Beschaufelung nach (14) ungefahr gleich 2717 Hz und mit Beschaufelung nach 
(18) ungefahr gleich 2630 Hz. Diese Werte unterscheiden sich in den angegebenen Ziffern nicht 
von den zugehérigen exakten Werten. 
Als zweites Beispiel sei dieselbe Scheibe mit der Abanderung untersucht, dafi ihr Dickenverlauf 
durch y = cr dargestellt wird, wo ¢ eine positive Konstante ist. Ohne Beschaufelung ergeben 
sich nach (14) folgende Torsionseigenfrequenzen, wobei in Klammern jeweils die relativen Fehler! 
gegeniiber der exakten Lésung angegeben sind: a, ~ 1478 Hz (0,16%), a» ~ 7970 Hz (1,3%), 
G3 13480 Hz (0,5%). Eine Beschaufelung mit dem Gewichtsfaktor # = 1/20 setzt nach (18) 
die Grundtorsionseigenfrequenz auf a} ~ 1394 Hz (—0,2%) herab. Wird diese Scheibe auBerdem 
am Innenrand elastisch mit einer Bettungsziffer ¢ = 3G gelagert, so ergibt (18) als noch tiefere 
_Grundeigenfrequenz aj ~ 1110 Hz (—1,6%)?. 
Als drittes Beispiel seien die drei ersten Torsionseigenfrequenzen der Scheibe mit dem ,,Sonder- 
profil nach Grammel (ebendort Seite 39, Abb. 34) berechnet. Hier ist der Dickenverlauf durch 
= 104 e"/!5 7-8 gegeben, wo y und rin cm zu rechnen sind. Ferner ist ry = 10 cm und r, = 50 cm. 
Ist die Stahlscheibe am Innenrand starr eingespannt und am Aufenrand frei, so ergibt Formel (14) 
folgende Werte: a, ~ 1090 Hz (—0,16%), a» = 5980 Hz (—0,12%) und a, ~ 10 112 Hz (—0,02%). 
-DaB diese letzten Naherungswerte so besonders genau werden, liegt daran, daf’ Formel (14) da- 
durch hergeleitet ist, daB das vorhandene Profil durch eines mit dem Dickenverlauf y = c, e%" r 
_approximiert wird, wo c, und ¢, geeignet gewahlte Konstanten sind. Fir ,sonderprofile® nach 
-Grammel wird die Formel (14) also besonders genau. 
SchlieBlich seien Naherungswerte fiir die Grundtorsionseigenfrequenzen von Scheiben mit dem 
r\P 
Profil yy == aie) berechnet, wo c = 2,517,, T/T, = 0,05 (0,05) 0,25 und p = | oder 2 ist. 
-Exakte Werte zum Vergleich finden sich in einer kiirzlich erschienenen Arbeit von Tameroglu®. 


| 1 Kin positives Zeichen bei dem Fehler bedeutet, daB der Naherungswert zu groB ist. 
2 Obgleich bei diesem Beispiel |Axp,/%p,| wesentlich groBer ist, als dies nach unserer Empfehlung sein 
sollte, ist die Genauigkeit von «* iiberraschend gut. 
3 §. Tameroglu, Ing.-Arch, 26 (1958), S. 212. 
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SchlieBt man die Profile mit r,/r, = 0,05 aus, so erfiillen alle iibrigen die Ungleichung (14a). Wen-, 
det man auf diese iibrigen die Formel (14) unter Beachtung von (14b) an, so erhalt man die Grund-; 
torsionseigenfrequenzen mit einem Fehler von absolut héchstens 2,1%. In den beiden wegen (14a) 
ausgeschlossenen Beispielen fiir r)/r, = 0,05 hatte man fiir p = 1 einen Fehler von —1,4% und fiir: 
p =2 einen Fehler von —10% erhalten. Die Beispiele von Tameroglu zeigen, daB Profile mit klei-; 
nem r,/r, am leichtesten auBerhalb des zulassigen Anwendungsbereiches der Formel (14) [mit oder 
ohne Beriicksichtigung von (14 b)] liegen kénnen. Fir r,/r, => 0,15 ist der Fehler bei den Grund-; 
torsionseigenfrequenzen fiir Scheiben mit den Profilen y = Y in den angefiihrten Beispielen stets 
kleiner als 1%, falls in (14) die Ersetzung (14.b) vorgenommen wird. 

Es sei noch angemerkt, daB die Dickenverteilungen y = y(r) in diesen Beispielen nur deshalb| 
in so einfacher analytischer Form vorgegeben wurden, damit es méglich war, die zum Vergleich. 
benétigten exakten Werte zu berechnen. | 


oe 


6. Schlu®bemerkung. Die in dieser Arbeit hergeleiteten Formeln sind Anwendungsbeispiele } 
fiir ein allgemeines Verfahren zur Aufstellung von Naherungsformeln fiir die Eigenwerte technischer } 
Eigenwertprobleme, das der Verfasser dargestellt hat!. Nach geeigneter Erweiterung kénnte man| 
mit dieser Methode bei Bedarf auch entsprechende Naherungsformeln fiir die Dehnungsschwin- | 
gungen einer Kreisscheibe aufstellen. In diesem Falle miiBte man die gegebene Differentialgleichung ; 
nicht wie in meiner friiheren Arbeit durch eine mit konstanten Koeffizienten approximieren, sondern | 
durch eine etwas kompliziertere Differentialgleichung, die jedoch in die Besselsche tibergefiihrt wer- | 
den kann. | 


(Eingegangen am 14, Juni 1958.) | 
Anschrift des Verfassers: Dr.-Ing. H. Wittmeyer, Linképing (Schweden), Valkebogatan 14 A. 


1 H. Wittmeyer, J. Soc. Industrial and Applied Math. 6 (1958) S. 111. 
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Biegeeigenfrequenzen eines gelenkig gelagerten Balkens 


Von H. Wittmeyer 


1. Einleitung. Die Biege- und Torsionsfrequenzen von geraden Balken bzw. Staben kénnen — 
wie bekannt— ohne Schwierigkeit mit den Verfahren von W. Ritz und R. Grammel oder durch 
Iteration berechnet werden. 

Fiir gewisse Randbedingungen kann man jedoch diese Frequenzen mit wesentlich geringerem 
Zeitaufwand berechnen, wenn man die einfachen Naherungsformeln benutzt, die der Verfasser in 
friheren Veréffentlichungen 12:4 hergeleitet hat. In Fortsetzung der durch diese Veréffent- 
lichungen begonnenen Untersuchungsreihe wird in der vorliegenden Arbeit eine einfache Formel 
fiir die angenaherte Berechnung der Biegeeigenfrequenzen fz,,., (n = 1, 2, 3, ...) eines an beiden 
Enden gelenkig gelagerten Balkens hergeleitet, dessen verdnderliche QuerschnittsgréBen zur Mitte 

zwischen den Auflagerpunkten symmetrisch sind. 

| Zur Herleitung der Naherungsformel wird im Prinzip die gleiche Methode angewandt wie in den 

oben erwahnten vier Veréffentlichungen. Diese Methode ist am allgemeinsten in der vierten dar- 
gestellt. In ihr approximiert man den ersten Eigenwert einer gegebenen Differentialgleichung 
durch den ersten Eigenwert einer Differentialgleichung mit konstanten Koeffizienten. Diese wird 

in geeigneter Weise aus einer ein- oder mehrparametrigen Schar solcher Differentialgleichungen 
ausgewahlt. 

Die in der vorliegenden Arbeit behandelten Randbedingungen weichen insofern von den damals 
-zugelassenen ab, als nunmehr an beiden Enden sowohl geometrische als auch dynamische (natiir- 
liche) Randbedingungen vorgeschrieben sind, wahrend friiher an einem Ende nur geometrische und 
am anderen nur dynamische Randbedingungen zugelassen waren. Bemerkenswert bei diesem An- 
-wendungsbeispiel ist, daB man in ihm die Differentialgleichung mit konstanten Koeffizienten aus 
"einer z weiparametrigen Schar auswahlen mu8, wahrend man bei den Beispielen der vier friiheren 
Arbeiten des Verfassers stets mit einer einparametrigen Schar auskam. 

) In Abschnitt 2 wird zunachst die in dieser Arbeit hergeleitete Formel in einer praktisch unmittel- 
bar anwendbaren Form dargestellt. Die Genauigkeit der Formel zeigen Beispiele in Abschnitt 3. 
Die Abschnitte 4 bis 6 enthalten die Herleitung der Formel. 


2. Die Naherungsformel. Bei dem beiderseitig gelenkig gelagerten Balken der Lange 2 / (Abb. 1) 
sei x der Abstand eines Querschnittes vom linken Auflager und € = x/I die zugehérige dimensions- 
lose Koordinate. Ferner sei 


| EI =K,(x) = Kgu Jj) (1) res ee 
‘die Biegesteifigkeit des Balkens, wo E der Elastizitatsmodul | | 


und IJ das Flachentragheitsmoment, bezogen auf die Biege- ty, A he 
at rox) ELM da 


-achse, ist und 

M(x) = My, m(6) (2) | 

die Masse des Balkens je Langeneinheit. Dabei besage der Index Abb. 1. Der beiderseitig gelenkig gelagerte, zur 
'M, daB die betreffende Gré8e in Balkenmitte (« = 1) zu nehmen ee Se 

ist. Im Intervall 0 < € < 1 seien j(€) und m(§) positiv und im 

Intervall 0 < & < 1 zweimal stetig differentiierbar. Ferner seij (1 + €&) =j (1 —) und m(1 + &) 
=m (1 —&) far 0 < £ <1. Die Biegesteifigkeit und die Masse je Langeneinheit seien also zur Bal- 
kenmitte symmetrisch. In Balkenmitte wird die Stetigkeit der Ableitungen von j(g) und m(é) 


nicht gefordert. 
Zur Frequenzermittlung berechne man zundchst mit 


pest es 
| /mei) pee ut PSO 3 
B, = ie)” &=7 (i= 0,1,.2; 3) (3) 


1 H, Wittmeyer, Ing.-Arch. 20 (1952), S. 331. 

2 H. Wittmeyer, Z. angew. Math. Mech. 36 (1956), S. 355. 

3 H, Wittmeyer, Forsch. Ing.-Wes. 24 (1958), S. 37. 

4 H. Wittmeyer, J. Soc. Industrial and Applied Math. 6 (1958), S. 111. 
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die Werte! . 
B = 75 (By + 4B, +2 By +4 By +1), (4) 
ean ae es fe 
Ep a i 9 (By ok B,) ’ | 
L, = {lg [m(Ex) 7(Es)]}? » (6)) 
j(9) : 
Ma 18 moo) ( 


Tabelle 1. Interpolationstabelle fiir die in (9) und (10) benétigten Funktionen 
F, = 0,62666 — 0,012556 Lp + 0,0001097 L2; F, = 1,25331 — 0,014522 Lp + 0,0000588 Lp . 


Lp F, te yAlind F, AF, 
0 | 0,6267 1,2533 
== i25 —145 
1 | 0,6142 | 1,2388 
| —122 —144 
2 | 0,6020 1,2244 
—120 —142 
3 0,5900 | 1,2102 
ill —140 
4 0,5782 | 1,1962 
| EAL lGual —140 
5 | 0,5666 822 
L113 = 139 
6) 018553 1,1683 
| WIE ——137 
7 0,9442 | 1,1546 
Im folgenden wird vorausgesetzt?: 
Le & 350 3 L105 (8}) 
Dann ermittle man mit FL) (t = 1,2) nach Tabelle 1 die GréBen 
Oo; = F(L,) — L, (0,073100 — 0,0011195 L,) + Lj (0,001806 + 0,0000610 L,) > (9})) 
Og = F(L,) + Lj (0,005275 — 0,0000181 L,) (10) 


und weiter fiir ungerade n > 1 


I 
Ce y v+L, (0,093758 — 0,119377 =| — L, (0,091617— 0,013185 L,), (11) 
sowie fiir gerade n > 2 
1 
One x ia one Bs (0,093758 — 0,238754 = + 0,013185 Lj. (12) 
Mit 


erhalt man schlieBlich als Naherungswerte f,, fiir simtliche Biegeeigenfrequenzen San, ex 48 Bal: 
kens die Werte i 


Le ees | 
fan = Boi Va (n= 2, SHS (14)} 


Wie sich spater ergeben wird, liefern diese Formeln die genauesten Werte fiir die Grundbiege-} 
eigenfrequenz (n = 1) und fiir die Kigenfrequenzen héherer Ordnungszahl n. Besonders gute Ge-+ 
nauigkeit ist zu erwarten, falls j(&) und m(é) monotone Funktionen sind, da die Massen- und Steifig-; 
keitsverlaufe, fiir die (14) exakt gilt (s. Abschnitt 7), ebenfalls monoton sind. 

Die Genauigkeit der Frequenzen fp, ist fir gerade n etwas unsicher. Die entsprechendent 
Formeln sind hier nur wegen des Fehlens besserer, gleicher Einfachheit angegeben. 


* Hier und im folgenden bedeutet lg den Logarithmus zur Basis 10 und In den Logarithmus zur Basis e 

* Der zweite Teil von (8) enthalt die Forderung, daB das Verhaltnis von Biegesteifigkeit zu Masse j¢4 
Langeneinheit in der Mitte nicht kleiner sein darf als an den Enden. Diese Forderung diirfte in der Praxi 
im allgemeinen erfiillt sein. 
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Ersetzt man bei einem Balken M durch K,~1 und Kz durch M~ so ersieht man aus (1) bis (14), 
dafs sich fz,, dadurch nicht andert. Dieselbe Eigenschaft hat auch der exakte Wert fp, x» Wie 
man nach dem Vorgang von Hohenemser! direkt aus der Differentialgleichung fiir den schwingenden 
Balken herleiten kann. 


3. Beispiele.2 In den Beispielen zur Anwendung der Formel (14) sind Naherungswerte fp, fiir 
die Biegeeigenfrequenzen fz,,,, eines beiderseitig gelenkig gelagerten Stabes der Lange 21 be- 
rechnet, dessen Querschnittswerte zur Mitte zwischen den beiden Auflagern symmetrisch sind. Die 
Frequenzbeiwerte §,, nach (13) sind mit den entsprechenden exakten Werten 


pees = fBn,ex i pe (15) 

verglichen und die relativen Fehler ig 
Pa ars iS ex 

.=- (16) 


n,ex 


Tabelle 2. Exakte Frequenzbeiwerte B,,,,,, und relative Fehler y,, der zugehérigen Ndaherungswerte f,, nach (13) 
fiir sieben beiderseitig gelenkig gelagerte, zur Mitte symmetrische Balken. 


Relative Fehler in Prozenten 


Dimensionslose | Dimensionslose Exakte Frequenzbeiwerte : 
r Sanaa Saas? ; - Para- SPSS nach (16) bei Anwendung der 
Seat ee snot nach @) meter ach (15) Naherungsformel (14) 
i® m(E) eS 
Bl, ex B2, ex | B3, ex BA, ex G1(%) | 2(%) G3 (%) | Me (%) 
1 ba : 2 |0,3872 1,554 3,55 | =20.05) 0321 20.30 
ESe ‘ : eae sad ae 
zm Gl 1+cé i Ge 
2 Se Eis ee a aes 4 | 0,3832) 1,539 | —0,11) 0,6 
pega] 1+ alee | mae 
3 ee 10 |0,3779| 1,516 | 3,54 | 6,26 |—0,13| 1,2) 1,4 1 
| 
' o _ 7 i : . a = a 
gags 
aS hse 
SB 4 —3(E—)4—3 E— 1) 
4.\823 3 | —___1+—_—_ — |0,3891| 1,555 | 3,52 | 0 0,5| 0,7 
= oo gd 4 4 | | | 
, N ky & 
dae 3 
Sie ee ae ; 
iz é (eA ee 
tk Py 
5 | && 5 2 |0,2828 0,954 | ==0,01] =1,0 | 
Begs | a+cs | 1+ce cra | 7 
6 |zeaz) a+ Cc) x 4 o2Hae 0,783 | 2-0: 0121-8 
fas] | 
sais eee ; ia oy | i i 
7 [Ag 10)" 1051997) 0,601L-121.70: | 2,83. |=, |—3,8| 0,3 |—4,3 


von f,, ermittelt. Die fiir verschiedene Balken errechneten Werte f,,,,, und ~, (diese in Prozenten 
-ausgedriickt) sind in Tabelle 2 aufgefiihrt. Die Genauigkeit der Naherungsformel erweist sich in 
diesen Beispielen fir den Grundton bei einem Fehler von absolut genommen héchstens 1,1% 
als recht befriedigend. Die berechneten Oberténe fz, niedriger Ordnungszahl n zeigen geringere 
Genauigkeit. Die bis zu n =4 berechneten Oberfrequenzen f,,, des Balkens Nr. 7 mit stark 
-veranderlichem Querschnitt weisen Fehler bis zu absolut genommen 4°% auf. Besonders schlechte 
Genauigkeit haben bei diesem Balken die Frequenzen mit gerader Ordnungszahl n. Daf dies kein 
‘Znufall ist, zeigt die Herleitung der diesbeziiglichen Formel in Abschnitt 6. Wie bereits erwahnt 
‘wurde und genauer aus den spateren Ableitungen in den Abschnitten 5 und 6 hervorgehen wird, 
‘kann man fiir sehr groBe n wieder eine bessere Genauigkeit der Naherungswerte fg, erwarten. 


| 2 Diese Beispiele sowie die hierfiir bendtigten Funktionen und Konstanten in den Abschnitten 4 b) und 
6 wurden von Herrn Ing. G. Andersson und Herrn B, Franzén berechnet. 
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Beispiele ohne weitere Rechenarbeit verdoppeln. Die so erhaltenen neuen Beispiele betreffen dann | 
Balken, die in ihrer Mitte eine kleinere Masse je Langeneinheit und eine kleinere Biegesteifigkeit | 
haben als an den Enden. Das Verhaltnis von Biegesteifigkeit zu Masse je Langeneinheit ist jedocl | 
auch bei den neuen Beispielen in der Mitte nicht kleiner als an den Enden. 

Die in Tabelle 2 aufgefiihrten ,,exakten Werte“ Ba -, wurden durch hinreichend oft wiederholea| 
Iteration mit anschlieBendem Ritzschen Verfahren berechnet, wobei alle vorkommenden Integra- 
tionen mittels der Trapezregel durchgefiihrt wurden. Um den HinfluB der hierbei benutzten Schritt- 
lange Aé von & auf die Resultate so weit wie méglich zu eliminieren, wurden die Frequenzen fiir | 
die drei Schrittlangen A& = 0,1; 0,05 und 0, 025 berechnet, als quadratische Funktionen von 
(AE)? dargestellt und nach (4)? = 0 extrapoliert. Der durch diese Extrapolation erhaltene Wert. 
wurde als der ,,exakte Wert‘ angesehen. Um zu zeigen, wie der durch numerische Integration er- 


haltene Naherungswert f, fiir Px, exakr VOn der Schrittlange Ag abhangt, ist in Tabelle 3 der Wert 


Os als Funktion von Aé fiir den Balken Nr. 1 aus Tabelle 2 aufgefiihrt. In der letzten Spalte von 
Tabelle 3 ist der durch Extrapolation gefundene ,,exakte Wert‘ angegeben. Man sieht, daB durch. 
diese Extrapolation bereits die dritte geltende Ziffer geandert wird. 


Tabelle 3. Abhdngigkeit des durch Iteration mit anschlieBendem Ritzschen Verfahren gewonnenen 


Ndaherungswertes 8, fiir den Frequenzbeiwert Px, ex von der Schrittlénge A§ bei den numerischen Inte- 
grationen, (Balken Nr. 1 aus Tabelle 2.) 


Schrittlange A& 0,1 | 0,05 | 0,025 | 0 (extrapoliert) | »xakter Wert“ 
———| | | | | 
Naherungswert /, | 0,38871 | 0,38760 | 0,38731 | mOss one Bi ex = 90,3872 


4, Herleitung der Naherungsformel (14) fiir die Grundbiegeeigenfrequenz. a) Formulierung 
der Aufgabe als VariationsproblemI. Die zux =/ symmetrische Durchbiegungsamplitude 
w = w(x) ist wegen der vorausgesetzten Symmetrie von K,(x) und M(x) Lésung des bekannten | 


Variations-Eigenwertproblems 
| 


5 | 5 [Ke (w")? 4 Mw] dx = 0, an) 
m 
wobei Striche Ableitungen nach x bedeuten, mit den geometrischen Randbedingungen | 
mi) (a0), (18) 
w == 0 C— (19} 


Die hier vorgenommene Beschrankung auf symmetrische Schwingungsformen werden wir zwecks 
Vereinfachung im ganzen Abschnitt 4 beibehalten und auch noch in Abschnitt 5 voraussetzen. 
Die ieeinere 4 =, hangen mit den exakten Kigenfrequenzen des Balkens der Lange 2 1 wie } 


folgt zusammen: 
= (20 fB ok —1,ex)" (k = 1,2, 3,...). eay| 


Durch die angefiihrte cies haben wir also zunachst die Kigenfrequenzen f,,, ,, mit ge-- 
rader Ordnungszahl n ausgeschlossen. Durch die Transformation 


w=, w= gq, w=, (21) 
fiihren wir statt w(x) und seiner Ableitungen die drei Funktionen q,, q, und q; ein. Sie hangen durck | 
die Beziehungen | 

qdi—gm=9, q2 — 43 = 0 (22)) 
zusammen. Um die Funktionen q,, q, und q, in dem transformierten Variationsproblem unabhangig / 
voneinander variieren zu kénnen, addieren wir unter dem Integralzeichen von (17) die linken Seiten 1 


der Gleichungen (22), noch multipliziert mit Lagrangeschen Faktoren p, bzw. p,. Damit ergibt die} 
jean ton (21) schlieBlich 


1 1 


; i 9 6 ! t 
3 | 3g (Kp a3 —2 M gi + 2 py (9, — 4%) + 2 Po (4s —95)] dx = 0 (23) 


| 
0 


mit den Randbedingungen 


n=0 (x=0), w=0 (=). (24)) 
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In einer kleinen Zwischenbetrachtung sei die physikalische Bedeutung der hier eingefiihrten 
GréBen p, und p, geklart. Dazu schreiben wir die Eulerschen Gleichungen von (26) auf, die durch 
Variation von q, und q; entstehen: 


Pi +p =0, (25) 
Kp qs — pp = 0. ; (26) 

Hieraus folgt bei Beachtung von (21) 
Pi =—(Kgw’), P2 = Kgw". (27) 


Bei geeigneter Vorzeichenfestsetzung stellt also die Gréfe p, die Querkraft und die GréBe Pz das 
Biegungsmoment dar. In der folgenden Herleitung arbeiten wir also gleichzeitig mit den Defor- 
mationen q,, g, und g; sowie den zu den beiden ersten gehérenden ,,generalisierten Kraften“ p, 
und p,. Der Gedanke, diese Gréfen bei der Berechnung von Biegeeigenfrequenzen gleichzeitig zu 
benutzen, findet sich bereits bei Myklestad! und geht auf einen Gedanken von Holzer? zuriick 3, 
der vor Mykelstad ein entsprechendes Verfahren bei Torsionsschwingungen anwandte. In neuerer 
Zeit wird dieser Gedanke in steigendem MaBe benutzt*>®, 

Nun transformieren wir das Variationsproblem (23), (24) in eine Gestalt, bei der nur noch 
»Weniger wichtige Koeffizienten von der unabhangigen Variablen abhangen. ,,Weniger wichtig“ 
sind dabei solche Koeffizienten, deren Abanderung das asymptotische Verhalten der Eigenwerte 
in erster Naherung nicht andert. Eine solche Transformation, die p,, Py, G1, qo, q3. x und A in P,, 
P,, Qi, Qo, Qs, 7 baw. wu tiberfiihrt, 14Bt sich sofort als Spezialfall einer allgemeineren angeben, die 

in Kapitel 4 der vierten zitierten Arbeit des Verfassers angegeben ist: 
” 


3a) (28) 


T 
0 
eae th (29) 
p= 0 Pi 6G = 1,2), (30) 
; qi = C; Q; (==. 3) (31) 
Hierin sind mit der Bezeichnung 
i 1 
M sory 
folgende Abkiirzungen benutzt: 
GG = Ves ? (33) 
M\z 
| T iy 4 
ose(a 0 
| Sen el M\ei;—t 
| 6. = JIE BF («| 2 (= 1, 2.3). (35) 
; — 
Fiihren wir noch die Bezeichnungen ein 
1 1 M\ 
h, = {ln (¢, Ke M), hg = ln (te), (36) 


worin e, und e, positive Konstanten sind, die nur so gewahlt werden miissen, da die obigen Loga- 
rithmen dimensionslose Argumente haben, so ergibt die Transformation (28) bis (31) das 
Variationsproblem I: 


. 
; ° it x : 
6[ [PQ + PQ + + @— P, Q,— Pa Qa — 5 HO —h, (P, O + Ps Q) 
0 


ier, OP 0) dy =0, (37) 


1 N.O. Myklestad, J. Aeronautical Sciences 11 (1944), S. 153. 

2 H, Holzer, Die Berechnung der Drehschwingungen, Berlin 1921. 
3 Herr Professor K. Liirenbaum machte mich freundlicherweise darauf aufmerksam, dal} Holzer diesen 
Gedanken wiederum von Giimbel (Z. VDI 1912, S. 1025 u. 1085) tbernommen hat. 
4 FE, Pestel, Abhandl. Braunschweigischen Wissenschaftl. Gesellschaft 6 (1954), S. 227. 
5 K. Marguerre, J. Math. and Phys., 25 (1956), S. 28 (hier weitere Literaturhinweise). 
8 §, Falk, Ing.-Arch. 26 (1958), S. 96 (hier weitere Literaturhinweise). 
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wobei Punkte Ableitungen nach 7 bedeuten, mit den Randbedingungen 


Q=9 (7 =9), (38) 
Q=9 (=1). (39) 
Damit die Lésungen eindeutig werden, fiigen wir noch die Normierungsbedingung 
1 
{ @dy=1 (40) 
6 


hinzu. Die nach aufsteigender GréBe geordneten Eigenwerte seien mit “4 = um, (k = 1, 2, 3,...) 
bezeichnet. 

b) Das dem Variationsproblem I ,,benachbarte“ Problem I. Der erste Kigenwert yu, 
des Variationsproblems I soll durch den ersten Eigenwert jZ, eines einfacheren Variationsproblems II 
angenahert werden. Problem II entsteht aus Problem I, wenn man hierin hy durch die Kon- 
stante b, und h, durch die Konstante 6, ersetzt. Die Eigenwerte des Problems II seien mit 4, 
(k =1,2,3,...) bezeichnet. Die zu dem kleinsten Eigenwert 4, = [, (b,, b,) gehérigen Eigen- 
funktionen seien P,, P;, Q,, Q, und Qs. 

Um die giinstigsten Werte von b, und b, zu ermitteln, bilden wir zunachst ein Zwischenpro- 
blem III, das zwei Parameter ¢, und €, so linear enthalt, dafB es fiir ¢, = ¢, = 0 in Problem II und 
fiir €, = €g = | in Problem I iibergeht: | 


| 


1 
s[[Ra+?, 6, +5 8— P, O— Pr — >» OF — Ths + es (ey — BN (P,Q + Pa O,) 


— [by + &5 (tg — by)] (P, Q, — Ps Q,) dy = 0. (41) 


Als Randbedingungen gelten (38), (39). Die Normierungsbedingung (40) wird ebenfalls gefordert| 


In erster Naherung wird 4, ~ /, sein, falls man b, und b, so bestimmt, daB | 
| 


Gh, (0, 03 by, be) : 
EE oh (i = 1,2) (42) 


ist. 
Um mit Hilfe von (42) die giinstigsten Werte bj, 6 von 6,, b, zu bestimmen, beachten wir, dail 
fiir den stationaren Wert A, des unter dem Variationszeichen 6 stehenden Integrals in (41) 

aA, 


=? (i = 1, 2) (43) 


gilt, was in ahnlicher Weise wie in Kapitel 2.3 der vierten oder in Kapitel 4.4 der dritten Arbeit 
des Verfassers gezeigt werden kann. Dabei sind bei der mit 0 bezeichneten partiellen Differentiatiori 
die Kigenfunktionen P,,..., Q; als konstant anzusehen. Aus (43) ergeben sich fiir ¢, =e, = 0 untes 
Beachtung von (42) die folgenden Bestimmungsgleichungen fiir die ,,giinstigsten Werte“ i b' 
und b, = 6} im Variationsproblem II: | 


bs a eee : (455 
J (PQ. — P,Q) dn 


Weil P,, Py, Q; und Q, in (44) und (45) die Eigenfunktionen des Variationsproblems IT sind|I 
hangen auch die rechten Seiten von (44) und (45) noch von b, und b, ab. Bevor wir b* und b% aud 
diesen beiden Gleichungen bestimmen kénnen, miissen wir die Lésung des Variationsproblems I/ 
untersuchen. Das soll im nachsten Abschnitt geschehen. | 


c) Untersuchung des Variationsproblems II. Die zu #, = 2, (b,, by) gehérenden Eigen 
funktionen P,, P,, Q;, Qz und Q; des Variationsproblems II, das sich aus (38) bis (41) fiir e, =e, =4 


| 


|! 

/ 
} 

| 
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ergibt, geniigen den folgenden Eulerschen Gleichungen: 


G——(, +) Q=0, 
G0, 05) 0, 0) 
P+h+h +h) B=0, | (46) 
P, + P, + (b, —b,) P, = 0, 
Q, — P, =0. 


and den Randbedingungen 


@Z=P=0 (=), 0, = Pee 0 Gj =1). (47) 


Be intern wir aus (46) und (47) die GréBen Pi pe Qo und Q, so erhalten wir die Differential- 
chung 


Q, —2 (02 + 82) 0, —[z, — (8? — 83] 0, =0 (48) 
mit den Randbedingungen 


me 
2 @ hen = 0); 
Q,—26,0,=0, 
- a (49) 
Felgen. | 
Q, — (6; + bs) Q%=0, 
Der iibliche Exponentialansatz fir Q, in (46) fihrt zu der Darstellung 
: Q, =asinan + beosan +e c Gin By + dCoj by (50) 


nit zunachst noch unbestimmten Konstanten a, b, c und d sowie 

o = [(u, + 4 BF 8)" — (bf + BS)" und B = [(, + 4.87 5)"? + BF + BR. (51) 
is Bestimmung von a, b,c und d erhalten wir durch Kinsetzen von (50) in (49) das lineare Glei- 
shungssystem 


a 0 1 0 1 \ («) 


—2b,a — 7 —26,6 Ge b 
[x cos a [—« sina [6B Co} [6 Sin B ¢ |=0. (52) 
— (b; + by) sina] — (b, + by) cos a] — (b, + b2) Gin B] — (b, + be) Coj B] 
[— a? cos a [a3 sin « [6? Coj 6 [6* Gin p | q 
ut (6 + bp) a? sina] + (b; +b.) a? cosa] — (b, + bs) 6? Gin B] —(b, + b,)6?Coj 6] 


Die Forderung, daB die Determinante dieses homogenen Gleichungssystems fiir die Unbekannten 
1, b,c und d verschwinden muB, fiihrt nach einigen Vereinfachungen zu der Gleichung 


b, B (4b? — ao? — ?) sin « Co) B + (63 — 62) (a? + B) sina Gin B + a B (a? + f?) cos a 

Cof B — bw (4 b? + a? + f?) cosa Gin B = 0. (53) 

Aus (53) und (51) lassen sich fiir ein vorgegebenes Wertepaar 6,, b, die zugehorigen Werte (4, 

x, 8 und damit unter Beachtung von (52) auch das Verhiltnis a:b:c: d bestimmen. Aus (50) und 

40), worin Q, durch Q, zu ersetzen ist, kann dann Q, berechnet werden. Die iibrigen Eigenfunktionen 

P,, P,Q, und Q; ergeben sich mit Hilfe von (46). Den spater benétigten Wert 117'*/(2 7) )”? kann 

man in dem praktisch ausreichenden Bereich |b,| <2, —1 Sb, 0 mit 1% Genauigkeit’ wie 
‘olgt darstellen: 


a, = “4"* — 0,62666 — 0,021315 i? + 0,0003161 bt + b, (0,006603 b2 — 0,253976) 


= (2 ge) !8 
+ b2 (0,001249 bf + 0,021801) . (54) 
Die aus den Integralen ; 
= { (P,Q, + P,Q) dy, ee | (55) 
1 —_—_ — —_ — — a —_ — 
K, = J (PQ + Pe) dn, K, = {n(B.G— P,Q) dy 


1 Die folgende Polynomapproximation (54) hat ebenso wie (84) fiir kleine Werte von |b,|, |b,| bzw. |b;|, 
b | wesentlich groBere Genauigkeit als angegeben ist. 
Qg* 
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gebildeten spater benétigten Quotienten K|/J, und K,/J, erweisen sich im Gebiete |b,| < 2 
—1<b,<0 als nahezu konstant. Zwecks Vereinfachung der spater herzuleitenden Formel| 
werden daher diese GréBen als konstant betrachtet. Wir setzen 

C= = = 0,625, Coe 


a 2 


== (ei (56 

d) Approximation des ersten Eigenwertes des VariationsproblemsI durch de) 
des Problems II. Um nun das Gleichungssystem (44), (45) angenahert nach b, und b, aufzv 
lésen, setzen wir fir hy und he folgende abgebrochene Taylorreihen mit noch unbestimmten GréBen «4 
und ¢, an: 


h(n) = h(e,) + hoa) y—a), 9 SqS1, | 57 
hg() = hy(c2) + hg(ea) (77 — 2) » {=o—f 
Durch Einsetzen von (57) in (44) und (45) erhalten wir unter Beachtung von (55) und (56) 
= hy(c,) a h,(c,) (Cj —4), (54 
by = ho(co) a hig(¢2) (C, — ¢) . (54 


Wir setzen nun c, = C, und c, = C, gemaB (56) und ersetzen auBerdem die ersten Differentia: 
quotienten durch Differenzenquotienten gemaB der Formel 
d 1 f 

Sj =e (fe Pate (64 

Darin beginnt das fortgelassene Glied erst mit der dritten Ableitung, die auch in (57) nicht be 
riicksichtigt ist. Wahlen wir in h,(y) als Schrittlange im Sinne von (60) die GréBe Ay, = 1 — Gi 
und in h,(7) die GréBe An, = 1 — C,, so erhalten wir, wenn wir noch (56) beachten, die N aherungy 
formeln / 


« — MQ)—h(2G—1) _ byl) — h,(0,25) | 
PLES Gea g Tee Targa 
i) he d 
by e 2 (1 ar ) = h(1) — h,(0) . (64 


Berechnen wir zu diesen Werten b, = bj und b, = 6; mit Hilfe von (54) den ersten EHigenwes 
fy = Pi des zugehérigen Variationsproblems II, so gilt nach Abschnitt b) | 

[iy [i « (64 
Fiir die rechnerische Auswertung von (61) ist es noch giinstig zu wissen, welches £ = (x/I)-Argi 
ment € = , zu dem 7-Argument 7 = 0,25 gehért. Da wegen (28) die Beziehung d = dx/l = di) 


gilt, so wird wegen (34) 
é 
r M \1/4 
n=y | (x) dé. (64 
0 


Bezeichnen wir den zu € = 0,25 gehérigen 7-Wert mit 7,, so wird 
2 


525 
I M \114 
a | (x, dé. (65 
6 


Werten wir hierin das Integral nach der Trapezformel aus und J nach der Simpsonschen Regt 
so wird mit den Bezeichnungen (1), (2), (3), (4) 


eee My, \"!4 | 
qeeres x) (66 
und damit weiter 
1/B,+ B 
m= 5(- ; ‘ é 0,25 ! (6% 
Interpolieren wir 7 als Funktion von é linear in dem Intervall 0 < é < 0,25, so erhalten wir 
LBB | 
po (Be a 


Hieraus errechnet sich der zu 7 = 0,25 gehérige Wert € = £,, wie er in (5) angegeben ist. Die | 
unter der Voraussetzung 0 < &_ < 0,25 hergeleitete Beziehung wird auch noch gelten, falls i 
etwas gréBer als 0,25 wird. 


} 
] 


| 
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j Nun nehmen wir noch einige Umformungen vor, um auf die in Abschnitt 2 angegebene Form 
fix die Berechnung der Grundbiegeeigenfrequenz zu kommen, Dazu setzen wir in (36) 


ey = (Kem My)", ey = My Kum (69) 

und fiihren die Bezeichnungen 
(3 by Ig e)? = De; (70) 
8brlge =L, (71) 
ein. Damit geht (61) nach Quadrieren in (6) iiber und (62) in (7), wenn wir (36), (69), (1) und (2) 
beachten und aufferdem uns erinnern, daB ny = 0,25 dem €-Wert € = €, nach (5) entspricht. Be- 


zeichnen wir die zu b, = bj, b, = bi gehérigen Eigenwerte des Variationsproblems II mit 7 und 
setzen zur Abkiirzung 


(uz 1/4 
Con =o (k = 1,2,3,...), (72) 


so erhalten wir fiir k = 1 bei Beachtung von (54), (70) und (71) 
a = 0,62666 — 0,012556 L, + 0,0001097 L; + L, (— 0,073100 + 0,0011195 L,) 


+ Lj (0,001806 + 0,0000610 L,) (+ 1%) . (73) 
Tabulieren wir noch den von L, unabhangigen Teil 
F\(L,) = 0,62666 — 0,012556 L,, + 0,0001097 Lz, (74) 


so kommen wir zu der in (9) angegebenen Darstellung von a,. Aus (29), worin A durch A, und u 
durch ju, zu ersetzen sind, sowie (20) erhalten wir 


| (2 70 fis, 2k 1,02) = An = % Me (k = 1, 2,3,...). (75) 
Daraus folgt fiir k = 1 wegen (63) und (72) 

— be — 
) Spies Ly x acs = x Cre (76) 
Bezeichnen wir den erhaltenen Naherungswert fiir fp,,,, mit fg,, so erhalten wir schlieBlich unter 
Beriicksichtigung von (33), (66) und (13) die Gleichung (14) fiir n = 1, womit die Giiltigkeit dieser 
Formel fiir die Berechnung der Grundbiegeeigenfrequenz gezeigt ist. 


5. Herleitung der Naherungsformel (14) fiir die Oberbiegeeigenfrequenz fj , ungerader Ordnungs- 
zahl n. Die aus den Eulerschen Gleichungen des Variationsproblems I nach (37) bis (40) her- 


a Differentialgleichung fiir Q, hat — wie man leicht zeigt — dieselben Faktoren bei 0, und 


4Q, wie die entsprechende Gleichung (48) fiir das Variationsproblem IT bei Q, und “4, Q,. Nach 
Lehmann! verhalten sich daher die zu den Problemen I und II gehérigen Eigenwerte yu, baw. [, 
i k — o in erster Naherung asymptotisch gleich, und zwar unabhangig davon, welche Werte 
die GréBen b, und 6, in dem Variationsproblem II haben. 

Setzen wir nun in Variationsproblem II die GréBen b, und b, gleich bj baw. by und bezeichnen 
die zugehérigen Werte mit yi (k = 1, 2,3,...), so liegt es wegen des genannten asymptotischen 
Verhaltens nahe, die Beziehung (63) in folgender Weise zu erweitern: 
by Wie (ee 12, Bs en) < (77) 
Diese Beziehung wird dann sowohl fiir k = 1 als auch fiir k + o0 recht gute Naherungswerte 
fiir die Eigenwerte jz; des Variationsproblems I liefern. Wegen der Einfachheit von (77) wollen 
wir diese Beziehung zur naherungsweisen Berechnung von ju, fiir alle k benutzen und dabei in 
Kauf nehmen, daB diese Beziehung fiir kleinere Werte k > 1 méglicherweise nur grobe Naherungs- 
werte fiir ju, liefert. Zur Berechnung von jij fiir k > 1 soll die folgende, aus (51) und (53) herleit- 
bare asymptotische Beziehung (mit 2 k — 1 = n) 


ee at ( 4 ini iea is 0(53] (78) 


296 nN It Dye 


benutzt werden, in der dann b, = bj und b, = b} zu setzen sind. 

Ebenso wie aus (63) die Beziehung (14) fiir n = 1 hergeleitet wurde, laBt sich aus (77) die 
Beziehung (14) fir n = 2k—1 (k = 2,3,4,...) herleiten. Dabei geht (78) mit Hilfe von (70) 
und (71) in (11) itber. Damit ist die Formel (14) fiir alle ungeraden n = | hergeleitet. 


1N, J. Lehmann, Wissenschaftl. Z. Techn. Hochsch. Dresden 2 (1952/53), Heft 3. 
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6. Herleitung der Formel (14) fiir die Oberbiegeeigenfrequenzen fp; » gerader Ordnungszabl n. D 
Eigenschwingungsformen, die zu den Kigenfrequenzen fz,,,,, mit geradem n gehéren, sind Z 
Balkenmitte antimetrisch. Daher sind sowohl das Biegungsmoment als auch die Durchbiegury 
in Balkenmitte, d. h. fir x =I, gleich Null, und der bei x =/ durchgeschnittene und bei x = Hi 


und x = 1 gelenkig gelagerte Balken hat Eigenfrequenzen ips fiir die 
pee = fp orex (r — 1, 2, 3, ft .) (77 


gilt. Es bleibt also die Aufgabe, Naherungswerte hee fiir die Eigenfrequenzen fz, ,, eines bei « =I 
und x =I gelenkig gelagerten Balkens zu bestimmen, der in 0 [x S/ dieselben Querschnitt 
gréBen hat, wie der von uns bisher betrachtete Balken. Die die antimetrische Schwingungsforr 
darstellende Funktion w = w(x) geniigt fir 0 < «<1 natiirlich auch der Variationsbedingum 
(17) sowie der Randbedingung (18). Die Bedingung (19) ist jedoch durch / 

0 (Cpe) (86 


zu ersetzen. Durch dieselben Transformationen wie in Abschnitt 4a) kann dies Variationsprobler 
auf die Form (37) gebracht werden. Als Randbedingung ist auf er (38) die Bedingung 
%=0 j@=v ) (83 
hinzuzufiigen. | 
Um nun zu einfachen Naherungsformeln fiir die Eigenwerte dieses Problems zu kommer 
setzen wir in diesem Variationsproblem | 
Reet ee Vivi (8 

wo bf und b3 die Konstanten sind, die sich als Lésung von (44), (45) ergaben. Die Eigenwerte | 


dieses vereinfachten Problems ergeben sich aus den folgenden Gleichungen, die in ahnlicher Weis 
hergeleitet sind, wie die Gleichungen (51) und (53): 


8 bi? a B — [4 bi? (0? — f%) + (a? + B2)] sina Gin B — 8 bia B cosa Coj f =0, 
fp? — 0? = 2 (by? SF b>?) ? 


a == (2 + B2)2 — 4 bf? B32. 


‘ 


Die Lésungen dieser Gleichungen kénnen im Intervall 


ers2, —1<b9<0 
wie folgt approximieit werden: 
47. \1/4 
i oon — 1,25331 — 0,024651 b%2 + 0,0001695 bx4 
at a rr + (0,063674 — 0,000370 b#2) b32 (+. 0,15%) a 
undvmit tree (reo ee 
ie (ey nt b*2 8 p*2 1 
tte Hig hel aah alae Gy 


Die Eigenfrequenzen Pe des durch (82) vereinfachten Systems berechnen sich damit, wie mag 
abnlich wie in den Abschnitten 4d) und 5 schlieBt, zu 


< a2 ay Rap we : 
for = Fie Vast (3 ar 12,94.) (84 


tipetes hate (87 
so gilt (87) exakt, wenn (82) exakt erfiillt ist. Allgemein gilt (87) wenigstens in erster Naheruni 
symptotisch. Um zu einer einfachen Formel zu kommen, wollen wii (87) als unsere Naherungs 
forme] ansehen. Fiihren wir in (84) und (85) noch die Transformation (70) und (71) ein, so i 


halten wir (10) bzw. (12). Damit haben wir, wenn wir noch (79), (86) und (87) zusammenfassen 
die Gleichung (14) auch fiir gerade n > 0 hergeleitet. | 


q 


. ae Balken, fiir die die Niherungsformel (14) exakte Werte ergibt. Zur Herleitung der Naheranes 
formel (14) dieser Arbeit wurde die Gleichung (37) durchweg durch die Annahme (82) mie 


Setzt man nun 


Bis auf die geringen Fehler, die durch die Polynomapproximationen der a, (n = 1, 2, 3,.. .) mitte 


der Gleichungen (9) bis (12) entstehen, liefert (14) also exakte Werte, wenn (82) fiir den Balk 
exakt erfiillt ist. 
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Die Massen- und die Steifigkeitsverteilung eines derartigen Balkens erhalt man durch Inte- 
gration von (82), wenn man noch (36), (28) und (34) beachtet. La®t man noch den Stern bei den 
nunmehr heliebigen Konstanten b, und 6, fort, so ergeben sich die QuerschnittsgréBen M und 
Kz der Balken, deren Eigenfrequenzen durch (14) (nahezu) exakt geliefert werden, wie folgt: 

a) Falls b, + 0 ist, 


by 
by 


Ky = Kgy e*—?%[1 —(1—e"2®) ]  *, M= My e—**—?% [1 — (1 — e— 7%) E] ms (88) 


und 


b) falls b, = 0 ist, 
Ky = Kaye 720-9). M = My, e—740—*) , (89) 
Es sei darauf hingewiesen, da alle diese Massen- und Steifigkeitsverteilungen einen monotonen 
Verlauf haben. — Die Berechnung der Eigenfrequenzen fiir die Balken mit den Querschnittswerten 
nach (88) oder (89) mit Hilfe von (14) kann noch dadurch vereinfacht werden, daB man die be- 
notigten a,-Werte nicht aus (9) bis (12) berechnet, sondern aus den Gleichungen (54) und (78) 
sowie aus (84) und (85), wo in den beiden letzten noch bf und b¥ durch b, bzw. b, zu ersetzen sind. 


8. Bemerkung iiber Balken mit unsymmetrischem Querschnittsverlauf. Bei Balken mit un- 


symmetrischem Querschnittsverlauf kann man zunachst ganz ahnlich vorgehen wie in Abschnitt 4 
) bei dem symmetrischen Balken. Dies wurde bereits in Abschnitt 6 angedeutet. Bei der Wahl der 
Parameter b, und 6, wurde jedoch eine Schwierigkeit bewuBt umgangen: Diese Parameter wurden 


in Abschnitt 6 einfach gleich den fiir den symmetrischen Fall berechneten Werten bj bzw. b% 
gesetzi. Das waren jedoch keineswegs optimale Werte fiir den unsymmetrischen Fall in dem 
Sinne, wie es bj und 63 fiir den symmetrischen Fall waren. Eine in diesem Sinne optimale Wahl 
ist im unsymmetrischen Falle leider nicht méglich, da die GréBen, die C, und C, nach (56) ent- 
sprechen wiirden, fiir den unsymmetrischen Balken auch nicht annahernd konstant sind. Fir 
den unsymmetrischen Balken ist offenbar eine Vereinfachung des Variationsproblems in der Art 
(82) zu einschneidend; eine Approximation der gegebenen Differentialgleichung durch eine solche 
mit nicht oder nur teilweise konstanten Koeffizienten ware wiinschenswert gewesen. Eine solche 
Approximation wurde jedoch wegen des zu erwartenden Arbeitsaufwandes nicht in Angriff ge- 
nommen. 

Bei dem symmetrischen Balken konnten wir diese Schwierigkeit dadurch umgehen, dal} wir 
nur eine Halfte betrachteten. Dadurch bekam Q, an dem einen Ende des betrachteten Balken- 


_teils seinen Kleinstwert und an dem anderen Ende seinen Gré®twert. Diese Eigenschaft von Q, 


mui vermutlich zumindest angendhert vorhanden sein, damit unser hier benutztes allgemeines 
Verfahren zur Herleitung von Naherungsformeln mit Erfolg angewandt werden kann. Das zeigt 
jedenfalls die Erfahrung sowohl in der vorliegenden Arbeit als auch bei den Herleitungen von Na- 


_herungsformeln in den vier friiheren Arbeiten des Verfassers. Der Grund ist kurz gesagt folgender: 


UmfaBt die durch Q, dargestellte Schwingungsform mehr als ,,eine Viertelperiode einer sinusformigen 


-Schwingung“, so treten in den ,,Gewichtsfunktionen“ bei h, und hy in (44) baw. (45) Vorzeichen- 


wechsel auf, die die Auflésung dieser Gleichungen nach b, und 6, erschweren oder gar unméglich 
machen. 


9. SchluB. Die in der vorliegenden Arbeit gegebene Herleitung einer Naherungsformel fiir die 
Biegeeigenfrequenzen eines an beiden Enden gelenkig gelagerten Balkens, dessen Querschnitts- 
gréBen zur Mitte zwischen den Auflagern symmetrisch sind, zeigte in einem weiteren Beispiele 
die Brauchbarkeit eines friiher vom Verfasser entwickelten allgemeinen Verfahrens' zur Auf- 
stellung von Naherungsformeln fiir die Kigenwerte selbstadjungierter Differentialgleichungen. 
Eine Fehlerabschatzung, wie etwa fiir die Naherungsformel fiir Torsionseigenfrequenzen”, wurde 
fiir die hergeleitete Formel nicht gegeben, Ihre Brauchbarkeit konnte jedoch an Beispielen gezeigt 
werden. 

Bemerkenswert ist, daB die erwahnte allgemeine Methode nicht ohne weiteres auf den nicht 
zur Mitte symmetrischen Balken angewandt werden kann. Aus diesem Grunde konnten in dieser 
Arbeit fiir die Biegeeigenfrequenzen f,,,,. gerader Ordnung n nicht so genaue Naherungsformeln 
angegeben werden wie fiir die Frequenzen ungerader Ordnung. 


(Eingegangen am 14. Juni 1958.) 
Anschrift des Verfassers: Dr.-Ing. H. Wittmeyer, Linképing (Schweden), Valkebogatan 14 A 


1 Siehe FuBnote 4 von Seite 117. 
2 Siehe FuBnote 1 von Seite 117. 
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Von W. Engl 


1. Einfiihrung in die Behandlung des ebenen Spannungszustandes mittels hyperkomplexer Funk--| 
tionen. Wir stellen zunachst kurz die Ergebnisse der fiir den zu behandelnden Gegenstand grund- + 
legenden Arbeit von L. Sobrero! zusammen. Wenn man Lésungen f(x, y) der Gleichung | 

ee hn fas 
Af = 0x2 1 oy? : 
sucht und den Ansatz f(x, y) =f (x + iy) macht, so findet man, daB die Konstante i der Glei-- 
chung 1 + i? = 0 geniigen muB. Man stellt fest, daB fiir analytische Funktionen gilt f(x + iy) => 

= u(x, y) +iv(x, y) mit du = Av = 0. 

Sucht man die Bipotentialgleichung | 

Le aif. of Caf eS, l 
AAS = oe + 2 asaya + aye =O ()) 
ebenfalls mit einer Funktion f(x, y) der linearen Kombination z = x +] y zu befriedigen: 
floy) =fle tins )) 
so findet man, daB die Konstante 7 der Gleichung | 
L+2j+jt=0 ()) 
geniigen muB. Fiir analytische Funktionen erhalt man 


ee ede ries hae ede ) 
AAa = AAb = AAc = AAd = 0. (5)) 


Die den Cauchy-Riemannschen Differentialgleichungen entsprechenden Regularitatsbedingungen \) 
lauten: 


mit 


a, = — d,,, by = a, , c, =b,—24,, d, =¢,. (6)) 


Weiter 1aBt sich der Cauchysche Integralsatz beweisen: Ist f(z) in einem einfach zusammen- 
hangenden Gebiet iiberall regular, so gilt: 


f f(z) dz =0. (0) 
c 


i 
Aus diesem Fundamentalsatz folgt wiein der Funktionentheorie die Cauchysche Integralforme! | 
und die Entwickelbarkeit einer in einem Kreisringgebiet definierten Funktion in eine Laurentsche } 

y ae Reihe. Vergleicht man die Gleichgewichtsbedingungen des 
ebenen Spannungszustandes bei Vernachlassigung der | 
Massenkrafte 


00x OF OT O0y q 


mit den Cauchyschen Bedingungen (6) fiir die Funktion 
f@) =o + bape Fad, 
so sieht man, dal die beiden Gleichungen, welche b nicht. 
enthalten, mit den Gleichgewichtsbedingungen identisch | 
Abb. 1. Gleichgewicht an einem Element des Randes. sind, wenn man nur o, fiir —e, o, fir a und ¢ fiir d 
schreibt. 

Die Elemente a,c und d erfiillen also die Gleichgewichtsbedingungen und sind auBerdem Bi- 
potentialfunktionen. Diese beiden Bedingungen sind gerade notwendig und hinreichend, um die | 
Elemente a,c und d als Spannungskomponenten zu charakterisieren. 

Besonders wichtig ist noch die physikalische Bedeutung des Integrals. Bezeichnen wir mit | 
dP, und dP, die Komponenten der duBeren Kraft, die lings des Randstiickes d8 angreift, so gilt | 


ds=dx+jdy, —dP,=tdx—o,dy, — dP, =o, dx —tdy. (9) 


* Auszug aus der gleichnamigen Dissertation des Verfassers, Miinchen 1953. Erster Berichterstatter. 
Prof. Dr. G, Hettner, zweiter Berichterstatter Prof. Dr. L. Féppl. 


1 L.Sobrero, Theorie der ebenen Elastizitat, Hamburger Mathematische Einzelschriften, Heft 17, Leipzig 1934, 
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Wenn man andererseits mit f(z) die hyperkomplexe Funktion bezeichnet, welche die elastische 
Beanspruchung des Punktesz darstellt und dz = dx +j dy das entsprechende hyperkomplexe 
Differential nennt, so hat man 

f(z) dz = (0, + jf b—po, + ft) (dx + j dy) = ey, dx —t dy) stoner J° de —a, dy) 

ON Gs ae er 1 dP.) , (10) 
wobei die Punkte einige Glieder ersetzen, die fiir unsere Betrachtungen nicht von Interesse 
sind, da sie die Gréfe b enthalten, der keine physikalische Bedeutung zukommt. Geht 
man von infinitesimalen zu den endlichen GréBen iiber, so sieht man hieraus: Das erste 


| B 
und das letzte Element des Integrals f f(z) dz stellen die Komponenten — P, und. — P, der 
A 


auBeren Kraft dar. Wenn die Integrationskurve geschlossen ist und in das von ihr berandete 
Gebiet keine Singularitat der Funktion f(z) fallt, so ist das Integral gleich Null. Andererseits ist 
schon jetzt einzusehen, daB Singularitaten fiir die Funktion f(z) nur dort existieren kénnen, wo 
‘Angriffspunkte der auBeren Krafte liegen. Man kann also schlieBen, daB der Satz von Cauchy 
folgendem physikalischen Sachverhalt entspricht: Wenn in irgendwelchem Bereiche eines ebenen 
Systems keine auBere Kraft angreift, so stehen die am Rande dieses Bereiches angreifenden Span- 
nungen im Gleichgewicht. Ist die Integrationskurve geschlossen und fallt in das von ihr um- 
‘schlossene Gebiet ein Angriffspunkt einer auBeren Kraft (und deshalb eine Singularitat fiir die 


B 
Funktion f(z)), so verschwindet das Integral f{ f(z) dz nicht mehr, und seine auBersten Elemente 
stellen eben die Komponenten — P,,und — P,, der angreifenden Kraft dar. Damit ist der Residuum- 
-begriff physikalisch gedeutet. Auch das Integral 
B B 
Bie et SC a | et) AE — 90) aE 2 (11) 


hat eine wichtige physikalische Bedeutung: z) = x9 ++ J Yo ist hier ein beliebiger Punkt des ebenen 
Systems. Es ist namlich das erste Element des Integrals M(z)) gleich dem Moment der Spannungen 
beziiglich des Punktes zp. 


2. Die hyperkomplexe Spannungsfunktion. Sobrero beschrankt sich in seiner Arbeit auf die- 
jenigen Spannungszustande, bei denen der Rand lastfrei ist. Die Belastungen greifen also stets 
im Innern der Bereiche an. In den Anwendungen ist jedoch der Fall, daB die Rander belastet 
sind, der haufigere. Allerdings tritt hierbei eine fiir unseren Problemkreis charakteristische Schwie- 
rigkeit auf. Betrachten wir zum Beispiel die Normalspannungen langs eines Randes, der an einem 
Punkt durch eine Einzelkraft belastet ist, so sind diese tiberall Null, wenn die Kraft die Richtung 
der Normalen hat, ausgenommen an ihrem Angriffspunkt, wo sie unendlich werden. Soll nun die 
Funktion f(z) diesen Spannungszustand beschreiben, so muB sie auf dem Rand ihres Regularitats- 


b 
gebietes unstetig sein. Der Cauchysche Integralsatz besagt nun, daB / f(z) dz unabhangig von 


dem die beiden Punkte verbindenden Integrationsweg C ist, wenn C ganz im Regularitatsgebiet 
von f(z) liegt, aber auch dann noch, wenn C mit dem Rand des Gebietes zusammenfallt, voraus- 
gesetzt, daB f(z) auf dem Rand stetig ist. In dem genannten Beispiel der Einzelkraft wird die 
Giiltigkeit dieses Satzes auf dem Rande durchbrochen. Hierin liegt die mathematische Wurzel 
fiir die Einfiihrung der Spannungsfunktion D(z), welche die Unstetigkeit am Rande hebt. 

Wir beweisen den folgenden Satz: Das erste Element einer hyperkomplexen Funktion P(z) 
stellt die Airysche Funktion F(x, y) eines ebenen Spannungszustandes dar. Die Spannungen 
erhalt man durch die zweite Ableitung der Funktion ®(z). Denn es ist 


Oi) =F +jb+je+jed 
und also LQ ae 


‘Es gelten nun die Regularitatsbedingungen 


danach ist C, = dy, cy = 6, —2d,, Fi, = 6, , F, =— d,; 
Wir erhalten damit C= dad = id = = FF, 


PD : : 5 ae 
FEE ee ety Per oad fe a he 1g = Oy “reg Pe Pr ed es 
ufolge der bekannten Beziehung fiir die Airysche Funktion 


Og = F yy Oy = Fy» 1 — =F, ,. 
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Wir fassen kurz zusammen: @(z) = F(x; y) +j.--; (1 | 
2D : : : 
(0) = ya = Op JF eet I tas (15 


Wir beniitzen im Folgenden das Symbol (c) fiir die Gesamtheit der Spannungskomponenter, 
und nennen den Ausdruck ,,hyperkomplexe Spannung‘. Man beachte auch die Analogie mit dey 
komplexen Behandlung ebener Potentialstrémungen, wo der Realteil des komplexen Potentials d 
Strémungspotential angibt. Die erste Ableitung liefert dort die konjugierte Geschwindigkeit. Bei unp 
tritt an Stelle der ersten Ableitung die zweite, da die Differentialgleichung von vierter Ordnung ist! 

Es miissen noch die Randbedingungen fiir die Funktion ® aufgestellt werden. Dabei kniipfe 
wir an Gleichung (9) und (10) an: 


(0) de = a(S) ——aP, +j--- —7°dP, = —py de +j-- —P Pads, | 
wenn p,(s) und p,(s) die Lastdichte als Funktion der Bogenlange der Randkurve (C) angeben: | 
d : ; a aes | 
(22 = — f p(s) ds —C, +j---—3?( J pls) ds + G)). (144 
dz] Rand 0 x 0 | 
(C) (C) 
Wir setzen zur Abkiirzung ‘ ; | 
PRE (p(s) dss Gy Peo) Ooi ae (15) 
i 
oer (C) 
= ue 7 > 
ze aaa taal ee a 
Der Zusammenhang zwischen (16) und der Airyschen Funktion ist leicht herzustellen: 
dD oF 5 OF 
Fae Po Ne aye (17 
OF oF 
— =P. — aly Bes 18 
(ae Ws 4 & Me ( I} 
Wir kénnen nun die beiden Komponenten P, und P, in (16) noch zu emem vom Koordinatenti 
system unabhangigen Ausdruck vereinigen: 
d ‘ : : 
(F)pngl@®® (0) 008 (Hy) = (BX), j++ =— P, 008 (na) + P, cos (ny) bj 
Rand >. oF { is H 
=a eee ( | 


Dabei haben wir die Gré®en P, und P, zu dem Vektor S$ zusammengefaSt. Das Zeichen | bey 
deutet, daB die bezeichnete Vektorkomponente zu einer Koordinatenachse gehért, die senkrech [ 
auf der (x, y)-Ebene steht. Um die Vektorschreibweise zu legitimieren, werden wir spater desi 
Nachweis fiir den Vektorcharakter der GréBen P, und P, fiihren. SchlieBlich errechnet sich dei 
Wert von ® auf dem Rande zu 


ID rong = (— Py + j---—j® P,) de Rona =P, dy — P, dx +j---=dF+).... | 
Die Punkte deuten die restlichen Elemente an, in denen P,, und P, zufolge der Produktbildung 
mit Gliedern gekoppelt sind, denen eine unmittelbare physikalische Bedeutung nicht zukommii 


Diese Elemente werden bis auf unwesentliche Konstante durch reine Quadraturen aus den Regu4 
laritatsbedingungen erhalten. 


Die F 1 _ : I 

ie Forme D Rand =i GA RM ea is GS oa aa: (204 
kann ebenfalls in eine vom Koordinatensystem unabhangige Form gebracht werden: j 
Prods SD 1) dS Cs yo ed tee (21) 


Der erwahnte Nachweis des Vektorcharakters von 5 kniipft an Gleichung (9) an | 

dP, = p, ds = 0, dy —t dx = [o, cos (nx) +t cos (n y)] ds, 

dP, = p,ds =t dy —o,dx = [rt cos (n x) + a, cos (n y)] ds, 29" 

Px = 0, ¢08 (n x) +7 cos (ny) , ( ; 

Py =T cos (n x) + 0, cos (ny). 

Die rechten Seiten der Gleichung (22) zeigen Vektorcharakter, sie stellen die Multiplikation de 

Spannungstensors ¢ mit dem Einheitsvektor der Normalen dar. Damit ist die Schreibweise der 
Gleichung (19) und (21) gerechtfertigt. 

Zu den in den Randbedingungen auftretenden Konstanten sei noch bemerkt, da® sie bei ein 

fach zusammenhingendem Rand Null gesetzt werden kénnen. Der analytische Charakter de 
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Airyschen Spannungsfunktion wird sehr schén aus der Tatsache sichtbar, daB die ihr zugeordnete 
hyperkomplexe Spannungsfunktion im gesamten Bereich (Gebiet und Rand) regular ist. 


3. Die Zerlegungsformel fiir hyperkomplexe Funktionen. Jede hyperkomplexe Zahl 
sau tjy titu tj 
z=a+of, (23) 


wo « und # gewohnliche komplexe Zahlen bedeuten. Diese Darstellung verdankt man Levi-Civita. 
Dabei wollen wir von der Bedeutung der komplexen Einheit als RichtgréSe in der Caufschen 
Zahlenebene absehen und nichts weiter darunter verstehen, als die Wurzeln aus —1l. In der Tat 


laBt sich in die Form bringen 


rechnet man leicht nach, daB =n eee oi 
2 
ist, also Ss a a ee rae ‘ 
2 4 4 = ; 
Die GréBe 1 + 7? definiert man als hyperkomplexe Nebeneinheit : 
l+j?=o. (24) 
_Zufolge der hyperkomplexen Grundrelation verschwindet deren Quadrat 
oe=(1+pPr=—1+4+277?+7=0. (25) 


Den Beweis fiir die Méglichkeit der Darstellung z =« + wf, den man durch Ausrechnen und K oeffi- 

_zientenvergleich leicht fiihrt, zeigen wir nur fiir den wichtigsten Fall der speziellen Zahl z = x +7 y: 
; s74 78 34473 

z=atop=(«+iy)+o (—)=* +- oe 


> 


2 
73 

| Q+t=<+sjy—G+PL+fy a2 4+ iy. 
Augenscheinlich bezweckt diese Darstellung eine Vereinfachung der Rechnung; denn statt 

der komplizierten Rechenoperation mit hyperkomplexen Einheiten (z. B. Division) rechnet man 
bequem und in bekannter Weise mit komplexen Zahlen und mit einer hyperkomplexen Neben- 
einheit, deren Quadrat Null ist, die sich also wie eine kleine GréBe verhalt.1 Dieser Vorteil zeigt 
sich vor allem, wenn wir eine Funktion einer hyperkomplexen Veranderlichen betrachten. Die 
Levi-Civitasche Darstellung der hyperkomplexen Zahl laBt dabei wichtige Verallgemeinerungen zu. 
@(z) sei im Gebiet G, welches den Nullpunkt enthalten soll, regular und im Nullpunkt selbsteven- 
 tuell auch singular; dann Ja4Bt sie sich auf eine und nur eine Weise in eine Laurent reiheentwickeln: 


D(z) 7, Ne oe S a, (« + wf)’ = > a, (a” + @ By ar—t + (a? (wf)? ++) 


= > aa’top D> arvo’—1=D(a)+oP O(a). 
Diese Zerlegungsformel besagt: Jede analytische Funktion ®(z) der hyperkomplexen Veranderlichen 
z=x+yjy+ju+yv—a +o 1aBt sich zerlegen in 


P(e) = O[s) + wf Oe). (26) 
_Dabei bedeutet D(a) dieselbe Funktion ®, jedoch mit der komplexen Veranderlichen 
& = (e—u) +i(y—v) 27) 
als Argument, w die hyperkomplexe Nebeneinheit 1 + j? und f die komplexe Zahl 
faut. (28) 


Ebenso beweist man durch Reihenentwicklung: Die Ableitung der analytischen Funktion 


D(z) = D(x) + w 6 O'(x) ist gegeben durch 


2 
MHS +o. (29) 
Das unbestimmte Integral von P(z) lautet 

{ D(z) dz = f D(a) dx + wB P(x) + C. (30) 
Man kann auszweikomplexen Funktionen ®(a) und (a) die allgemeinste hyperkomplexe Funktion 
D(z) + w Ye) (31) 

1 Haufig gebraucht werden 

io=jo=j+j® und ij=—5t+yP=—14+9 : 


Sy W. Engl: Weiterfiihrung der Sobreroschen Theorie des ebenen Spannungszustandes Ingenieur-Archiy /_ 


erzeugen. Der Beweis erfolgt durch Zerlegung in die Elemente, welche die Form der allgemeinsten |, 
Bipotentialfunktion haben. Die in diesem Abschnitt abgeleiteten Formeln ergeben den Zusammen- 
hang zwischen den Funktionen einer hyperkomplexen Verdnderlichen und den gewéhnlichen | 
komplexen Funktionen. 


4, Polynome zur Lésung verschiedener Belastungsfalle des geraden Balken. Bei der Anwendung || 
der in Ziff. 2 aufgestellten hyperkomplexen Spannungsfunktion geht man zundchst von einfachen |/ 
Funktionen aus und sieht nach, welcher Spannungszustand von ihnen dargestellt wird. Dabei | 
erhalt man das einfache Ergebnis, daB die bekannten Lésungen! fiir Belastungsfalle des geraden || 
Balkens durch Polynome dargestellt werden. C bedeutet jeweils eine hyperkomplexe Konstante, | 
welche aus den Randbedingungen bestimmt wird. 


Reiner Zug: IMT Oe ole (32) ) 


reine Biegung: D sa (Oe5) (33) )| 
4 i 
Last am eingespannten Balken: Oz S Gz, (34) | 
A=2 
beiderseits frei aufliegender, gleichmaBig belasteter Balken: 


5 
ey Cras (35) | 
A=2 
beiderseits frei aufliegender Balken mit Dreiecksbelastung: 


6 
b=>G,2, (36) 


A=2 


beiderseits frei aufliegender Balken mit parabelférmiger Last- | 
verteilung: | 


4p 
‘ Py et (x — 1), 
# al 
3 3 7 
Abb. 2. Beiderseits frei aufliegender Balken mit —- A 
parabelférmiger Lastverteilung. ®@ eas C, axe 
=2 


(37) | 


Nachdem der letzte Fall bisher noch nicht gelést wurde, soll er ausfiihrlich beschrieben werden. 
Wir stellen zunachst die Randbedingungen fiir ® auf. Den Nullpunkt der Randkurve legen 


wir an die Stelle x = 0, y = —/A und durchlaufen von hier aus dem Rand im positiven Sinne. || 
Die erste Teilstrecke ist lastfrei, demnach gilt?: | 
(zr) Bk (0<x«<) (38) || 

z=x—jh Ey 5c | 

D(x—jh) =0 +) ei (39) | 


es folgt die Gerade x = 1; —hSySh. 

Bekanntlich haben wir keine strenge Lésung in dem Sinne zu erwarten, daB die Normal- 
spannungen fiir x = 0 und x =! verschwinden, sondern nur eine Lésung im Sinne des de Saint-Venant- 
schen Prinzipes. Danach ergeben die Normalspannungen keine resultierende Normalkraft und 
kein resultierendes Moment; das Integral iiber die Schubspannungsverteilung wird dabei gleich | 
der Auflagekraft. 

Das hei®t nun aber, daB wir langs der Geraden nichts vorschreiben, sondern nur an deren | 
Endpunkt x =/, y =h die fiir das Fortschreiten in mathematischem Umlauf erforderlichen | 
Konstanten errechnen miissen. 

Am belasteten Rand lauten die Randbedingungen 


Cian aoe (y = h) se || 
Oe -1) alg eae ae a) 


Fir das letzte Stiick des Randes x = 0, h < y < —h brauchen wir nichts vorzuschreiben. 
Den Ansatz (37) fiir ® miissen wir noch durch eine Jineare Funktion in z und eine Konstante 


ergdnzen, da wir die Konstanten in den Randbedingungen des lastfreien Randes gleich Null gesetzt 
haben: 


7 
D= > Cz, (42) 
ee oe A=0 
1 L. Féppl, Drang und Zwang, Bd. 3, Miinchen 1947, | 
® Siehe (19) und (21) und beachte den Faktorj auf der linken Seite von (19). 
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wobei die C, von der hyperkomplexen Form 


d C, a Cy aaa Cre eae Cys ey as C4 (43) 
sind. 
Wir haben also zunichst die erschreckende Zahl von 8- 4 = 32 Konstanten. Nun kann man 
b 
i? Cye = Css = Cg4 Cre Crs Coe = 0 (44) 


setzen, da sie zum ersten Element von @ keinen Beitrag liefern. Weiter iiberlegt man sich, daB 
das erste Element von @ keine héhere Potenz in x als die vierte enthalten darf, da auch in den 
Randbedingungen fiir ® keine héhere vorkommt. Zufolge dieser Bedingung erhalt man 


Cy, = Cy) = Cy = Gy = Gg = Cy, = 0. (45) 
Somit haben wir die Zahl der Konstanten auf 20 reduziert. Wir zihlen jetzt die Anzahl der aus 
den Randbedingungen folgenden Gleichungen ab: Schreiben wir mit unserem Ansatz die Glei- 
chung (41) an und fiihren von x‘ bis x° einen Koeffizientenvergleich durch, so erhalt man 5 Glei- 
chungen. Ebenso verfahrt man mit Gleichung (38), (39) und (40); das ergibt zusammen 20 Glei- 
chungen. Die Auflésung dieser 20 Gleichungen 1a48t sich nun in endlicher Zeit bewaltigen, und 
man erhalt die hyperkomplexe Spannungsfunktion 


2 $2 gt DL lis a ae an 
D = Fn laan? ont mito! ool” Ce TY 


h21 h21 
Att | yA 
ra |= 
Ue ae ein nh PE AT 
—( +4500 Td 3) St aed a4} Se. 
Die hyperkomplexe Spannung errechnet sich zu 
Ere) j y 1 HO oof SS | 1bsT Ps) -g dh? 
(0) = “ga = pin” 7 Hip tells + [ae i a as ieh 2 |# 


227k! |. ht 2h . fh? ‘hel 
- aD Ds oe aa 3) a ee es (= —7} (47) 


Durch Zerlegen in ihre vier Elemente erhalt man 


opal pt Bayt Bley — Sy 2B ys ty $18 y OF ey 
See 
OBE HA tla Taw ay A 
—3W ley +c y +2 hea? — 2h 121, ra 
3 a eee 
ee eee 
Der Vollstandigkeit halber geben wir noch die Airysche Spannungsfunktion an: 
_ esa ree 


5. Flamants Lésung fiir die unendliche Halbebene, deren Rand durch eine Einzelkraft belastet ist. 
Zu Anfang miissen wir darlegen, was wir unter einer Einzelkraft verstehen wollen, da in unseren 
bisherigen Betrachtungen iiber die Randbedingungen dieser Begriff noch nicht aufgetreten 
ist. Dabei kniipfen wir bei (15) an und stellen uns zum Beispiel eine konstante Normalbela- 
stung der Intensitat p langs der Strecke von —¢ bis +-¢ vor. Zur Einzellast kommen wir durch den 


Grenziibergang p,—> co und ¢ -> 0 derart, daB ii py dx endlich bleibt und gleich dem Betrag der 


Einzelkraft P wird. In den Randbedingungen (19) und (21) stellen dann P,, und P, keine Abkiir- 
zungen fiir Funktionen langs des Randes mehr dar, sondern sind direkt die Betrage der Kompo- 
nenten der Einzelkraft, wenn wir die Konstanten C, und C, aus (15) gleich Null gesetzt haben. 
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Nach dieser Vorbemerkung kénnen wir die Randbedingungen formulieren: 


(ey es ty, [Of 0 ey. 8 (—o<x<0), (50) 
(ey = FQ, WL, =— Pati... OS¥<ex). 1) 
% 


Abb, 3. Grenziibergang zur Einzellast. Abb. 4. Die unendliche Halbebene mit Einzellast. 


Zum Aufsuchen der Lisung gehen wir zweckmafig von der zweiten Ableitung der Funktion ®, | 
also der hyperkomplexen Spannung (a) aus. Zufolge des bei der Definition der Einzellast gemachten | 
Grenziiberganges wissen wir, da die Spannungen im Nullpunkt unendlich werden. Wir denken | 


und (c) in eine Laurent-Reihe um den Nullpunkt entwickelt: 


()= D 4,2’. (52) 
Da die Spannungen im Unendlichen verschwinden miissen, fallen alle Glieder mit positivem Expo- 
nenten fort und es bleibt 22D oo l 
p= 0) = Sas. (53) 


Die Reihe konvergiert auBerhalb von z = 0 gleichmafig, wir diirfen also gliedweise integrieren und 


erhalten gg co 1 
> ee Inz + >i a_, (—») a! +C, (54) 


y=2 


as 1 
@ =a_,2(mz—1 a a=, 
a_,2(Inz )+a_,Inz ae i) a 


Die Lésung folgt nun unmittelbar aus der Tatsache, daB die Spannungsfunktion ® am Rande stetig ist. 
Sie muB daherim Nullpunkt existieren und nimmt dort den Wert Null an. Demnach miissen alle a_, 


mit y = 2 und die Konstante D verschwinden. Hier sehen wir die Bedeutung der hyperkomplexen || 
Spannungsfunktion fiir die Behandlung von Aufgaben mit belasteten Randern. Wir erhalten somit || 


® =a_,2(mz—1)+ Cz, Pains +C. (56) 
Die Konstanten bestimmt man aus den Randbedingungen. Man beachte, daB 
[In-z],—. = {Ina + oB | 
p=0 
(hava yea ae [anit], =e Nan |e ae 


Dabei ist der Verzweigungsschnitt des Logarithmus langs der positiv reellen Achse gefiihrt. Weiter 
kénnen wir bei | 


C=Q4iG+tPQtPC, 
die Konstanten C, = C; = 0 setzen [siehe (44)]. Wir erhalten 


a,==(-jP47Q), C=—(P+/*Q). 


Die in z linearen Glieder, welche durch das Nullsetzen der Konstanten auf der lastfreien, negativ | 


reellen Achse hereinkommen, kénnen, da sie fiir die Spannungen ohne Belang sind, auch wieder | 
weggelassen werden. Die hyperkomplexe Spannungsfunktion lautet dann 


We A : 
== (—j P+ j*Q)ehn: (57) 
und damit die hyperkomplexe Spannung 


2D 2 ; p 
(0) = “Ge = IP +P QRS. (58) 


| 


Te aedsy (55) | 


| 


| 


q 


1 
' 
I 
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Durch Anwendung der Zerlegungsformel erhalt man leicht die bekannte Airysche Spannungs- 
funktion und die reellen Spannungskomponenten 


Pa(T2—Sy)arotg=. (59) 


4 7 


Wahrend sonst in der Literatur! iiblicherweise die Airysche Spannungsfunktion als gegeben an die 
Spitze gestellt wird und daraus die Spannungen hergeleitet werden, wird hier ein einfacher Weg 
aufgezeigt, wie man diese selbst findet. Dabei wird lediglich von der Tatsache Gebrauch gemacht, 
dafs die hyperkomplexe Spannungsfunktion eine im ganzen Bereich regulare Funktion ist. 


6. Einfiihrung krummliniger Koordinaten in die Spannungsfunktion und Zerlegung des Span- 
nungstensors in Hauptspannungen und krummlinige Komponenten. Bei den bisherigen Betrach- 
tungen tiber die hyperkomplexe Spannungsfunktion @(z) war stets z = x +j y, so da® auch die 
Airysche Funktion F als Funktion von x, y erschien. Es ist aber auch erwiinscht, F als Funktion 
pe uamliniger Koordinaten &, 7 zu betrachten. Man kénnte natiirlich einfach 
e=sin), y=MEn) mit F720 
in F(x, y) einsetzen und hatte dann F(g (&,7), h(E, )) = F*(é, 7). 
_ Fur den am meisten interessierenden Fall orthogonaler, isometrischer Koordinaten erméglicht 
die hyperkomplexe Zerlegungsformel einen praktischeren Weg. Wir gehen zuriick auf (26) und (27): 
Orthogonale isometrische Koordinaten werden durch Betrachtung der konformen Abbildung einer 
y-Ebene mit y = € + 17 auf die a-Ebene erhalten, bei der die Geraden § = konst. und 7 = konst. 


in der a-Ebene die gewiinschten orthogonalen Kurvenscharen ergeben: 
| a= py), x+iy=pE+in)=un) +ivgn, B=—-F=—F. (60) 


Setzt man (60) in ®(z) ein, so erhalt man 


| D(z) = D(ply)) — 23° By). (61) 
Schreibt man 
: P(p(y)) = Fy), (62) 
mmt iv ’ 
a : De) = Yo) ie Gy = Plumb OE M)) bie (63) 
a 


Wir wollen noch darauf hinweisen, da zwischen der Einfiihrung orthogonaler isometrischer 
Koordinaten durch a = g(y) bei einem komplexen Potential und bei einem hyperkomplexen Bi- 
= ein wesentlicher Unterschied besteht. Betrachtet man namlich ®(«), so kann man dafir 


is (0) =O (x + iy) = Poly) = Vy) = PE +in)s 
das heiSt, durch Einfiihrung von « = g(y) in die analytische Funktion ®(x) entsteht wieder eine 
analytische Funktion Y(y) = Y ( + in). Dagegen laBt sich die aus D(z) = P(x + jy) durch 
Einfiihrung von «a = ¢(y) entstehende Funktion nicht als Funktion von € + j 7 schreiben; sie ist 
also keine im hyperkomplexen Sinne analytische Funktion von é, 7. 
Der hyperkomplexe Spannungstensor hat die Komponenten 
d2 . 6 : 

(0) = FS Ho, + jb—jo, + j*tey. (13) 
Wir betrachten im folgenden den Ausdruck (1 + @/2) (c) und kénnen ihn wie jede hyperkomplexe 
Zahl in zwei komplexe Zahlen und die Nebeneinheit w zerlegen: 


(1+ 3) =, +4.) +i6—t.5) Poe ue (64) 


Der Realteil der ersten komplexen Zahl ist die Spannungssumme, fiir die wir wegen der Invarianz 
der Tensorspur auch schreiben kénnen 


G, +0, = 0, + 03. (65) 
Die zweite, mit w verknipfte, komplexe Zahl schreiben wir, wie folgt: 
_ tare tg sey 
Oy — Ox ite = Ve, oye At, ¢ Oy — Fg (66) 


Der Betrag gibt gerade die Hauptschubspannung an: 


GSS ee = 
T= 5 = VO, —9,)2 + 4 ty, 


1 Siehe FuBnote 1 von S. 132. 
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und der Winkel die Richtung der Hauptspannungen mit der x-Achse 
2Txy 


tg2Q,,2= 


Oy — Ox é 
Wir kénnen also zusammenfassen: 


(1+ 5)@ = +o) +i Sg BG 2 (6) 


Die Transformationsgleichungen der Spannungen 


Gy Ti Gx Oy — Ox s 
ie ee cos2@ +17,,8in20, 


Cb SA SEE 
2 


a ‘ 
= cos2O —t,,sin20, 


Ov — O. 7 
Tr, = 7 sin 2 O +7,, cos 20 


lassen sich zusammenfassen zu o; +0, =o, +0,; (6€ 


ms its, = “x * ¢08 2 O —Tezy sin 20 od ard ** sin2O +t, cos 2 6) 


=(2 52 +in,,)e°. (4 
Es seien nun krummlinige isometrische Koordinaten wie in (60) gegeben, dann gilt 
Ox . OY 
Ox oy OE ia OE @ 
OS aS oa ene = |p'(y)| (cos O +i sin O). 
(i) +(@) 
Wir haben somit die Beziehung ; 
bie = (7 
lp’y)/? 
Multipliziert man die Gleichung (64) mit (70), so erhalt man 
“(y)? oO i On — OE : 
(1 f- 3) 0) On a + 0 (5 Tas) e? Sar en (a +- ite) : 
ZusammengefaBt gilt 
. On — 4, : aa 
(c) =o; +0, +i---+@-:-- ifcas Z)IperrO) = Datel ees + itg,). (7 


Mit Hilfe von Gleichung (71) ist nun der Ubergang auf krummlinige Tensorkomponenten leich) 
méglich, wobei noch beliebig ist, ob die krummlinigen Komponenten Funktionen von x, y ode} 
nach Gleichung (63), die natiirlich nicht nur fiir die hyperkomplexe Spannungsfunktion, sondersj 
fiir jede hyperkomplexe Funktion gilt, Funktionen der krummlinigen Koordinaten &, 7 sind. 

7. Zusammenfassung. Wir gehen von der Sobreroschen Theorie der Beschreibung des ebenex 
Spannungszustandes mittels hyperkomplexer Funktionen aus. Wahrend Sobrero sich auf die Ani 
wendungsfalle beschrankt, bei denen Krafte nur im Innern des Materials angreifen und die Randev 
durchweg lastfrei sind, erweitern wir diese Theorie’ durch Kinfiihrung der hyperkomplexen Span) 
nungsfunktion (12), (13) und behandeln damit“ auch. belastete. Rander.'| Die hyperkomplexe Span 
nungsfunktion hebt die von den Belastungen herriihrenden Unstetigkeiten am Rande auf und stell: 
somit die Giiltigkeit des Cauchyschen Integralsatzes im ganzen Bereich wieder her.' Ferner ist fi 
das praktische Rechnen mit hyperkomplexen F tnhronee die Zerlegungsformel (26), (27) und (28} 
von Bedeutung. . Sie deckt den Zusammenhang der hyperkomplexen mit den komplexen Funk: 
tionen auf. In den Anwendungen wird gezeigt, das man die iiblichen Belastungsfalle des geraderi 
Balken mit Hilfe von Polynomen lésen kann: (32) bis (37). Dabei ist es méglich, weitere bis jetzt! 
unbekannte Faille in das aufgestellte Schema einzuordnen und zu lésen. In dem Falle des Balken 
mit parabelférmiger Belastung ist die Behandlung bis zu Ende durchgefiihrt: (49). Fir die unend4| 
liche Halbebene mit Einzellast wird eine Ableitung angegeben, welche nur die Stetigkeit der Span: 
nungsfunktion und das Verschwinden der Spannungen im Unendlichen benutzt: (57). Die Zer? 
legungsformel fiir die hyperkomplexe Funktion erméglicht in einfacher Weise die Einfihrung 
orthogonaler, isometrischer Koordinaten in die Spannungsfunktion (63) und die Zerlegung i 
krummlinige Tensorkomponenten: (71). 


(Eingegangen am 21. Juli 1958.) 
Anschrift des Verfassers: Dr. W. Engl, Karlsruhe, Germersheimer Str. 20. 
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